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Summary 

This study examines the use of embedded shape memory alloy (SMA) actuators for adaptive 
control of the themomechanical response of composite structures. Control of static and dynamic 
responses are demonstrated including thermal buckling, thermal post-buckling, vibration, sonic 
fatigue, and acoustic transmission. A thermomechanical model is presented for analyzing such 
shape memory alloy hybrid composite (SMAHC) structures exposed to thermal and mechanical 
loads. Also presented are ( 1 ) fabrication procedures for SMAHC specimens, (2) characterization 
of the constituent materials for model quantification, (3) development of the test apparatus for 
conducting static and dynamic experiments on specimens with and without SMA, (4) discussion 
of the experimental results, and (5) validation of the analytical and numerical tools developed 
in the study. 

The constitutive model developed to describe the mechanics of a SMAHC lamina captures 
the material nonlinearity with temperature of the SMA and matrix material if necessary. It is in 
a form that is amenable to commercial finite element (FE) code implementation. The model is 
valid for constrained, restrained, or free recovery configurations with appropriate measurements 
of fundamental engineering properties. This constitutive model is used along with classical 
lamination theory and the FE method to formulate the equations of motion for panel-type 
structures subjected to steady-state thermal and dynamic mechanical loads. Mechanical loads that 
are considered include acoustic pressure, inertial (base acceleration), and concentrated forces. 
Four solution types are developed from the governing equations including thermal buckling, 
thermal post-buckling, dynamic response, and acoustic transmission/radiation. These solution 
procedures are compared with closed-form and/or other known solutions to benchmark the 
numerical tools developed in this study. 

Practical solutions for overcoming fabrication issues and obtaining repeatable specimens are 
demonstrated. Results from characterization of the SMA constituent are highlighted with regard 
to their impact on thermomechanical modeling. Results from static and dynamic tests on a 
SMAHC beam specimen are presented, which demonstrate the enormous control authority of the 
SMA actuators. Excellent agreement is achieved between the predicted and measured responses 
including thermal buckling, thermal post-buckling, and dynamic response due to inertial loading. 

The validated model and thermomechanical analysis tools are used to demonstrate a variety of 
static and dynamic response behaviors associated with SMAHC structures. Topics of discussion 
include the fundamental mechanics of SMAHC structures, control of static (thermal buckling and 
post-buckling) and dynamic responses (vibration, sonic fatigue, and acoustic transmission), and 
SMAHC design considerations for these applications. The dynamic response performance of a 
SMAHC panel specimen is compared to conventional response abatement approaches. SMAHCs 
are shown to have significant advantages for vibration, sonic fatigue, and noise control. 
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Chapter 1 
Introduction 


1.1 Background and Motivation 

The effects of noise and vibration transmission and propagation proliferate throughout our 
everyday experiences. These effects range from annoyance or interference with other activities 
to mechanical failures, loss of hearing, and even loss of life. Attempts to attenuate or otherwise 
remedy noise and vibration problems date back as far as our acknowledgment of them. Methods 
of noise/vibration reduction typically involve increases in one or more of three things, mass, 
stiffness, or damping. Conventional noise/vibration abatement attempts typically consist of 
post-design treatments to reduce the noise/ vibration at the source or reduce the perception of 
the noise/vibration. Examples of conventional treatments include sound-absorbing blanket or 
foam, lead-vinyl or viscoelastic damping layers, viscoelastic vibration isolators, tuned vibration 
absorbers, etc. These add-on treatments can be adequate for some applications where weight is 
not an issue and system parameters do not vary widely. However, there are other applications 
where weight-efficiency and adaptability are all-important. One class of such applications arises 
in the area of high-performance aerospace structures. 

Excitation levels for many aerospace vehicle structures are high due to engine noise and 
turbulent boundary layer fluctuating pressures. The ever-increasing need for weight-efficient 
structures in these applications leads to significant sonic fatigue and interior noise issues. Con- 
ventional aerospace structures typically employ passive treatments such as structural stiffening, 
constrained-layer damping, and acoustic absorption materials to reduce the structural response 
and interior acoustic levels. These treatments typically suffer from substantial weight penalty 
and are often limited to relatively low temperatures and/or high frequencies. Smart materials 
may enable new structural design paradigms for high-performance, durable structures. 

A significant amount of research in the last several years has greatly advanced the field of 
control of flexible structures through the use of smart materials. The field of smart structures 
involves a broad range of technologies that enable the realization of structural systems that are 
able to sense and to control their own behavior, such that the range of operational performance 
may be extended over conventional structures. Smart structures have the ability to change 
stiffness, shape, natural frequency, damping, and other mechanical characteristics in response 
to changes in temperature, electric field, or magnetic field. The most common smart materials 
are piezoelectric materials, shape memory alloys, magnetostrictive materials, electrorheological 
fluids and magnetorheological fluids. All of these adaptive materials have great potential in a 
variety of applications. Shape memory alloys (SMAs) appear to have significant advantages for 



high-performance aerospace vehicles in the areas of structural response control, structural shape 
control, and damping enhancement. 

Shape memory alloys have been investigated for a variety of applications since their 
discovery. Early potential for applications in self-erecting structures, thermally actuated devices, 
and energy -con version systems was identified by Buehler and Wang [1], Other potential 
applications including damping or energy absorbing devices, thermally actuated couplings and 
fasteners, and biomedical devices were briefly described by Wayman and Shimizu [2], Later, 
Wayman [3] and Otsuka and Shimizu [4] gave more detailed descriptions of specific devices such 
as pipe couplings, thermostats, a robot hand, and various biomedical systems. Other researchers 
have reported progress on similar applications, but an entire new field of applications was 
created when Rogers and Robertshaw [5] introduced the idea of embedding SMA actuators in 
a composite laminate. A structure of this type has been termed a shape memory alloy hybrid 
composite (SMAHC). 

Numerous researchers have investigated the use of SMA actuators, either external to a 
structure or embedded within the structure, for static and dynamic structural control. Studies 
have considered SMA actuators, external to the structure, for active vibration control [6, 7, 8] 
and shape control [9, 10] of cantilevered beams. Other studies have proposed SMA actuators for 
vibration control of space structures [11] and presented analytical formulations to demonstrate 
the vibration/structural acoustic control of SMAHC panels [12, 13]. Active vibration control 
of flexible linkage mechanisms using SMAHCs has been investigated by Venkatesh et al. [14]. 
Studies of thermal post-buckling [15] and random dynamic response [16] suppression of SMAHC 
panels were presented using a thermoelastic finite element (FE) formulation. A different approach 
was offered by Ro and Baz [17, 18, 19], where thermal, static, and dynamic analyses were 
developed for the case of SMA actuators passing through sleeves in composite plates. Birman 
[20] investigated optimal distributions of SMA actuators in sleeves to improve the buckling 
performance of sandwich panels under mechanical loading. Theses have also been devoted to 
studying the dynamic response tuning capability of Nitinol actuators [21] and predicting the 
structural acoustic behavior of SMAHC panels [22]. 

Shape memory alloys are a class of materials that exhibit a dififusionless martensitic transfor- 
mation when cooled from the higher-temperature austenitic state. The interaction of temperature 
and applied stress in driving the martensitic and reverse transformations can be used to exploit 
phenomena such as the shape memory effect (SME) [23] and pseudo-elasticity. The SME can be 
described in simple terms in the following manner. A SMA can be easily deformed in the low 
temperature martensitic condition and can be returned to its original configuration by heating 
through the reverse transformation temperature range. This type of SME is termed free recov- 
ery. Conversely, in a constrained recovery configuration, the SMA element is prevented from 
recovering the initial strain and a large tensile stress ( recovery stress) is induced. A situation in 
which the actuator performs work (deforms under load) is called restrained recovery. Extensive 
work has been done to characterize these materials, both qualitatively through theoretical models 
[24, 25, 26, 27] and quantitatively for particular alloy compositions [28, 29, 30]. Compilations 
of papers have also been published [31, 32, 33, 34], which give a good overview of the char- 
acteristics of SMAs and their applications. Birman [35] gave a comprehensive review of work 



done in the areas of alloy characterization, constitutive modeling, and applications. Some detail 
on thermomechanical characteristics of SMAs will be given in Chapters 2 and 3. 

Although there are several alloy systems that exhibit the shape memory effect, Nitinol 
has the most engineering significance. The name for this family of alloys was derived from 
Nickel, Titanium, and Naval Ordnance Laboratory and was patented by Buehler and Wiley [23]. 
The SME for this family of alloys is limited to nearly equiatomic compositions, i.e.. Nickel 
content from 53 to 57 percent by weight. Nitinol is considered to be the shape memory alloy 
with the most potential for engineering applications because of its ductility at low temperature, 
high degree of shape-recovery capability, large pseudoelastic hysteresis, corrosion and fatigue 
resistance, biomedical compatibility, and relatively high electrical resistance. 

The literature referred to in this introduction is by no means comprehensive, but should 
serve as a guide to work that has been done and the state of the art in this area. Other 
references to specific works in the literature appear at relevant places throughout the remainder 
of this document. The scope and objectives of the present work, including a description of the 
contributions made to the state of the art, are given in the next section. 


1.2 Scope and Objectives 

A Nitinol alloy was chosen for this study because it is the most readily available alloy 
and because of the favorable shape memory and electrical resistance characteristics. The 
relatively high electrical resistance is beneficial for experimental situations involving activation 
of the actuators by resistive heating. Because of the ability of SMA actuators to induce 
large tensile forces and recover large strains, they are particularly attractive for applications 
involving adaptive structural stiffening (constrained recovery) and shape control (restrained 
recovery). The environment associated with many aerospace applications includes inherently 
elevated temperatures and broadband surface pressure excitation. Embedded SMA elements are 
particularly suited to such applications because they generate global structural stiffening and 
obviate control electronics and auxiliary power. It is for these reasons that the present work 
will focus on SMAHC panel-type structures with embedded actuators operating in a constrained 
recovery mode for adaptive structural stiffening. 

The purpose of the present work is to develop and experimentally validate an analytical 
tool that captures the material nonlinearity of SMAs and complex thermomechanical behavior 
of SMAHC laminates. Constitutive theories for SMAHC structures could be developed from 
the previously mentioned models for SMA thermomechanical behavior [24, 25, 26, 27], as was 
done by Boyd and Lagoudas [36] and Birman [37]. However, these models are difficult to 
use in practice. An alternative approach is to employ a constitutive model which makes use 
of experimental measurements of fundamental engineering properties. Two such models are 
developed as part of the present work. Differences between the two models and reasons for 
development of the refined model are discussed in Chapter 2. The model is used in mechanics- 
of-materials approach to develop the constitutive relations for SMAHC material systems. The 
model captures the material nonlinearity of the SMA, as well as that for the matrix material 
when appropriate, and is shown to be quite general in that it can be used to model constrained, 
restrained, or free recovery behavior. 
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The complex nature of SMAHC structures and the goal of experimental validation dictate the 
use of the FE method to facilitate accurate thermomechanical analysis. Because no commercial 
FE codes have the capability of modeling the constitutive behavior of SMAHC laminates, a 
research FE code was developed as part of this work. However, it is a goal of this work to 
cast the constitutive model and thermomechanical FE formulation such that it can be readily 
implemented in a commercial FE code. 

The experimental configurations that will be considered for demonstration of adaptive 
structural stiffening by thermal activation of embedded SMA actuators involve immovable in- 
plane boundary conditions. This condition leads to a potential for structural instability due 
to thermal buckling and demands that the FE analysis be capable of predicting static, thermal 
buckling and post-buckling behavior in addition to analysis of the response due to practical 
dynamic excitation types. To this end, the constitutive relations are used in conjunction with von 
Kdrman strain-displacement relations to allow for the generality of large, out-of-plane deflection 
due to thermal post-buckling. 

A brief description of the work performed in this study and contributions made have been 
given in this section. The manner in which the detailed discussions of various topics are organized 
in this document is given in the following section. 

1.3 Research Organization 

Some detail of the general thermomechanical behavior of SMAs is given in Chapter 2 with 
an emphasis on phenomena observed in the mechanical testing performed in this work. That 
is followed by development of a constitutive model to describe the behavior of SMA actuators 
embedded in a composite matrix. The constitutive model is used with the von Karman strain- 
displacement relations and classical lamination theory to develop the resultant force and moment 
relations for a SMAHC laminate. The particular material system used in this study is described 
in Chapter 3. An overview of the processes involved in fabricating SMA hybrid composite 
specimens is followed by results from thermomechanical property testing of the constituent 
materials. 

Chapter 4 begins with a description of the rectangular plate element used in the finite element 
formulation. This is followed by development of the equations governing the thermomechanical 
response of a SMAHC laminate to steady-state thermal and dynamic mechanical loads. The 
governing equations are separated into static and dynamic components by assuming small (geo- 
metrically linear) dynamic motion about a potentially large (geometrically nonlinear) thermally 
deflected shape. 

Description of the solution procedures begins in Chapter 5 with the static response. Equations 
governing static stability and large static deflection are first developed and then used to present 
thermal buckling and thermal post-buckling solution procedures. A overview flowchart is shown 
in each case. The solution procedures are then validated against a series of known solutions. 
The chapter is concluded with experimental validation of the model for the case of a SMAHC 
beam clamped at both ends and heated uniformly by DC electrical current. 

Continuous and discrete frequency domain solutions are shown in Chapter 6 for predicting 
the dynamic response. The dynamic response solution procedure is validated against known 
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solutions and experimental measurements for isotropic and conventional composite beam spec- 
imens excited by random inertial loading. The analysis is then validated with experimental 
measurements of the dynamic response performance of a SMAHC beam specimen subjected 
to steady-state thermal and dynamic inertial loads. Chapter 6 concludes with correlation of 
predicted and measured dynamic responses of a clamped aluminum panel subjected to random 
acoustic pressure, which introduces the topic of acoustic radiation in preparation for Chapter 7. 

Continuous and discrete frequency domain solutions are shown in Chapter 7 for predicting 
the acoustic radiation from a known panel response. The acoustic radiation solution procedure is 
validated against a classical analysis method. Chapter 7 is concluded with comparisons against 
experimental measurement from a fully clamped aluminum panel subjected to spatially uniform, 
temporally random acoustic pressure. 

A variety of application and demonstrations are discussed in Chapter 8 which show the 
complex thermomechanical behavior of SMAHC structures and the effectiveness of embedded 
SMA actuators in attenuating dynamic response as compared to conventional response abatement 
approaches. Finally, a detailed summary and concluding remarks are given in Chapter 9. 
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Chapter 2 

Background on SMAs and Constitutive Modeling 


A general discussion of the thermomechanical behavior of SMAs begins this chapter. This 
discussion will serve as an overview and motivation of the SMAHC constitutive modeling 
section. Material property data for the material system used in this study will be presented in 
the next chapter, which will be used to quantify the constitutive model for the specific materials 
used. The constitutive model will be used in Chapter 4 to develop a finite element formulation 
for the thermomechanical response of SMAHC structures. 

2.1 Background on Shape Memory Alloys 

A very brief introduction to SMAs and the SME was given in Chapter 1 . An overview of 
some essential characteristics of SMAs will given in this section to facilitate understanding of 
results and discussions to be presented later in this document. 

The thermoelastic martensitic and reverse transformations are shown schematically in Figure 
2. 1 in terms of the martensitic fraction versus temperature. At high temperature, a shape memory 
alloy exists in a usually cubic phase called austenite. Martensite is the low-temperature phase 
and has a lesser symmetry. The reference temperatures indicated in Figure 2.1 mark the start and 
finish temperatures for the forward transformation to martensite and the reverse transformation 
to austenite, respectively. The thermoelastic martensitic transformation occurs when the alloy in 
the austenitic parent phase is cooled through the transformation temperature range M s ‘ -A/y * 
producing martensite. The reverse transformation occurs when the alloy is heated through the 
reverse transformation temperature range .A, — + -A j and the material structure returns to that of 
austenite. It is known that a shape memory alloy specimen shows no macroscopic deformation 
during the martensitic transformation, described by Wayman and Duerig in [32] for example. 
Because the martensite is of a lower-order symmetry than the parent phase, the crystalline 
structure change in the forward transformation produces several martensite variants that are self- 
accommodating (called twinning ) so that no macroscopic deformation results. However, the 
crystalline change to austenite is unique. When a stress is applied to the twinned martensite, 
the variant in the sense of the stress grows at the expense of other variants (called detwinning). 
Thus, inducing the crystallographically unique reverse transformation causes a strain recovery, 
and the material reverts to its original configuration. This sequence is the basis of the SME and 
is shown schematically in Figure 2.2. 

In Nitinol for example, strains as high as 8% can be completely recovered by heating the 
alloy above A f (Cross et al. [28]). Recall that this type of SME phenomenon is termed 
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Martensite (twinned) 



Austenite Martensite (deformed) 


Figure 2.2 Schematic of the shape memory effect, after Wayman and Duerig in [32], 


free recovery. This study is primarily concerned with constrained recovery, where the SMA 
element is prevented from recovering the initial strain and a large recovery stress is induced. 
The recovery stress is a function of initial strain and temperature, as shown in Figure 2.3 for 
example. The recovery stress will approach, but obviously not exceed the austenitic yield 
stress due to plastic slip or production of stress induced martensite, which will be discussed 
subsequently. The relationship between these two yield stresses will be discussed subsequently. 
For example. Cross et al. [28] found that the maximum yield stress for an alloy in a 100 mil 
diameter rod form was 90 ksi, while the corresponding maximum recovery stress was 80 ksi 
at 8% initial strain. The austenitic yield stress is approximately a factor of 10 greater than 
that of the corresponding martensitic phase. The Young’s modulus also increases by a factor 
of three to four under the same conditions (see Figure 2.4). Note that following a process of 
free recovery, no shape change is produced in the alloy when it is cooled again below 3/ y and 
the specimen must be strained again to repeat the SME. This type of behavior is termed one- 
way shape memory. Two-way shape memory, where the alloy “remembers” a low-temperature 
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Recovery Stress, MPa 



Temperature, °C 


Figure 2.3 Nitinol recovery stress versus temperature and initial strain, after Cross et al. [28]. 
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Figure 2.4 Nitinol Young’s modulus versus temperature, after Cross et al. [28]. 


martensitic configuration and a high-temperature austenitic configuration, can be imparted to 
a SMA with appropriate training, as described by Perkins and Hodgson in [32]. One-way or 
two-way shape memory behavior can be trained (rendered stable and predictable) into a SMA 
by a number of means. A description of one simple method follows. 

One-way training can be accomplished by imparting a specific amount of plastic (detwinning) 
deformation to the SMA in the low-temperature martensitic condition, heating the material to 
induce the austenite phase transformation and allowing the alloy to recover the initial strain, and 
cooling the material to induce the transformation to martensite. This procedure is repeated until 
the hysteresis becomes stable. Two-way shape memory training can be performed by restraining 
the SMA element, after being plastically deformed at low temperature, and subjecting it to a 
number of thermal cycles, alternating between a specific low and high temperature. The benefit 
of this training method is that the shape memory element “memorizes” two stable shapes (a 
low-temperature shape and a high-temperature shape) at the expense of some degradation in 
performance. 

There are other means by which a one-way actuator can be returned to a strained state 
following recovery. Two examples include operating actuators in agonist-antagonist pairs or 
opposing the actuators with a bias spring. In either case, the active actuators perform work 
(deform under load). Recall that operation of an actuator in this type of configuration is 
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termed restrained recovery. Obviously, the amount of force generated in restrained recovery 
is diminished relative to constrained recovery. 

The martensitic transformation can also be stress induced in a shape memory alloy at 
temperatures above Ms. An increase in the applied stress produces an effect analogous to 
a decrease in temperature. However, only one martensite variant ( stress induced martensite, 
SIM) grows in the direction most favorable for the applied stress. This phenomenon is 
exemplified in the stress-strain curves for an alloy at various temperatures shown in Figure 

2.5. For temperatures in the range A/, < T < .4„ SIM will be formed and retained after 
unloading. Partial reversion to austenite will occur upon stress release for temperatures in the 
range .4., < 7 < Aj. At temperatures approaching or greater than Af, the complete reverse 
transformation ensues immediately upon stress release as the martensite is completely unstable 
without stress assistance. These phenomena are all classified as pseudoelastic effects. The 
latter effect of full reverse transformation upon stress release is an example of a subset of 
pseudoelasticity called superelasticity. 

The pseudoelasticity discussed here is called transformation pseudoelasticity. Twinning 
pseudoelasticity exists as well, which gives rise to rubber-like behavior. This phenomenon 
occurs when a SMA is allowed to age in the martensitic phase and is then stressed to induce 
prestrain. Instead of maintaining the strain, the strain is recovered upon stress release (in a very 
similar manner to superelasticity, but entirely in the martensitic phase). This behavior indicates 
that if one is going to prestrain a SMA actuator, it must be done shortly after inducing the 
martensitic transformation on cooling, i.e., immediately after removing strain induced in the 
material during packaging. 

Additional insights can be afforded from the stress strain curves shown in Figure 2.5. When 
the material is in the purely martensitic phase, the yield stress to induce martensite detwinning 
and produce appreciable strains is quite small. If loading were continued as shown in Figure 

2.6, the martensite would eventually become completely detwinned and would deform elastically 
until plastic slip is induced. Obviously, the amount of reversible strain is limited by the onset 
of plastic slip in the martensitic phase. Rather complex behavior is involved at intermediate 
temperatures M / < T < Af depending upon the state of the alloy. The stress required to yield 
the material in the austenitic phase is much higher than that for martensite. At temperatures 
in the vicinity of .4 /, yielding is due to the production of SIM. The stress required to induce 
SIM continues to increase with increasing temperature above 4 j until a critical temperature is 
reached, at which point plastic slip will be induced in lieu of SIM. 

The transformation temperatures of an alloy can be measured by monitoring practically any 
property of the alloy, such as electrical resistance, but a particularly accurate and insightful 
means is differential scanning calorimetry (DSC). The transformation temperatures for Nitinol 
can be altered by changing the composition of the alloy. For example, Buehler and Wang [1] 
report that .4, for Nitinol can be varied from -50° C to 166°C through changes in the Nickel 
content from 54.5% to 56.5% by weight (see Figure 2.7). The transformation temperatures are 
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Strain 

Figure 2.5 Stress-strain curves for a shape memory 
alloy at various temperatures, after Funakubo [33]. 
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Figure 2.6 Stress-strain curves for a shape memory alloy at 
various temperatures, after Wayman and Duerig in [32]. 
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Figure 2.7 Variation in austenite start temperature with 
weight % Ni in Nitinol, after Buehler and Wang [1]. 


also strongly influenced by stress. It is known from early work, such as Cross et al. [28] for 
example, that the best shape recovery performance for fully annealed Nitinol is achieved with an 
anneal temperature of approximately 930 °F. Goldstein et al. [38] found that partial annealing 
of the residual stresses due to cold work can be used to change the transformation temperatures 
and induce a rhombohedral R phase. 

All of the discussions to this point have assumed that the alloy is in a fully annealed state. 
Some annealing is necessary to remove some of the cold work dislocations which impede the 
twin boundary movement, but it is usually desirable to only partially anneal the alloy and leave 
some of the dislocations for strength and durability. Production of the R-phase during cooling 
from austenite is a direct result of partial annealing, as observed by Saburi in [34]. The R-phase 
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is itself a product of a thermoelastic martensitic transformation and may be of interest for some 
applications because of its lesser hysteresis and reversible strain. The R-phase is simply a by- 
product of circumstances for the focus of the work presented here and does not significantly 
affect the performance of alloy. The effects of the R-phase are measurably manifested in all 
aspects of the alloy characteristics, however. For example, stress-strain curves may show an 
additional knee due to R-phase rearrangement under certain conditions, and DSC plots exhibit 
an additional endothermic peak due to the formation the R-phase. 

The Nitinol actuator material used in this study is partially annealed in the as-received 
condition, so the material will exhibit the presence of R-phase. The actuators will be used in a 
one-way mode of constrained recovery operation, but characterization of the dependency of the 
material on thermal cycle will be quantified in lieu of one-way training. The Nitinol actuators 
will be embedded within a composite matrix, along the direction of the matrix fibers, and bonded 
within the matrix during cure. Constrained recovery operation will be insured by mechanical 
constraint on the resulting SMAHC specimens. 

2.2 Constitutive Modeling 

A constitutive model of the SMAHC material system is needed in order to predict the 
thermomechanical behavior of structures manufactured from the SMAHC material. Constitutive 
theories for SMAHC materials could be developed from micro-structural shape memory models 
(e.g., [24, 25, 26, 27]), as was done by Boyd and Lagoudas [36] and Birman [37]. However, 
these models are difficult to use in practice because they are qualitative in nature and depend 
upon determination of impractical parameters/material properties. An alternative approach is to 
employ a constitutive model which makes use of experimental measurement of basic engineering 
properties. Two such models were developed as part of the present study and will be discussed in 
this section. The models are developed by a mechanics-of-materials approach (see Jones [39]). 

A representative volume element of a SMAHC lamina is shown in Figure 2.8. The element 
is taken to be in the plane of the plate. Note that the matrix material is a composite material 
itself, for example graphite-epoxy, with principal material directions 1 and 2, and the SMA 
actuator is embedded in the 1 -direction. 

Consider the uniaxial state of stress in a SMA actuator along the actuator axis. One form for 
the thermoelastic constitutive relation of the SMA actuator in the 1-direction was proposed by Jia 
and Rogers [5] and was extended to include thermal expansion by Turner, Zhong, and Mei [16]: 

v\ u = E „( i +<r r T > .'I, 

cr la = E a {ei - a la AT) T < A„ 

where E„ is the Young’s modulus of the SMA, 1 1 is the strain in the 1 -direction, a T is the recovery 
stress of the activated SMA, and ai„ is the thermal expansion coefficient for the inactivated SMA. 
Note that the Young’s modulus and recovery stress of the SMA are functions of temperature. 
The recovery stress is also a function of the initial strain. The thermal expansion coefficient can 
be treated as essentially constant for temperatures below the phase transformation temperature. 
A detailed derivation of the constitutive relations for a SMAHC lamina using this approach, 
referred to as the recovery stress model (RSM), is given in Appendix A. A complete finite element 
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Figure 2.8 Representative volume element for a SMAHC lamina. 


formulation was developed by Turner, Zhong, and Mei [16] for the thermomechanical response 
of SMA hybrid composite panel-type structures based upon the RSM constitutive relations. The 
RSM constitutive relations and the resulting finite element governing equations are in a form that 
is physically intuitive, but not consistent with conventional thermoelastic analysis and, therefore, 
not easily adapted to commercial finite element codes. It will also be shown that the RSM 
misses some of the mechanics of the material system. Thus, the constitutive relations will be 
re-cast in a more familiar and robust form. 

The uniaxial thermoelastic constitutive relation for the SMA can be written in terms of an 
effective coefficient of thermal expansion (CTE): 


Via 


= E a 


T 



Ola(r )(It 


( 2 . 2 ) 


The way in which this relation is used depends upon the application. A measure of the nonlinear 
CTE o 1( ,(r) over the temperature range of concern would be appropriate for free or restrained 
recovery applications. For constrained recovery applications, however, one must resort to 
measurement of the recovery stress. 

Comparison of equation (2.2) with equations (2.1) shows the functionality of the effective 
CTE for constrained recovery applications; it has the usual effect for temperatures below the 
austenite start temperature 
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and is related to the recovery stress at temperatures above /l. s by the equation 

T T 

a r = -E„ J (\u,(t)(It or j (\ h ,{t)(It = -y- (2.4) 

To To 

Note that in this case, the nonlinear thermoelastic nature of the SMA is still captured, albeit in 
a different way, because measurements of recovery stress and modulus versus temperature are 
inherently cumulative (integrated). Also, note that an experimental approach could be devised 
to measure appropriate recovery stress <r r quantities for restrained recovery applications as an 
alternative to direct measurement of O| 0 (T). 

Thus, the general form for the thermoelastic constitutive relation in the 1 -direction for a 
prestrained SMA actuator can be written as 

/ 

<7|« = E„ (i - j o u,(t)(It and sgn(oi„) = {^j ^ < J' (2.5) 

where o | (l is due to thermal expansion only for T < .4, and is determined by direct measurement 
or from equation (2.4) for T > 4,. 

In order to derive the thermoelastic constitutive relation for the 1-direction of a SMA hybrid 
composite, it is assumed that a stress a \ acts alone on the element shown in Figure 2.8 {a> = 0) 
and that the SMA and composite matrix are strained by the same amount, tj (i.e., plane sections 
remain plane). The 1 -direction thermoelastic constitutive relation of the SMA actuator is given 
by equation (2.5). Similarly, the one-dimensional, stress-strain relation in the 1 -direction for the 
composite matrix can be expressed as 

T 

<7\m = Elm tl - j ay m {T)(lT (2.6) 

to 

Temperature dependence of the matrix thermal expansion coefficient has been included for 
generality. 

The resultant force in the 1 -direction (ao = 0) is distributed over the SMA and composite 
matrix and can be written as 

<71,4] - rr Ul A a + rr | ,4 (2.7) 

where ( <r , , .4 1 ) , (a\ u . ,4„), and ) are the (stress, cross-sectional area) of the entire 

element, SMA actuator, and composite matrix, respectively. Thus, the average stress cr\ is 

(7\ — (?]„ v u -f- o'\ m v rri (2.8) 

where r u = .4„/.4] and v,„ = A m /Ai are the volume fractions of SMA and composite matrix, 
respectively. Substituting for o\, <ti„, and < t\„ , results in 
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So, the Young’s modulus and thermal expansion coefficient for the hybrid composite in the 
1 -direction have the form 
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A similar constitutive relation may be derived for the 2— direction by assuming that the 
applied stress a> acts upon both the SMA and the matrix (<T| = 0). The uniaxial thermoelastic 
constitutive relations in the 2— direction for the SMA actuator and the composite matrix are 
analogous to those in equations (2.2) and (2.6), i.e., 
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respectively. The CTE in equation (2.12) for the SMA does not depend on the recovery stress, 
but is still nonlinear due to the differing martensitic and austenitic properties. 

The total elongation is due to strain in the composite matrix and the SMA actuator and may 
be written in the form 
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Thus, the total strain becomes 
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Substituting for o, f2 0 , an d (2m results in 
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Therefore, the Young’s modulus and thermal expansion coefficient for the hybrid composite in 
the 2-direction have the form 
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and 
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Expressions for the hybrid composite Poisson’s ratios and shear moduli follow from similar 
derivations: 
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The thermoelastic constitutive relations in principal material coordinates for an orthotropic 
material have the following form 
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where { t } is the total strain, {t} w is the mechanical strain, {<} 7 is the thermal strain, and [77] 
is the compliance matrix. These relations can be inverted to give 
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where [f] is the stiffness matrix. The constitutive relations, for plane stress, reduce to 

/{- 


( 2 . 21 ) 


1 


'Qu 

Qu 

0 

fM 1 


Qxx 

Qu o 

,2 

> = 

Qu 

Q 12 

0 

** 

> — 

Qu 

Qn 

0 

Ju j 


0 

0 

1 

Cy 

111-’ j 


0 

0 

Q 66 


>>(r) )dr 
0 


( 2 . 22 ) 


where [(,)] is the reduced stiffness matrix of the composite lamina. The reduced stiffnesses are 
related to the engineering properties through the following equations 
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The expressions for the hybrid lamina engineering properties in terms of the properties of the 
constituent materials are reiterated here for convenience. 
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The thermoelastic constitutive relations for a general k-th orthotropic layer of a laminate 
with an arbitrary fiber/SMA actuator orientation angle #, relative to the 1 -direction in Figure 
2.8, follow from transformation of the relations in principal material coordinates, equation (2.22). 



(2.25) 


or 


{<*}* - “ {QtQ (2.26) 

The transformed reduced stiffness matrix [Q] k is related to the reduced stiffness matrix in 
principal material coordinates by 

[Q], = rn _ 1 [Q],rn - 7 < 2 - 27) 

and [7 ] is the transformation matrix given by 

cos" # 

[T] = sin 2 # 

—2 sin 0 cos # 

The constitutive relations given by equations (2.22) or (2.26) will be referred to as the effective 
coefficient of thermal expansion model (ECTEM). 

The constitutive equations in principal material coordinates, equations (2.22), will now be 
considered further to clarify some advantages of the ECTEM approach. The expressions for the 
thermal strains in equations (2.24) allow the constitutive relations to be rewritten as 
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where the thermal effects of the SMA and matrix material have been separated and 
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are constants. These constitutive relations have a form similar to those in the RSM [16], which 
are rewritten here for convenience of comparison: 
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While the RSM only allows the recovery stress to affect the stress in the 1 -direction of a lamina 
(term 2 of equation (2.31 )), the analogous term in the ECTEM (term 2 of equation (2.29)) shows 
an effect of the 1 -direction effective CTE on the stress in the 2-direction of the lamina, through 
the reduced stiffness Q i_>. The ECTEM also allows for the fully nonlinear functionality of the 
2-direction CTE for the SMA while the RSM neglects it altogether for 7 > A s . As a result 
of these differences, the thermal forces calculated by the two formulations are different and 
the ECTEM more accurately captures the behavior of the material system. The form for the 
ECTEM constitutive relations given in equation (2.29) will be useful for reference in the thermal 
buckling discussion of Chapter 5. 

In summary, the ECTEM essentially reduces to the RSM for constrained recovery applica- 
tions, but captures the two-dimensional behavior of SMA actuators and their interaction with 
a composite matrix more accurately than the RSM. The ECTEM is also capable of modeling 
restrained or free recovery behavior in a similar manner with experimental measurements of 
SMA recovery stress, modulus, and/or nonlinear CTE. 


2.3 Deformation Behavior and Stress Resultants 

A model for the deformation behavior of a SMAHC laminate will be necessary, along with 
the constitutive relations, to develop the equations governing the static and dynamic response 
of such structures. Recall that one demand on the formulation is to model the thermoelastic 
response in the presense of immovable in-plane boundaries, i.e., thermal buckling and post- 
buckling analysis. Thus, the midplane strain and curvature vectors are defined from the von 
Karman strain-displacement relations [40]: 

j = {<"}+-» 

= {C) + M) + --M <2-32) 



to allow for potential membrane stretching effects due to geometrically nonlinear deflections 
associated with thermal post-buckling. Note that the midplane strain {t 0 } consists of two 
components; the midplane strain due to membrane displacements { < ° m } and the nonlinear 
midplane strain due to bending displacements 

Substitution of the strain-displacement relations in equation (2.26) results in 

W* = + ~{ / ' }) - [Q] * } Tk (2.33) 

The resultant force and moment vectors of the SMA hybrid composite laminate are defined as 
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where the laminate stiffness [ 4], [B] and [/)] and the thermal resultant force and moment vectors 
{ \ A7 } and {Mat} are all temperature dependent. The laminate stiffness matrices have the 

usual form 
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and the expression for the thermal resultants is as follows. 
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Expanding the resultant force and moment relations, equation (2.35), in terms of the midplane 
strain and curvature components, equation (2.32), results in 

{A-} = (/ilK„} + |.'t]K} + [«])») -(Vii'l 

= (/V,„) + {A«} + {Ai.)-!Aa7') 


(2.37) 
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a {M} = [#]{<;;,} + m<i) + [«H«} - W*r} (7 39) 

= {M,„} + {Me} + {M h } - {Mat} 

where the subscripts m, 0, b, and A T indicate components induced by membrane strain, nonlinear 
membrane strain, bending strain, and thermal expansion, respectively. The resultant force and 
moment relations and the strain-displacement relations presented in this section will be useful 
for reference in formulating the equations governing the thermomechanical response of SMAHC 
panel-type structures in Chapter 4. 

A review of the behavior of SMAs was given in this chapter followed by development 
of a constitutive model for SMA actuators embedded in a composite matrix and used in a 
constrained-recovery model of the SME for adaptive stiffening. In order to use this model to 
predict the behavior of the SMAHC material system used in this study, the thermomechanical 
properties of the constituents must be determined. Results from such characterization testing will 
be presented in the next chapter, following motivation of the selected materials, description of the 
SMAHC material system configuration, and discussion of the processes involved in fabricating 
the SMAHC specimens. 
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Chapter 3 

Specimen Fabrication and Material System Characterization 


A constitutive model was developed in the previous chapter to describe the mechanics of a 
SMAHC material system. That model requires measurement of a number of material properties 
from the constituent materials in order to quantify the behavior of a structure fabricated from such 
a material. In this chapter, a brief motivation is given for the selected material system, followed 
by details of the SMAHC specimen fabrication processes and hardware. Processes involved in 
and results from thermomechanical tests on the constituent materials will be discussed, and the 
chapter will conclude with disclosure of the final material property data set used to generate the 
numerical results for SMAHC specimens in this study. 

3.1 Material System Selection 

In the typical application being considered, such as adaptive stiffening of aerospace struc- 
tures, the SMA elements would be activated through an inherently elevated temperature of the 
service environment. Thus, activation of the SMA by surface heat in the planned experiments 
would be the most accurate simulation. Surface heating could be accomplished through the use 
of quartz radiant heaters or a number of other sources. These heating methods involve practical 
problems, however, such as the need for thermal insulation of the test article support structure and 
other equipment, interference with structural/acoustic/thermal measurements, thickness-direction 
asymmetry of the thermal load, and a complicated thermal control system. Alternatively, the 
SMA can be activated electrically through resistive heating. This method has the advantages of 
simple operation, relatively uniform heating, and excellent controllability. Electrical isolation of 
the SMA elements may be required for electrical activation of laminates involving carbon-epoxy 
matrices. Electrical isolation can be achieved through coating techniques, but this could present 
significant complications. It was decided that a glass-epoxy matrix would be used to circumvent 
these problems as experiments performed in this study are intended for demonstration and theo- 
retical validation only. This material system also affords the advantage of visual flaw detection 
due to its translucency. 

Most recent SMA hybrid composite structural designs have focused on SMA wire-type 
actuators distributed uniformly throughout the laminate or individual laminae. There are several 
disadvantages associated with the use of SMA wire actuators. Integration of wire-type actuators 
into a composite leads to complicated fabrication procedures because of the relatively large 
number of actuators usually needed. The small surface area of the actuators makes them very 
sensitive to actuator/matrix interface flaws because voids can be of significant relative size. Small 
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size has also been blamed for a relatively high rate of actuator breakage during cure because of 
sensitivity to stress concentrations at the mechanical restraints. Finally, it is difficult to achieve a 
desirable overall volume fraction of SMA when trying to optimize the integration of the actuators 
by placing them in only selected layers. Therefore, the present work will focus on SMA hybrid 
composite laminates involving ribbon-type SMA elements to simplify the fabrication procedure, 
allow for more flexibility in fabrication, and hopefully desensitize the actuators to interface voids 
and stress concentrations. For example, the ribbon-type SMA elements can be placed in strips 
cut out of particular matrix laminae to result in a rather simple lay-up operation, while keeping 
the volume fraction high in the desired locations. Another potential benefit of this type of lay-up 
procedure is that it could potentially be extended directly to e-beam consolidation procedures 
where spools of SMA ribbon would strategically replace spools of matrix material ribbon. 

In summary, the material system used in this study consists of a glass-epoxy matrix material 
with embedded Nitinol ribbon. Details of the SMAHC specimen fabrication procedures/hardware 
are given in the following section, followed by a discussion of the thermomechanical testing of 
the constituent materials. 


3.2 Specimen Fabrication 

Although a panel-type specimen is closer to the configuration of intended applications, a 
more simple test specimen will suffice for the objectives of this study; demonstration of the 
thermomechanical response of a SMAHC laminate and experimental validation of the theoretical 
model developed in this work. A schematic of the selected beam-type specimen is shown in 
Figure 3.1. Various constraints, other than the use of ribbon-type actuators, influenced the design 
of this laminate. These constraints include fabrication tooling geometry and dynamic range of 
anticipated experimental configurations. The specimen design is by no means optimized, but the 
predicted performance is in a range that is suitable for the expected loading conditions and to 
demonstrate significant improvement from the SMA reinforcement. 

Note that this specimen design calls for a SMA actuator width of 0.0127 m (0.5 inch) to be 
embedded within each 0° glass-epoxy layer, i.e., replace a portion of each 0° layer. Material 
availability limited the procured SMA ribbon cross section to 0.0023x1. 5e-4 m (0.09x0.006 
inches). This was considered acceptably close to the glass-epoxy unidirectional tape thickness 
(~1.27e-4 m, ~0.005 inches) and it was planned that five widths of the ribbon would be placed 
side-by-side to develop a width near the half-inch in the specimen design. Alternatively, in the 
case that the adjacent SMA actuators prove detrimental because of a weak interface between 
actuators, the ribbons could be separated within laminae or between laminae. 

The SMAHC specimen fabrication procedure used in this study consists of the following 
steps; prestrain the SMA actuators, lay-up the composite laminate, cure the laminate while 
restraining the actuators, and machine beam specimens from the consolidated SMAHC laminate. 
The Nitinol material was packaged on a spool, which subjected the ribbon to significant out- 
of-plane bending strain. Some in-plane bending, axial, and torsional strain was also imparted 
to the ribbon during the spooling operation. This packaging strain must be removed prior to 
imparting the actuator prestrain. Packaging strain removal can be conveniently done by inducing 
the reverse (austenitic) transformation. An apparatus was developed to perform this operation, 
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22"x1 "x0.080" G lass- E poxy/S M A Hybrid Beam (6 Count) 
(45/0 SMA /-45/90) 2sym Lamination 
SMA in 0° Laminae Only 



Cut Beams from Larger Laminate 



SMA Ribbon: 

0.5'‘x0.005" Cross Section 
6" Overhang at Beam Ends 


Figure 3.1 Schematic of beam specimens cut from a SMA hybrid composite laminate. 

shown in Figure 3.2. It is desirable to induce the martensitic transformation just prior to 
prestraining in order to avoid aging effects, such as the rubber-like form of pseudoelasticity. 
Inducing the transformation is also accomplished during the packaging strain recovery process 
as the ribbon cools to ambient temperature. It will be shown in the next section that, for this 
particular Nitinol composition in an unrestrained condition, the martensitic transformation is not 
complete at ambient temperature. However, the martensite will be stress induced during the 
prestraining operation. Stress— inducing the martensite results in essentially the same state as 
would be achieved by cooling the ribbon to sub-ambient temperatures prior to prestraining. The 
main difference is that there may be some austenite remaining in the ribbon after prestraining 
by the former method, but should be of no consequence because any subsequent thermal cycling 
(e.g., composite cure) will render the state from the two methods equivalent. The former method, 
without cooling to sub-ambient temperatures, was selected for simplicity. 

The prestraining operation needs to be very precise and repeatable. It was decided that 
the actuators would be elongated on a mechanical loading machine operated in stroke-control 
mode. Alignment of the ribbon during elongation is critical to avoid bending strain and to 
obtain repeatable results. Handling of the ribbon is cumbersome during installation in the test 
machine, and it is desirable to prestrain multiple strands at once to expedite the process. In 
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Figure 3.2 Packaging strain recovery rig. 


order to rectify these issues, an apparatus was developed to align and maintain proper geometry 
on three ribbon lengths at once during installation in the test machine. A close-up view of the 
alignment/gripping apparatus is shown in Figure 3.3. The entire assembly (alignment/gripping 
apparatus, spacer, and three ribbons) is shown installed in the test machine in the left half of 
Figure 3.4. The Nitinol ribbons are shown ready for elongation, subsequent to spacer removal, 
in the right half of Figure 3.4. 

The particular spacer length shown was used to maintain a SMA ribbon gage length of 
0.635 m (25 in) for manufacture of the beam specimens shown in Figure 3.1. The gage length 
of the ribbons was elongated to a length of 0.6604 m (26 in), which corresponds to a prestrain 
of 4%. This gage length was used in order to ensure uniform prestrain over the 0.6604 m (26 
in) span between mechanical restraints on the composite lay-up tooling, which will be discussed 
subsequently. A total of fifty-four ribbon lengths were prestrained 4% using this procedure to 
prepare for SMAHC laminate fabrication. 

A mechanical drawing of the lay-up and cure tooling plate is shown in Figure 3.5. This 
apparatus was designed to provide a flat surface (± 0.005 inches) and mechanical grips to restrain 
the SMA during an autoclave cure procedure. Recall that the beam specimen design calls for a 
lamination of (45/0/-45/90h s with the Nitinol in the 0° layers. This lamination was accomplished 
by cutting widths of the 0° layers to size during assembly and filling the void with five ribbon 
widths laid side-by-side. A picture of the lay-up procedure at the stage of finishing the first 
0° layer is shown in Figure 3.6. Note that Kapton release film, fiberglass breather cloth, 
and teflon-coated fiberglass bleeder cloth are layered between the tooling plate surface and the 
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Figure 3.3 Nitinol ribbon gripping apparatus for prestraining and mechanical testing. 
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SMA Hybrid Composite Tooling Plate 


26.000- 


5.500- 



■ Knurled surface between edge of plate and dowels 
O 0.375 -diameter X 0.750-length dowels 
o Tapped holes for 1/4-28 bolts 


Material 

Plate, 17-4 PH stainless 
DoWels, 300 series stainless 
Instructions; 

Top surface of plate should be ground 


Figure 3.5 Mechanical drawing of the laminate lay-up 
and cure tooling apparatus (all dimensions in inches). 
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Figure 3.7 Completed laminate with cure control thermocouples and NiTi ribbon grips torqued. 

first glass-epoxy layer, in that order. The laminate after complete assembly is shown in Figure 
3.7, where the cure cycle control thermocouples have been installed near the laminate and the 
Nitinol ribbon ends have been secured within the knurled mechanical restraints. The restraints 
consisted of three parts (tooling plate surface, a thin spacer with knurling on both sides, and 
the top clamping bar) to minimize bending of the actuators between the laminate boundary and 
the restraint and to maximize grip for each layer of SMA by minimizing the number of layers 
gripped together. The mechanical restraint screws were torqued to 11.30 N«m (100 in*lbf). The 
release paper was subsequently removed from the top layer of glass-epoxy and the laminate was 
topped with bleeder cloth, Kapton film, a caul plate, bleeder cloth, and breather cloth, in that 
order. The entire assembly was vacuum bagged as shown in Figure 3.8 and drawn to an internal 
pressure of 0.0475 atm. The assembly was subjected to the autoclave cure cycle recommended 
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Vacuum Port 


Perimeter Seal 


Figure 3.8 Laminate after vacuum-bagging and applying vacuum. 



Figure 3.9 Consolidated S MAFIC laminated panel. 



Figure 3.10 SMAHC beam specimens machined from panel. 


by the vendor (Fiberite, Inc.) and the resulting consolidated part is shown in Figure 3.9. Details 
of the cure cycle will be discussed in the following section. Two SMAHC beam specimens were 
machined from this laminate using the respective centerlines of the embedded Nitinol ribbon. 
The beam dimensions are shown in Figure 3.10, where it can be seen that the beam marked 
number 3 has been readied for installation in test fixtures. The beams have dimensions of 
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Table 3.1 Nominal properties of a Fiberite E-glass/934 epoxy lamina (50-55% fiber volume). 


Property 

Room Temperature 

Young’s Modulus, E| 

41.37-55.16 GPa (6-8 Msi) 

— 

Poisson’s Ratio 

0.28-0.32 

Density 

1.919-l.980e+3 kg/m 3 (1.795-1.852e-4 lbf s 2 /in 4 ) 


Cure Temperature 

176.7°C (350°F) 

Glass Transition 

193.3°C (380°F) 


0.5588x0.0254x0.0019 m (22x1x0.078 in). The SMA leads of the test-prepared beam were 
trimmed to an overall length of 0.6604 m (26 in). The overall volume fraction of the SMA 
within the glass-epoxy dimensions is approximately 13.8%. 


3.3 Material System Characterization 

The material system and processes used to fabricate SMAHC beam specimens was described 
in the previous sections. In order to predict the thermomechanical behavior of such a structure 
using the constitutive model developed in the previous chapter, a number of material properties 
from the constituents must be determined. Characterization testing that was performed on the 
glass-epoxy matrix material and Nitinol alloy are described in the following sections. 

3.3.1 Glass-Epoxy Matrix 

Glass-epoxy unidirectional prepreg tape was obtained from Fiberite, Inc. The material 
is Hye 8 End/934, an E-glass/934 epoxy resin material system. The resin material is for- 
mulated for autoclave processing and the following cure procedure is recommended; heat at 
1. 1-2.8° C/minute (2-5°F/minute) to 121. 1°C (250°F) and hold for 15 minutes, apply 689.5 kPa 
(100 psig) pressure and continue to hold at 121.1°C (250°F) for an additional 45 minutes, heat 
at 1.1 -2 .8° C/minute (2-5°F/minute) to 176.7°C (350°F) and hold for 120 minutes, and cool 
at 2.8°C/minute (5°F/minute) (maximum) to ambient. The material and processing properties 
shown in Table 3.1 were available from the Fiberite Materials Handbook [41]. Although these 
properties are useful as estimates, thermal expansion and more precise elastic property char- 
acterization are needed for accurate prediction of thermomechanical performance. Therefore, 
significant thermomechanical testing of the glass-epoxy material system was performed. 

Two laminates were fabricated from the glass-epoxy material system to provide specimens 
for thermomechanical testing; a unidirectional lay-up (Oho and an angle-ply lay-up (±45)4. Tab 
material was post-cured to the laminates and the resulting structures were machined into uniaxial 
tensile specimens with dimensions of 0.2286x0.0254 m (9x1 inches) (gage length of 0.1524 m, 
6 inches). The tensile specimens were instrumented with a “T-strain gage" at the mid-span 
and tested on a MTS machine. The unidirectional specimens were tested according to ASTM 
standard D 3039/D 3039M for estimates of the Young’s moduli and Poisson’s ratios; E \ , E>, i/i >. 
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and !'■> i . The angle-ply specimens were tested according to ASTM standard D 35 1 8/D 35 1 8M for 
estimates of the shear modulus (> \ > . Although the ASTM standards call tor modulus estimates 
to the determined over strain ranges of ~ 1000-3000 fit and ~ 2000-6000 //(, respectively, the 
transverse (90°) uniaxial and ±45 angle-ply stress-strain relations showed discernible nonlinearity 
over those ranges. The maximum normal strain and shear strain for the test configurations 
considered in this study are not expected to exceed 500 and 250 //t, respectively. Therefore, E> 
and (i i > estimates were determined from ~50— 1000 //( and ~50— 500 //<-, respectively, to achieve 
more representative estimates with ample headroom. Results from these tests/calculations are 
given in Table 3.2. 


Table 3.2 Mechanical properties of Fiberite E-glass/934 
epoxy material system in principal material coordinates. 


Temperature 

E, 

e 2 

"12 

G 2 i 

°C 

GPa 

GPa 

— 

GPa 

23.9 

49.30 

20.00 

0.29 

9.65 

37.8 

49.16 

19.44 

0.29 

9.24 

65.6 

48.75 

18.20 

0.29 

8.41 1 

93.3 

48.61 

16.20 

0.29 

7.58 

121.1 

48.54 

14.00 

0.29 

5.58 

148.9 

48.82 

11.38 

0.29 

3.45 


Specimens for measurement of the coefficient of thermal expansion (CTE) were also ma 
chined from the unidirectional laminate (0)2o described above. Two CTE specimen geometries 
were machined for testing by two different methods. Rectangular 0.1524x0.006 m (6x0.25 inch) 
specimens were tested according to ASTM standard E 228 at Materials Research and Engineer- 
ing, Inc. A summary of those test results is given in Arvidson et al. [42]. Specimens with the 
geometry shown in Figure 3.1 1 were tested by an interferometric dilatometer method at NASA 
LaRC. Both quadratic and cubic interpolation of the raw thermal expansion data are used in 
practice to generate linear and quadratic estimates of the CTE with temperature, respectively. 
Cubic interpolation was employed with the ASTM E 228 test method while quadratic interpola- 
tion was used with the interferometer method. The consequence of this difference can be seen 
in Figures 3.12 and 3.13, where the longitudinal (fiber direction) and transverse CTE estimates 
versus temperature are shown from the two methods. The two methods are in agreement 
for the overall trend in the CTE data, but the cubic interpolation of the data collected with the 
ASTM standard captures higher order effects. Data resulting from both methods are tabulated in 
Table 3.3, but the quadratically varying CTE data from the ASTM E 228 method will be used 
in the predictions reported in this study. 
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Inter f eroneter Dilatometer Specimen 
LaRC CTE Specimen 



Notes: 

1) Specimen length is critical 

2) All dimensions in inches 

Figure 3.11 Schematic of the CTE specimen geometry for 
measurement by the interferometric dilatometer method. 
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Figure 3.12 Fiberite E-glass/934 epoxy longitudinal 
(fiber-direction) CTE vs. temperature by two methods. 


300 
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Figure 3.13 Fiberite E-glass/934 epoxy transverse CTE vs. temperature by two methods. 
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Table 3.3 Fiberite E-glass/934 epoxy CTE in principal 
material directions as a function of temperature. 



Ol 

0 2 

Temperature 

ASTM E228 

Interferometer 

ASTM E228 

Interferometer 

°C 

U 

•<U 

// e/°C 

/ie/°C 

/if/ 0 C 

23.9 

— 

6.17 

— 

21.04 

26.7 

5.83 

6.21 

13.81 

21.13 

37.8 

6.95 

6.34 

21.02 

21.49 

48.9 

7.22 

6.44 

24.14 

21.85 

60.0 

7.00 

6.55 

24.66 

22.23 

71.1 

6.62 

6.68 

23.90 

22.59 

82.2 

6.34 

6.79 

22.79 

22.95 

93.3 

6.25 

6.91 

22.00 

23.31 

104.4 

6.39 

7.02 

21.92 

23.67 

115.6 

6.77 

7.13 

22.66 

24.05 

126.7 

7.18 

— 

24.01 

— 

137.8 

7.42 

— 

25.47 

— 


Tensile specimens with a geometry of 0.1524x0.0127 m (6x0.5 inches) were machined from 
scrap (glass-epoxy-only) material remaining subsequent to machining SMAHC specimens from 
the larger quasi-isotropic laminate shown in Figure 3.9. A schematic of the specimen machining 
organization from the SMAHC panel is shown in Figure 3.14. The tensile specimens were 
tabbed and instrumented with a “T-strain gage” at the mid-span prior to testing. These quasi- 
isotropic specimens were tested according to ASTM standard D 3039/D 3039M for estimates 
of the extensional stiffnesses; An and A 22 • A brief discussion of these measurements is given 
in the following paragraph. 

Recall the stress resultant relations given in equation (2.35). In performing the uniaxial tensile 
tests on the quasi-isotropic specimens, the applied force N r , axial strain t T and transverse strain 
t were measured. The transverse and shear force components are zero, A y = A J(/ = 0. 
For a symmetric quasi-isotropic laminate, the shear-extension coupling stiffnesses are zero, 
/l 1(i :: A 26 = 0 and the extensional stiffnesses are equal A\\ = A_>j. With these simplifications, 
the stress resultant equations reduce to 

A. i = A \\(° r + A i 

0 = Ai2t” + A\\t°y (3.1) 

0 = Am-,% 
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Clean Edges for overall dimensions of 22±0.005'xl2±0005' 
22±0.005'xl ±0.005' beans Wi th 0.5' ribbon Width centered on CL 
22±0.005'xl±0.005' bean specimens without ribbon 
6±0.005'x0.5±0.005' longitudinal mechanical test specimens 
6>±O.OO5'xO.5±O.0O5' transverse mechanical test specimens 
3.017+0. 005'-0.0000'x0.75±0.005' longitudinal CTE test specimens 
3.017+0. 005'-0.0000'x0.75±0.005' transverse CTE test specimens 

1) Machine two 22x1 inch beams from panel using centerline 
of metallic ribbon 0.5'-width as centerline of each beam. 

DD NDT CUT OFF RIBBGN ENDS PROTRUDING FROM PANEL EDGES. 

2) Machine two 22x1 inch beams w/o ribbon from two outer 
str ip scraps just removed from panel. 

3) Machine tensile specimens and thermal expansion specimens 
from interior scrap remaining from i>. See tensile and CTE 
specimen detail drawings. 


Figure 3.14 Schematic of the machining process for removing tensile, 

CTE, and SMAHC beam specimens from the SMAHC laminate. 

The relationship in the second equation .4 12 = - .4 n <'!, / <0 , allows the first equation to be rewritten 



Thus, the measured quantities are sufficient to determine the extensional stiffnesses of the quasi- 
isotropic laminate. 

The uniaxial and angle-ply properties discussed previously can be used with classical 
lamination theory (CLT) to calculate extensional stiffnesses for the quasi-isotropic laminate. The 
extensional stiffnesses measured from the quasi-isotropic specimens are shown in comparison 
to those predicted with CLT from the uniaxial and angle-ply properties in Table 3.4. It can 
be seen that good agreement is achieved. Thus, usage of the properties in Table 3.2 is deemed 
appropriate for predictions of the glass-epoxy material system behavior. 
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Table 3.4 Comparison of measured and predicted extensional stiffnesses 
for a Fiberite E-glass/934 epoxy quasi-isotropic (45/0/-45/90) 2s laminate. 



Extensional Stiffness, An 

Temperature 

Measured 

Predicted 

°C 

MN/m 

MN/m 

23.9 

69.2 

63.5 

37.8 

63.3 

62.4 

65.6 

60.1 

60.1 

93.3 

57.5 

57.3 

121.1 

53.8 

53.2 

148.9 

50.7 

48.9 


A variety of test procedures and results have been presented in this section. A summary of 
the properties to be used in the numerical predictions later in this study are shown in Table 3.5. 
Note that slight extrapolation of the CTE data and linear interpolation of the CTE and elastic 
property data was performed to result in a common temperature resolution. 

3.3.2 Nitinol Alloy 

Nitinol ribbon material was provided by Shape Memory Applications, Inc. A Nitinol alloy 
was specified with an appropriate composition to give a martensite finish Mj temperature above 
the expected ambient temperature while maintaining an austenite finish temperature well below 
the glass-transition temperature of the Fiberite 934 epoxy matrix. The achievable ribbon cross 
section was 0.0023x0.00015 m (9x0.006 inches). The nominal transformation temperatures of 
the provided alloy, in a fully annealed state, were reported by the vendor as shown in Table 3.6. 
These temperatures were reportedly determined by differential scanning calorimetry (DSC). As 
the provided material is not in a fully annealed state, it was expected that all of the transformation 
temperatures would be somewhat lower than those in the table due to dislocations remaining 
from the cold work performed to generate the final geometry. Conversely, mechanical restraint 
in the planned test configurations will increase phase transformation temperatures. In any case, it 
was decided that additional DSC measurements should be performed on samples of the provided 
material to determine the baseline phase transformation temperatures. 

Small samples (clOmg) were cut from the ribbon material in the as-received condition for 
the DSC testing. Initial scans showed irregularities in the endothermic peak during heating 
and a disproportionately small-sized exothermic peak during the cooling, which indicated a 
potential for unexpected additional activity below ambient temperature. As a result of these 
initial findings, two more-rigorous scanning procedures were defined in an effort to differentiate 
between several effects seen in the initial scans. The first scan procedure, termed the heat-first 
procedure, entailed heating from ambient temperature to 150 C at 5 C/minute, cooling from 
150°C to -20°C at 5°C/minute, and repeating this entire profile once the sample was returned to 
ambient from -20°C (also at 5°C/minute). For future reference, the first cycle will be called ,/zrsr 
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Table 3.5 Summary of the thermomechanical properties of Fiberite 
E-glass/934 epoxy material system in principal material coordinates. 


Temperature 

Ei 

e 2 

A' 12 

G 21 


°2 

°C 

GPa 

GPa 

— 

GPa 

/u/°C 

/if/ 0 c 

21.1 

49.30 


0.29 

9.65 

5.22 


26.7 

49.30 

19.86 

0.29 

9.58 

5.83 

13.81 

37.8 

49.16 

19.44 

0.29 

9.24 

6.95 

21.02 

48.9 

49.02 

18.96 

0.29 

8.89 

7.22 

24.14 

60.0 

48.82 

18.48 

0.29 

8.55 


24.66 

65.6 

48.75 

18.20 

0.29 

8.41 

6.80 

24.28 

71.1 

48.75 

17.79 

0.29 

8.27 



82.2 

48.68 

17.03 


7.93 


22.79 

93.3 

48.61 

16.20 

0.29 

7.58 

6.25 

22.00 

104.4 

48.61 

15.31 

0.29 

6.76 

6.39 

21.92 

115.6 

48.54 

14.41 

0.29 

6.00 

6.77 


121.1 

48.54 


0.29 

5.58 

6.97 

23.33 

126.7 

48.61 



5.17 

7.18 

24.01 

137.8 

48.68 

12.41 

0.29 

4.27 

7.42 

25.47 

148.9 

48.82 

11.38 

0.29 

3.45 

7.63 

26.28 


Table 3.6 Transformation temperatures for the provided Nitinol alloy in a fully-annealed state. 



Temperature 

Transformation Marker 

°C 

°F 

Austenite Start 

67.2 

153 

Austenite Peak 

80.0 

176 

Austenite Finish 

85.0 

185 

Martensite Start 

46.1 

115 

Martensite Peak 

40.0 

104 

Martensite Finish 

32.2 

90 


heat and the repeat cycle will be called second heat. The second procedure (cool-first) consisted 
of cooling from ambient temperature to -20° C at 5° C/minute, heating from -20° C to 150°C at 
5°C/minute, cooling from 150°C to ambient temperature at 5°C/minute, and repeating this entire 
profile once the sample was returned to ambient. For this case, the first cycle will be called first 
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Figure 3.15 Heat-first DSC signatures for first- and second-heat of as-received material. 

cool and the repeat cycle will be called second cool. Repeatability tests were also performed on 
each sample type and with each scan procedure. 

Representative results from heat-first DSC tests are shown in Figure 3.15. The shape 
of the first-heat endotherms were ill-behaved and not repeatable, while the analogous second- 
heat endotherms were repeatable. Similar results with regard to the first-heat and second-heat 
endotherms are evident in the representative cool-first signatures shown in Figure 3.16. These 
effects are attributable to first-heat recovery of the mechanical deformation induced by shearing 
the small samples needed for DSC measurements from longer lengths. 

Note that both cycles of the heat-first test show two distinct exothermic peaks during the cool- 
down portion of the cycles; one peaking at approximately 39°C and attributable to the R-phase 
transformation, the other peaking at approximately 9°C due to the usual twinned martensitic 
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Figure 3.16 Cool-first DSC signatures for first- and second-cool of as-received material. 


transformation. The R-phase is a by-product of a partially annealed condition and does not 
adversely affect performance, as observed by Otsuka and Wayman [34]. Note that the cool-first 
signature does not show the martensitic transformation on the first-cool portion of the cycle, 
indicating that the alloy was delivered in a fully transformed state. From these DSC tests it can 
be seen that the transformation temperatures for the Nitinol in the as-received condition and in a 
free configuration are approximately /1 S =45°C (11 3°F), .4/=60°C (140°F), A/,=17°C (62.6°F), 
and Mf= 0°C (32°F). Note that the transformation start and finish temperatures indicated here 
mark the onset/termination of substantial phase transformation activity. Measurable activity can 
be seen in the tails of the exothermic and endothermic peaks of the DSC signatures. This 
observation will be revisited later in this section. The DSC results give significant insight into 
the state of the alloy, but mechanical properties are needed to model the behavior of the material. 
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Data pertaining to the material property and performance characteristics of the alloy were 
acquired through mechanical testing performed at Lockheed Martin Astronautics (LMA) [43] 
and at NASA LaRC. Although rigorous training could have been performed on the alloy to 
generate a thermomechanically stable actuator material, the elevated cure temperature of the 
glass-epoxy matrix material ( 176.7°C, 350°F) was considered to be high enough to risk degrading 
the expensive training. Thus, thermomechanical testing was performed to quantify the material 
characteristics as a function of thermal cycle. Discussion will focus on testing performed at 
Lockheed Martin first, followed by the testing performed at NASA LaRC. 

Isolength (constant strain) and isothermal tests were performed on nominal 0.254 m (10 
inch) lengths of the material. Isolength tests were performed to measure the maximum recovery' 
stress generated by the actuator material as a function of thermal cycle and initial strain. The 
results of this test are shown in Figure 3.17 for prestrain values of 2%, 3% and 4%. It can be 
seen that in all cases the recovery stress diminishes rapidly out to about 10 thermal cycles and 
approaches an asymptotic value at 50 cycles. Additional isolength tests were performed with 
the 2% and 4% prestrain samples after the thermal cycling tests to measure the variation of the 
stabilized recovery stress with temperature. These results are presented in Figure 3.18. Note 
that the recovery force starts to increase at a temperature very near ambient (~24°C, 75°F). This 
phenomenon appears to be in disagreement with the DSC results, but closer inspection reveals 
that it is not. The initial increase in the recovery force is associated with the transformation 
activity in the low -temperature “tail”, below the so-called .4.,, of the exothermic peak in Figures 
3.15 and 3.16. The abrupt change in the slope of the recovery stress vs. temperature curve, 
near 40°C (104°F), is associated with A s . Also note that the higher prestrain level generates 
greater recovery stress as expected. A couple of more subtle points can also be seen in this plot. 
Although the increased recovery stress at 4% prestrain does inhibit completion of the reverse 
transformation more than at 2%, as evidenced by the higher temperature as which it approaches 
an asymptotic value, the maximum slope of the 2% recovery stress vs. temperature curve is 
not as great as that for 4%. This trend is interesting in contrast to the data presented by Cross 
et al. [28], where reduction in the peak values of the recovery stress at lower prestrain levels 
is accompanied by higher maximum slope of the recovery stress vs. temperature curve. These 
insights will be revisited later in this chapter as well as in Chapter 8. 

A couple of options exist, with the available data from LMA, for modeling the behavior of 
the SMA material in the hybrid composite for comparison to experiment. The specimens could 
be thermally cycled enough (e.g., 50 times) to render the structure’s properties stable and simply 
model the SMA recovery stress behavior according to Figure 3.18. Alternatively, an estimate 
of the variation of recovery stress with thermal cycle number and temperature can be calculated 
from the available data. Construction of such a data set can be accomplished by normalizing the 
recovery stress versus temperature data to the corresponding maximum recovery stress versus 
cycle data set. The resulting carpet plots are shown in Figures 3.19 and 3.20 for the samples 
with 2% and 4% prestrain, respectively. Usage of these data sets allows modeling of the material 
system without any requirement on thermal cycling of the SMAHC specimens. 

Isothermal tests were conducted to measure load versus deflection on all of the previously 
described specimens, as well as a specimen in the as-received condition, at ambient and elevated 
temperature ( 120°C) to obtain estimates of the Young’s modulus. Representative plots, converted 
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Figure 3.17 Max recovery stress versus thermal cycle number 
for Nitinol ribbon samples with 2 %, 3 %, and 4% prestrain. 
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Figure 3.18 Recovery stress versus temperature for Nitinol 
ribbon samples with 2 % and 4% prestrain, after 50 cycles. 
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Figure 3.20 Recovery stress versus thermal cycle temperature for sample with 4% prestrain. 
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to units of stress versus strain, are shown for the cases of 0% (as-received) and 4% prestrain 
at ambient and elevated temperatures in Figures 3.21 and 3.22. The martensite variant 
re-orientation yield stress and initiation of linear elastic behavior of the martensitic phase are 
evident in the ambient temperature curves. Note the contrast between the as-received and 4% 
ambient temperature results. The 4% curve does not exhibit the long martensite variant re- 
orientation plateau because the prior 4% prestrain has diminished the available variants for 
re-orientation. However, the curves associated with the unique austenitic phase are linear and 
practically identical, as expected. 

Ambient temperature Young’s modulus estimates were calculated from only the initial 
(0-25.5 MPa, 0-3.7 ksi) portion of the ambient temperature curve, while the overall slope 
of the elevated temperature data was used for hot Young’s modulus estimates. Recall that the 
maximum normal strain is not expected to exceed 500 /<e in the configurations considered in 
this study. Therefore, evaluation of the ambient temperature modulus over the small range is 
considered appropriate. The resulting modulus data and values summarizing the isolength tests 
are shown in Table 3.7. 

Recall that in the SMAHC specimen fabrication procedures, described earlier in this chapter, 
the Nitinol ribbon was prestrained 4%. Although data was collected at LMA for this prestrain 
level, there is a fundamental inconsistency. The ribbon used in the fabrication procedures was 
necessarily released from tension after prestraining and before embedding in the composite, 
whereas the recovery force tests at LMA were performed immediately after prestraining and 
without regripping. Also, there was some suspicion that aging, after prestraining, may have 
some effect on the actuator performance. Therefore, additional isolength and isothermal tests 
were performed on excess ribbon with 4% prestrain remaining from the fabrication process. 

Isolength tests were performed on 0.127 m (5 inch) gage lengths to measure the recovery 
force versus temperature and thermal cycle. An adaptation of the prestraining apparatus was 
used to configure the Nitinol ribbon for measurement of recovery force. The apparatus, shown 
in Figure 3.23, allowed for convenient and repeatable installation of ribbon samples in the test 
machine. A picture of a sample installed in the test machine, after spacer removal, is shown in 
Figure 3.24. The ribbon was subjected to a slight preload ~4.45 N (1 lbf) to keep the sample 
taut prior to fixing the test machine grips (stroke control) and exposing the sample to thermal 
cycling. The ribbon was electrically heated to expedite the thermal cycling process. The recovery 
force and ribbon temperature at three locations were measured at 2-second intervals during each 
thermal cycle, each of which was completed in approximately 15 minutes. The raw data was then 
processed to estimate the recovery stress with a common temperature resolution and averaged 
over three samples. The resulting recovery stress versus temperature data corresponding to the 
second through fourth thermal cycle are shown in Figure 3.25. The data for the first thermal 
cycle are not included as that information is applicable only to the SMAHC cure cycle. The 
fiftieth cycle data is included for comparison with the LMA data. 

It can be seen by comparing the maximum recovery force with the data in Figure 3.17 that 
the performance of the Nitinol ribbon that was prestrained 4% and released is similar to that 
of the 3% prestrain data collected at LMA. Analogous curves were generated from the LMA 
data by averaging 2% and 4% prestrain data sets in Figure 3.18 and scaling the resulting load 
vs. temperature data by the maximum load for 3% prestrain at cycle numbers 2, 3, 4 and 50 
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Figure 3.21 Applied stress versus strain for 0.09x0.006x10 inch sample 
with 0% prestrain (as-received) at ambient temperature and 120°C. 
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Figure 3.22 Applied stress versus strain for 0.09x0.006x10 inch 
sample with 4% prestrain at ambient temperature and 120°C. 
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Table 3.7 Summary of Nitinol characterization tests. 

Maximum Stress, MPa Young's Modulus, GPa 
Sample Prestrain, % # of Cycles , 120 o c Ambient 



Figure 3.23 Apparatus for installation of a Nitinol ribbon sample 
in a tensile test machine for recovery force measurement. 
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Figure 3.24 Recovery force measurement in progress. 



Figure 3.25 Recovery stress versus temperature of the Nitinol ribbon 
(prestrained 4% and released) for thermal cycles 2-4 and 50. 

in Figure 3.17. The resulting curves are shown in Figure 3.26. Note that the 3% prestrain 
case “constructed” from the LMA data achieves roughly the same maximum recovery force, but 
over a much more narrow temperature range. This trend is a real physical phenomenon and is 
apparently attributable to aging effects. Also note the shape change with thermal cycle of the 
recovery stress vs. temperature curve, evident in Figure 3.25. This effect is not reflected in 
the data set constructed from the LMA data and expresses the importance of quantifying the 
recovery stress functionality with temperature vs. thermal cycle. 

Isothermal tests were conducted using 0.127 m (5 inch) gage lengths with the same gripping 
apparatus as shown previously. The entire assembly, including high-temperature MTS grips, 
was contained within a box furnace. Modulus tests were conducted on 3 ribbon samples at 
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Figure 3.26 Recovery stress versus temperature of the Nitinol ribbon for 3% 
prestrain constructed from the LMA data for thermal cycles 2-4 and 50. 

each of 6 temperatures; 21.1 1°C (70°F), 37.78°C (100°F), 65.56°C (150°F), 93.33°C (200°F), 
121. 1°C (250° F), and 148.9°C (300°F). Estimates from the three samples were averaged and the 
resulting modulus versus temperature behavior is shown in Figure 3.27. Stress-strain behavior 
similar to that shown in Figure 3.22 was observed, but the modulus estimates are somewhat 
different. These differences are also probably attributable to aging effects, akin to the rubber- 
like pseudoelasticity behavior. 

Note that the temperature range over which the modulus exhibits the most variation in 
Figure 3.27 is more narrow than the .4,-4/ range interpreted from Figure 3.25. This obervation 
is consistent with the data collected by Cross et al. [28]. The low-temperature transition 
point is associated with the temperature at which there is significant austenite and high-enough 


56 




Figure 3.27 Modulus for the Nitinol ribbon prestrained 
4 % and released as a function of temperature. 


temperature to not allow immediate transformation to martensite upon stress application. The 
high-temperature transition point corresponds to the temperature at which the austenite volume 
fraction begins to dominate. 

Published data will be used for Poisson’s ratio (0.3) and the transverse direction CTE for 
the martensitic 6.6e-6 /°C (3.67e-6 /°F) and austenitic 1 le-6 / C (6.11 e-6 / F) as measure- 
ment of these properties as a function of temperature would require significant additional ef- 
fort/development. In light of the previous discussions concerning the differences between the 
data collected at LaRC and LMA, the data collected at LaRC are more representative of the 
thermomechanical properties of the Nitinol ribbon actually embedded in the composite and is 
more accurate. Therefore, the properties interpreted from this data were used in making the 
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predictions shown later in this document. These properties are tabulated in Table 3.8. Note that 
the modulus data has been linearly interpolated within the more coarse temperature resolution. 
The 2-direction CTE has been interpolated between two transition points noted in regard to Fig- 
ure 3.27, 43.33°C (1 10°F) and 82.22°C (180°F), to approximate the variation within the reverse 
transformation temperature range. The “in-service” austenite start and finish temperatures, as 
interpreted from the recovery stress versus temperature data, are 26.7°C (80°F) and 148. 9°C 
(300°F), respectively. However, the temperature range of property variation is better modeled 
by the temperatures at which the recovery stress vs. temperature curve rapidly changes slope, 
as discussed previously. 

In this chapter, the glass-epoxy and Nitinol ribbon material system selected for this study 
was described. Discussion of the procedures for fabricating the SMAHC beam specimens was 
followed by results from thermomechanical characterization tests on the constituent materials. A 
full set of temperature dependent material properties were generated to quantify the constitutive 
model developed in the previous chapter. This constitutive model can now be incorporated in 
a thermomechanical model to predict the behavior of structures fabricated from this material 
system. One such model is developed using a finite element approach in the following chapter. 


58 



Table 3.8 Summary of the thermomechanical properties of the Nitinol ribbon. 



Recovery Stress 



Temperature 

Cycle 2 

Cycle 3 

Cycle 4 

E 

O 2 

°C 

MPa 

MPa 

MPa 

GPa 


21.1 

0 

0 

0 

27.17 

6.606 

26.7 

2.130 

1.276 

2.130 

24.82 

6.606 

32.2 

7.233 

5.530 

7.233 

22.41 

6.606 

37.8 

16.18 

13.20 

14.05 

20.06 

6.606 

43.3 

31.92 

28.95 

28.52 

25.72 

6.606 

48.9 

70.65 

70.23 

78.74 

31.37 

7.236 

54.4 

135.4 

137.9 

151.9 

36.96 

7.866 

60.0 

192.0 

199.6 

217.9 

42.61 

8.496 

65.6 

240.5 

256.6 

269.4 

48.27 

9.108 

71.1 

297.1 



54.88 

9.738 

76.7 

341.8 

341.8 

341.3 

61.43 

10.37 

82.2 

374.5 

373.7 

372.0 

64.19 

11.00 

87.8 

402.2 

399.2 

399.2 

63.16 

1 1 .00 

93.3 

425.6 

419.2 

418.8 

62.06 

11.00 

98.9 

449.9 

438.8 

437.5 

63.92 

11.00 

104.4 

469.0 

451.6 

448.2 

65.78 

11.00 

110.0 

481.0 

464.8 

458.8 

67.64 

11.00 

115.6 

491.6 

475.8 

467.3 

69.50 

11.00 

121.1 

500.1 

485.6 

472.9 

71.36 

11.00 

126.7 

508.6 

490.3 

478.8 

70.81 

11.00 

132.2 

515.0 

496.7 

481.8 

70.33 

11.00 

137.8 

520.5 

500.5 

485.2 

69.78 

11.00 

143.3 

523.5 

502.2 

487.3 

69.29 

11.00 

148.9 

527.8 

503.9 

490.7 

68.74 

11.00 
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Chapter 4 

Finite Element Formulation 


Previous models have been limited in their ability to predict the behavior of SMAHC 
structures because of inaccurate/impractical constitutive models and/or classical formulation 
approach. The present formlulation incorporates the ECTEM and will be developed using a finite 
element approach. This combination will later be shown to accurately capture the mechanics of 
the system and to afford the generality/versatility to allow comparison with experiment. 

The governing equations for the geometrically nonlinear, static thermal response and the 
geometrically linear dynamic response of a rectangular panel-type structure subjected to com- 
bined thermal and mechanical loads is derived in this chapter. It is necessary to account for 
geometrically nonlinear out-of-plane deflection because of the potential for large static deflec- 
tions due to thermal buckling. Therefore, a finite element with cubic out-of-plane interpolation 
functions is deemed necessary. 


4.1 Element Displacement Functions 

A rectangular 24 degree-of-freedom (dof) plate element [44] was chosen for the subsequent 
finite element formulation. The element has four nodes with two in-plane and four out-of-plane 
dof per node, and is shown schematically in Figure 4.1. The in-plane displacement functions 
are bilinear and that for the out-of-plane displacement is bi-cubic: 

»(.r, y) = a\ + «>x + a.\y + a+xy ^ 

v(x. y) = dr, + «<;■*' + a-y + 


and 


w(x, y) = ci + c 2 .r + c.yy + c 4 .r + cr,xy + cay 2 + c r x A + c 8 .r“// + cyxy- + no;/ 3 
+ cn x'*y + ci>.v 2 y 2 + envy 2 + ci 4 .r' 5 //“ + + no^'V* 


(4.2) 


where the a, and c, are generalized coordinates. 

The displacement functions can also be written in a form particularly suited to numerical 
implementation in terms of the nodal degrees of freedom and appropriate bilinear and bi-cubic 
interpolation functions of the element coordinates : 


4 1 

V :(.(',;</) = L i V)«* = -f[(« ~ J ’)( 6 “ V ) u ' + • r ( / ' _ 

f=i U ’ (4.3) 

i j 

v(.v.y) = ^ Li{£. y)v, = —[(a - x)(b - y)v \ + x{b - y)v> + xyv . :$ + (a - x)yv 4 ] 

j=l 
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and 


21 ) 


tr(.r. n) = ) Li(:r, ij) q, 


i=4 




(«•'* + 2.r :t - :to.r'“)(/> 3 + 2// :1 - 3/>r)«’l + (3fl.r 2 - 2.r :i )(/> ! + 2//' J - 3/»/-)«’.» 


+ (3o.r‘ - •>.)•■*) (:%' - 2»/ ; * )•«’;* + (r/' ! + 2,r 3 - 3rt.r) (3/«/ J - 2y' ! )«M 

.)•(.(• - r/) 2 (/) 3 + 2y/‘* - 3 by-) ^ + «(■'■' - «•»•") (/>'* + 2// ! - 3/>.i/') ^-r— 

+ a(.v' ] — a.r 2 ) (Miy 2 — 2i/ 3 


+ « 


flit’ 


+ h(a i + 2.r 5 - 3 ax 2 )y(i/ - />) 


<7.r 7 .s 
> /<7u 

<% 


-r </.r( .)■ - o)' (:!/»/- - 2i/’) (— ) | 

1 ) 


rler 


+ /)(3«-r* — 2.r f ) (y by , ^ j ^ 

+ abvifix - of(y - , + < lh - r y(-' J ~ «•'■)(•'/ - 


| + /)(o :! + 2.r 3 - :k/.r J ) (i/ 3 - (^) 4 

>/c) 2 ic\ . . > , .>/() 2 w 

I'dij' 

+ - **> (|£), + - *»)(«;), 


(4.4) 


where u, and v, are the in-plane dof, q, denotes the out-of-plane dof («m — + ( 7777) )* ar| d a ^ } 
are the element length/width, respectively. It will be advantageous for discussions that follow 
to cast these equations in matrix form 



[L, t \ 0 

LM 0 

0 [L w \ 



] { o m b } 


(4.5) 


where the nodal degrees of freedom have been collected into membrane {«„,} and bending {a&} 
components 


L«mJ = L«1 11 ~ U -i U 4 i’l v 2 <’3 v 4 J ^ 

[ttf, J = [u»i ... W4 W, r \ ... It) ,*4 W,pi ... 10, y 4 W, X y 1 ••• ™rry 4 J 

Also, the displacement shape functions have been collected into vectors 

[L u \ = [L] L> L-i U 0 0 0 OJ 

[/,,,] = [0 0 0 0 L\ L> L :i L 4 J (4.7) 

[L w \ — |_i 5 ^6 • • • L>\) J 

Recall the von Karman strain displacement relations, equation (2.32). The membrane 
midplane strains and midplane curvatures can be written in terms of the nodal degrees of freedom 
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and interpolation functions by substituting from equations (4.5) and (4.7): 


a -j 


w;„} = { 
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\ u 

U'-rr 
ir 
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L f'd 'V J "f 

[l w ,,x\ 

U «■•»*] } {«;>} = [0]{«i} 


{<1 in } — [f in ] { (l in } 


(4.8) 
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The bending midplane strain requires some manipulation to obtain a useful form. One approach 
follows from recognition that the second-order (nonlinear) terms can be expressed as the product 
of first-order terms: 


It will be useful later to make the definition 
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(4.12) 


The vector { } can be written in terms of nodal degrees of freedom as before in equations 
(4.8) to give 


where 


{6-} = 


[<*] = 




{«/>} = [0]{«i} 


l- r ,-x l-li-x ■■■ 1-2(1. x 

. L, 


_l- 2 >.y 1 -li-y ■■■ i- 20 <y 

Thus, the (nonlinear) bending midplane strain can be written as follows 

K) = 


(4.13) 


(4.14) 


(4.15) 
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where [0] is a constant matrix whose entries are evaluated based upon the nodal degrees of 
freedom from the previous iteration in an iterative procedure. Such a procedure will be described 
in Chapter 5. 




i+i — 


0 


0 




(4.16) 


The form for the von Karman midplane strain and curvature can be summarized as follows. 


{«} = {<“}+-{*} 

= K,} + K} + ~~M (4.17) 

= b]{«b} 

Finally, it will be useful in subsequent discussions to write the resultant force components from 
equation (2.38) in terms of the nodal degrees of freedom: 

{A',,,} = [-4][r,„]K,} 

)A’,} = i[.4]|«l[r.Hn} < 4I8 > 

|.v 4 ) = {»][<■*]{<**} 


4.2 Euler-Lagrange Equations 

The equations of motion for a mechanical system can be derived following a Newtonian 
approach or a variational approach. Variational methods have proven to be very effective tools 
in deriving the equations of motion for complex mechanical systems. The governing equations 
for the mechanical system considered in this study will be derived using the extended Hamilton’s 
principle, as described in Meirovitch [45] for example. Hamilton’s principle states that the path 
satisfying Newton’s law extremizes the time integral of the Lagrangian with respect to all possible 
paths between the limits of the time interval defining prescribed initial and final configurations. 
The extended Hamilton’s principle simply extends this idea to systems with non-conservative 
forces: 

( 2 

b f (T - V + \Y)fH = 0 (4.19) 

h 

and 

T= I 

I) 

V = J vdD (4.20) 

D 

IT = J {b} r {u s -„ y }dD+ j {t} r {us- u g }JS 
L) 's 
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where f> is the variational operator, 7 is the kinetic energy, 1 is the elastic potential (strain) 
energy, IT is the work done by non-conservative forces, p is the mass density of the structure, v 
is the strain energy density, {It} is the body force vector, {/} is the surface traction vector, { ».,} 
is the displacement vector of the structure, {«,,} is the displacement vector of the grip/boundary 
(assumed uniform over the boundary to allow for base excitation of the structure through a 
fixture), H indicates the domain (volume) of the structure, and .s’ indicates the surface of the 
structure. The displacement of the structure relative to the boundary can be expressed as 

{(/) = {».,- I/,,} (4.21) 


such that the energy terms can be rewritten in terms of the relative displacements: 


/' = 
r = 


J 7 } p{'< + 11 y) 1 {■» + ii g }dD 

D 

j vdD 
'[) 


ir 


I {b} 1 {u}dl) + j {l} r {u \(is 

l) s 


(4.22) 


Consideration is restricted to a panel-type structure subjected to a steady-state thermal 
excitation and mechanical excitation consisting of surface pressure, out-of-plane concentrated 
forces, and/or inertial loading. The following assumptions are made for the development of the 
governing equations: 

1. Body forces other than inertia are negligible. 

2. The structure is in a state of plane stress. 

3. von Karman deformation relations are valid. 

4. Linear elastic material properties are valid. 

Assumptions 1 and 4, and the above-mentioned restrictions allow the potential energy and non- 
conservative work terms to be simplified as follows. 


'= = J'p.ilWhiD 

D D 

Nf (4-23) 

ir = y [pic + Y f-w H {.(■} - {*},■)]</>' 


where {a} and { (} Aj are the stress and mechanical strain vectors, respectively, w is the out- 
of-plane displacement of the structure, p is the surface pressure, /' are concentrated forces 
applied normal to the plane of the structure, Nj is the number of concentrated forces, { .?■ } is 
the coordinate vector in the plane of the structure, {j }, designates the point(s) of concentrated 
force application, and <’>({.<•} — {•* },) is a spatial Dirac delta function defined by 


*({■«•} = 0 for (4.24) 
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Assumptions 2 and 3 allow the constitutive relations for a single lamina to be written as shown 
in the previous chapter, equation (2.26): 

Ma = [<?]*({«}* - M Tk) 

\ .. (4.25) 

= [0] A- ( { e ° } + ~{ / '} — j {°( T )}k ( M ) 

f 0 

Recall that the mechanical strain is equal to the difference between the total strain and the 
thermal strain ({f} \/ = { * } — {<}y). Substituting this expression into the form for the strain 
energy in equation (4.23) gives 


y = J j ^(Wa ~ _ { ( }rk) (lI) 

S -7 (4.26) 

h_ 

= j j f ) r _ -Ma [0] t { f }rt ” Mn [Q]k-{ ( hk) (,I) 

S ItL 
> 

Commuting the variation and integration operations in equation (4.19) results in the following 
expression for Hamilton’s principle. 

f_> 

f (W - t>V + 6\Y)df = 0 (4.27) 

h 

Consider the first variation of the kinetic energy, from equations (4.22): 

ST = S j \p{i< + u g } T {u + v a }dD (4.28) 

I) 

Commuting the variation and integration operations again and integrating through the thickness 
results in the expression 

ST = j ^({6«} r {u + ity} + {v + it g } 1 {Su}}dD 

[) 

= J p{Su} T {it + u g }dD (4.29) 

I) 

= J pit {bit} 7 { <i + ii g }dS 

s 

where f>u g has been omitted because the instantaneous constraints are assumed to be held constant 
during application of a virtual displacement to achieve the varied path. Similarly, the first 
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variation of the strain energy in equation (4.26) can be written as 

h_ 

= (S j j ^(i ( _ ^ [^M f hi ~ W 71 hi) 


M 


(II) 


S _ II 


I j It h- + { f }t M k\^ ( )k ~ “{^ }* [0] A {< }rk) 


(II) (4.30) 


.s' _!l 


I j (j^lt [^<?] / { ( ) a _ {Mif [QLMvi) 


(ID 


.s' _ !l 


Incorporating the von Karman strain-displacement relations, equations (2.32), gives 

7 7 

M ' = / / [()^ [, ) 7 +c)Hi'')[0],.({< ll )+.-{^))-({^ O } I +--{*/.l Y )[0] 1 . f {ofrDjrfr 


.S' 


(ID 

(4.31) 


Expanding the integrand results in 


.s’ -1 


J 1 

{^°} 7 [Q] a . j {tt{T)} k dT - z{Sk } 7 [Q] a I {tt(T)} k dT (ID 


Jo 


7o 


(4.32) 

and integrating through the thickness generates the stress resultant notation from equations (2.35). 
The following form for the first variation of the strain energy results. 


(V= f + {*<") 7 '[B]M + M T [B|(<‘') + 

.s 1 

- {^+"} 7 {A'a/ } - {a) r )A/ A 7'|)(/.S’ 


(4.33) 


Finally, referring to equation (4.23) the first variation of the non-conservative work has the form 


Oil 


' = *> f [/>“’ + 51 f'i w 
i «'=» 

N.t 

«’ + 51 f'* Sw “ {•*■}» )] (IS 

(=1 


(4.34) 
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Recall that <S({.r} - {.r},) is a spatial Dirac delta function and shouldn’t be confused with the 
variational operator. 

The form of Hamilton's principle given in equation (4.27) can now be rewritten by substi- 
tuting from equations (4.29), (4.33), and (4.34): 


/[/ + <<»} - [«]W - IM'[«H<") - fM r [»](-'> 

/i is 

+ {^°} 7 {.Va/ } + { ^ /■• } 7 { 3 / A /■ } + pb w + f'i *«’*({ •' } - 

(=1 

The first term can be integrated by parts to give 


dt = 0 


(4.35) 


/ 


ph{b »} 7 | {n + 

S .*>' 


I> 




- {bn) T [B]{t Q } - {6 k} 1 ' [/?]{/.-} + {*< } {-Vat } 

.v, 

+ {/>k } 7 { M&t} + P^ tr + } — {- r I* ) } ^ ,s 


(It = 0 
(4.36) 


Note that the first term in the latter equation is equal to zero because the integrand vanishes as 
the virtual displacements { <*>(/} are defined to be zero at times l\ and /_*. 

The expression given in equation (4.36) resulted from application of extended Hamilton s 
principle to the energy terms associated with the entire structure. The same approach is applicable 
to any portion of the structure, e.g., finite element, due to physical constraints. Therefore, the 
following development will focus on the element level. Recall that expressions for the midplane 
membrane strains { t 0 } and curvatures { k } were given in terms of element degrees of freedom in 
equations (4.17). The variation of the midplane strain within the element would then be written as 


{«<“} = [('„]{*«„} + T , MK «]<“»} + (4.37) 

or 1 ! 

{Or"} = [C m ){6a m ) + ^6}{G} + ~[9]{bG} (4.38) 

by substituting from equation (4.13). The components of the second and third terms can be 
examined as follows 


-[(>()][( + -[0][f 


i 
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Sw. T n\j 1 
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U\ T bu\ t 'j 

< Sw,y (('■;/ | 

) +-< 

u\ v bu\y | 

1 C. r +fnc. r J 


[u\ y fnv,r +U\ r b-Uhy ) 


(4.39) 


69 



It can be seen that 


= W\\h(;\ ' (4.40) 

= [«][( «]{*«»} 


So, the final forms for the membrane strain and curvature vectors in terms of element degrees 
of freedom are as follows. 


{^ f °} — + [$][f V;] { } 

{*>*} = [G>] {(>«/,} 


The structural displacement-related quantities in equation (4.36) can be written in terms of the 
element degrees of freedom from equation (4.5): 


{M 7 {" + ''</} 


r t 

{^b/ife} [f'UiNc] [I- (U'lrJ { (hub ~t~ My } 


dai, 

tui m 

da b 


{La} {L v } 0 

o 0 {/-«■} 


[Lu\ 0 

LGJ o 
J L 0 U-J J 


«»» o 

(<b il<i 


Tr 


{L u }{L u y +{L„}{L r }‘ 


0 


('m 0 

(>b i’y 


0 {L W }{L W }‘ 

{t>«,„} 7 ({L U }{L U } T + + + 


and 


tnc = {L w } 7 {<^ 4 } = {ba b } T {L w } 


(4.42) 

(4.43) 


where it has been assumed that the base acceleration component u g is normal to the plane of 
the structure. The expressions in equations (4.41), (4.42), and (4.43) can be substituted into 
Hamilton's principle, equation (4.36), to give 


/[/ + {L,.){L„i r ){a 


>} + {hab} { L u , } { L w } {('ii + 


r}) 


A, 


+ + [#] [Q] { <^4 } ) 7 [ A ] ( [f ] { (I rn } + ^[0][("«l] {«(,}) 

+ ([(',„]{£«„,] + {V}[('v}{b(lb}) T [B][(-'b]\ (' bl + {^aft} 7 [f 6 ] / [^](K , to]{«to} + ^][(’’o] {«/>}) 
+ [G] 7 [«]K "fc]{"*} ~ ( [('raH^m } + [0][Q]{^«6}) 7 {Aa7 } ~ {^«b} T [('b] 1 {Mat} 


— pSir 


ilf 


£./>■ - {•*■},■) }*lAt 


J=] 


(It = 0 

(4.44) 


where the spatial integration has been limited to the area A ( of the finite element in concern and 
the number of concentrated forces in the summation has been limited to that number u f within 
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the element area. Expanding the integrand results in the following expression. 


1 1 


+ {^/ w } T [r,„] ;r [,i][r' rM ]{« w } + ^a m } T [( \ B ] 7 MMO]{« 6 } 

+ {bdb} 1 [C 'h] 1 [^] 7 [^]K tn]{ fl m} + — } K tf] [$] [-4][^]K o]{ f, fc} 

+ {b<hn} T {C ni ] 1 [B}[(' h \{a,,} + {iu b } T [C e ] T [0] 7 [B][C b ]{a b } + {Sot,} 1 [O/Wm ■]{«».} 

+ \{fia b } r [G,] T WHO]W} + {^b} T [(' h ] T [D}[(),}{a h } - {^„,} 7 '[G„] 

- [0f{X±T} - {t<>b} r [Cbf{U±T} 

(11 = 0 


”/ 1 
{f>a b } l {L w }p- {f>a b } T {/.„,} ^2 fair 6({x} - 

?= 1 


(4.45) 

Consider the product [#] 7 { A }. Expanding the product in terms of the matrix components 
results in the following identity. 


[oy 




U’.r 0 


0 ir.„ w 


y “ J 



' A x A, # 

A XI/ A y J 


■::} 


(4.46) 


= [A r ] KV]{«6> 

This identity is useful for further manipulation of equation (4.45). Note that the terms 

{^«/,} 7 KV] r [^] T [.4][C r „]{« m } and {Sat} 1 [('o] 1 [9] 1 [B][Cb\{ab] (4-47) 

in equation (4.45) do not have symmetric counterparts. However, noting that from equations 
(4.18) [.4][C m ]{a„,} = { A,,,}, the former of the two terms in equation (4.47) can be split as 

follows 


{^6} 7 'Ktf] 3 ’[0] T [.4][r TO ]{« Mt } = ^{^«6} r KV] 7 ^] r [.4][c m ]{a,„} 


+ 4{<^«6} 7 KV] r [0] r {iV w } 


(4.48) 


and the identity in equation (4.46) can be employed to give 

)<s<i t ) 7 [r s ] r [«] 7 [/i][o,, ]{«,,,) = i{n« t } r [c»] r [«] I '[.4]!<%„ ])«.»> 




(4.49) 


Now, the first half of the term has a symmetric counterpart in equation (4.45) and the second 
half is symmetric itself. Similarly, noting that [#][0]{ rt f>} = {Aft}, the latter of the two terms 
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in equation (4.47) can be split as follows 


(M,) 7 [r B ]'' (<)]■' [B][a [r # ] 7 [«] 7 '[«](r t ]{»4 

, (4.50) 

+ ;{*«* 

and the identity in equation (4.46) can be employed to give 

I . (4-51) 

So, both terms in equation (4.47) have been transformed into symmetric forms. Finally, the 
identity in equation (4.46) can be used to rewrite another term in equation (4.45) as follows. 

[flJ^A'AV } = {< i, «fe} 7 [("«] 7 [ V A /][r w ]{o fc } (4.52) 


Substituting the expressions from equations (4.49), (4.51), and (4.52) into Hamilton’s principle 
for the element, equation (4.45), and separating terms involving {ba,„ } ! and {hub} results in 
the following expression. 


f If { [/</>({ i. ){/..}' + ) 7 ){S,„) + [(•„,]■' [,1][C, „](«,„} 

/i A, 

+ i[c„,]'[.il(»][r,Hn} + K',„] 7 W,. ]{«,,} - |c„,] , {.\a/')1 

+ {(’"(.) ' {/,„.} { i „,} 7 {( if , + ( i , ; [ + -[( ',] 7 [ 0] 1 [4][( } + — [( \ i '/ [ A ,„ |[( 1 

+ 5|r,) T w r I.'tlWlOH«*» + K»] r tB]K-K«-.) + |r,) r M r W,]K) 

+ + |[(y 7 [«]['’][c.H«»} + [C S ] 7 (/)][(',](«,} - [r»] r i.Y Ar ][f„K„ ( ,} 


»/ 


- Vtfi'Ur} - < - {/.„.} 53 *({.<■} - {.<}, )]}, m , 


(If = 0 


»=i 


(4.53) 


The virtual displacements {<)«,„ } and { } are arbitrary and independent, so the corresponding 
components of the integrand in equation (4.53) must each equal zero. The two resulting equations 
can be written in matrix form: 
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Pb(l) 

0 


+ 


nl h + {nib) n »h,„ 

II 1 , n b 0 



+ 


fbi f ) 

0 


+ 


PbAT 

,PmAT 


m h 

0 


0 

m m 



(4.54) 
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[m w ] = j ph( < {Lii}{Ltt} r + {L v }{Lr} 
A, 

[^*r»] — I [( m] [*^]K 

A, 

mb\ 

[»!*]= I K »] T [A m] [( '»]-'- V 


where the element matrices are as follows. 

M= f pl l {l.„.}{l.„.} T <IA, 

A, 

M = f ia] T '[«][(y,/.4, 

[/■(,„, ] = f ((■ ( ] 7 '[«][C,„]rf/l, 

A,, 

[A-a/] = j [0] 7 ’[A>][0R4< 

A. A ' 

[(»lfe)/i] = j ([O] 7 [^] 7 [//][( ft] + V v] [Aft]K #] + K ft] [AA][^][f ^]) < A.l( 
A, 

= f [r»i r [«] , '[.4][f„K4 r 

A, 

[«2»] = '-,[[< '»] J [*] r [.-l][»]K 

A, 

{Pb(t)} = j p{-r.y.i){L w }(lA ( 

A, 

{/>fcA7’} = J [Cb] 7 {M*r}d Ac 


[l } 1 mb. — [ /; 1 but ] 


= £ //'(/) H{-r] - {.v},){L t A<IA ( 
•=' A. 

{PwAt} = J [Cb ] 1 


(4.55) 

The element area integrations can be conveniently and efficiently evaluated though Gauss- 
Legendre numerical integration, as demonstrated in Appendix B. The governing equations for 
each element can be assembled to form the Euler-Lagrange equations for the entire system. 


Mb 0 

0 M„ 


A b 



+ 


l 

+ T 
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(4.56) 


fA(0 


I J F b(t) \ , \ Aa7 
10 1 0 \{PmAl 





A/ft 0 
0 A A, 



or 


A} + ([A'] - [A at] + r } [ A'l] + ^r[A2]){.4} = {P(t)} + {F(f)} + { Pat} - WWti 1 )} 


(4.57) 


where [A/] and [A] are the usual system mass and linear stiffness matrices; [A at] is the 
geometric stiffness matrix due to the thermal in-plane force vector { A A7'} ; [A l] and [A 2] are 
the first- and second-order nonlinear stiffness matrices which depend linearly and quadratically 
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upon displacement {.1}, respectively; {P(J)} is the surface pressure excitation load vector, 
{/•’(/)} is the concentrated force excitation load vector, { Pat] is the thermal force vector, and 
{/ ,,(/)} is the inertial excitation base acceleration vector. The subscripts h and w denote bending 
and membrane components, respectively, and the subscript B indicates that the corresponding 
stiffness matrix is due to the laminate bending-membrane coupling stiffness matrix [B], Note 
that the stiffness matrices [A - ], [A'a7']» [A 1], [.Y2] and the thermal force vector { P^j } are all 
temperature dependent. 

4.2.1 Separation of Static and Dynamic Components 

Recall that the system is restricted to be thermally steady-state. It is assumed that the structure 
will undergo small dynamic motions about a stable static shape (flat or thermally buckled) so 
that the system degrees of freedom will consist of static and dynamic components 

{.1} = {!),+ {•!),, (4.58) 


such that 

{.!),, «(.)). {,!}; so (4.59) 

Consider the effect, on the element level, of equation (4.58) on the form of [6 1 ] in equation (4.12): 

(4.60) 

It can be seen from equations (4.55) that the nonlinear stiffness matrices, [n 1] and [//2], are thus 
affected. For example. 
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ir. H , x m 


U'd«y U'thJc _ 


["4ml = 


Similarly, 


= / [(«]>]' 
s 

= J [oi i '(I*]J' + [*]3')[^]K»k‘.' 

= / [f -. 1 ' [<’]ri--'][r,„]rfS' + 1 Micyrf.s 

.s’ S 

= ['" 4m ] s "F ["■ 4m] ( / 

= T ["lmfc],/ 


(4.61) 


(4.62) 


Note that the resultant forces due to membrane and bending strain, equations (4.18), are also 
affected. So, these forces can be written as 


or 


{A,,,} = {A/,,}, + {Am},/ and {A/,} — {A/,}, + {A/,} (/ 

[Am] = [Am], + [A ,„],/ and [A/,] = [A/,], + [A /,] ( / 


(4.63) 

(4.64) 
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Therefore, 


(4.66) 


[»4] = [nib],, + and [(»!/,)«] = [("U)/j]., + [("4)«] ( / ( 4 - 65 ) 

Finally, the second-order nonlinear element stiffness matrix takes the form 

[«2»] = \j [C,] T [*] T W[*][f,]</.‘.- 
s 

= * [ K'.fuC +[«]'/ PM, + [»],! )K«K S ' 

s 

= ^ j ([c#] 7 [o]1 [4][0],[( V] + [( 0 ] [o}^ [- 4 ][0],/K #] 

+ K^] 7 [0]t [-4][^],KV] + [O] 7 [0}!, MM/KV])^ 

= [»-&]., + [«2*]w + [«*L + [n'hh 

This distributive property conveys to the system matrices: 

[•VI] = [V 1]., + [A'l],/ (4 57) 

[-V2] = [-V2], + [V2] s(/ + [.V2] (/s + [A2] rf 

and the system equations of motion, equation (4.57), become 

[M]{A} d + ([A] - [I<at])({A}, + {A} d ) + :j(|-Vl]„ + [ VI ],/)({■!}, + {A},,) 

+ \([ V2] s + [.V2], (/ + [V2] ( ,, + [•'V2] (/ )( {.4} s + {A},,) (4 ‘ 68) 

= W)} + {FV)} + {P±t} -[MW r g} 

It is desirable to separate the equations governing the static and dynamic response. However, 
there are terms which couple the static and dynamic displacements. The term [N2] d can be 
neglected because it is second order with respect to {,4} (/ . The form of other coupling terms 
can be clarified as follows. 

Consider the product of the first-order nonlinear dynamic stiffness and the static displacement 
vector. On the element level, referring to equations (4.55), the four components of this matrix 
can be manipulated as follows in discussions designated by numbers 1) through 4). 

1) The product of the first-order dynamic bending stiffness and static bending nodal degrees of 

freedom has the form 


[»H] ( /{«fc}« 






(4.69) 


This integrand can be rewritten by noting the following 

[A ro ] (l [( '#] { U /) } ,s — [^L {-V/"},/ 

= [01 [-4][( mli' 7 ™}, 


(4.70) 
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so the expression given by equation (4.69) can be rewritten as 


m,/], y (4.7D 

s 

The form of the integral in equation (4.71) can be immediately recognized as the product of 
the first-order static coupling stiffness and the dynamic membrane displacements. Therefore, the 
expression in equation (4.7 1 ) can be written as follows. 

[»UL/ { ( 'b}, = ["UmU"™},/ (4.72) 

2a) The product of the first-order dynamic bending stiffness due to the [B] matrix and static 
bending nodal degrees of freedom has the form 

[<»U)«], ,<«»}„ = j ([('»]' [»]'! W») + [f/[.v»ur»] + [(»]' (4.73) 

■S’ 

but recall from equations (4.18) and (4.46) that the substitutions 

[er»H"4,= ).V l >L and [<{.V,,1, = [.V s ],[C«]{«i),, (4.74) 

can be made. It follows that the expression given in equation (4.73) can be rewritten as 

[{»h)n}A a b}s = I W} ( / + j [('a} 1 [- v /»],/[( 

S f (4.75) 

+ J [c.fmouc^M, 

s 

2b) Conversely, the product of the first-order static bending stiffness due to the [B] matrix and 
dynamic bending nodal degrees of freedom has the form 

[i«4)bU« 4 },, = f ([f.iWlBHO] + ir«f [.\y,[c»] + [oftBMJo] <4.76) 

.s' 

but equations (4.18) and (4.46) allow the substitutions 

WdWI,,= )AU,) and MltV,)., = [ V.],,K-V] Wl, (4.77) 

It follows that equation (4.76) can be rewritten as 

' s ' 5 (4.78) 


76 



The first two terms of equations (4.75) and (4.78) are obviously identical, but the third term 
requires further investigation. Recall from equation (4.13) that [#][( #]{«?>} = Therefore, 

term three in equation (4.75) can be rewritten in terms of out-of-plane displacement derivati\es. 

!<•,/.„ t/'s-.y \ (4.79) 

K'.w + 


~Wd*r 

0 


1 

0 

"'(/•V 

U-J 

' 1 

Jl'thy 



l 


It is clear that the previous expression is the same as the following. 

[WH«6}rf = [»U«}rf = 


M’.S ^ 


\ U's'f ll \hr ) 


1 U'd-I 1 

1 V 


{ ( = 

M s- // l i (h y / 


l »’</•</ J 

[ U's.y U\hx + W* *.r ) 


(4.80) 


Thus, the third term in equations (4.75) and (4.78) are identical, and the final result gives the 
following equality. 


[(" 4)/?], /{«*}., = [(" U,)n} s { u b}, 


(4.81) 


3) The product of the first-order dynamic coupling stiffness and the static membrane displace- 
ments has the form 


s — j [('h] 1 [Q)d 


[» 4m], /{«»»}., = 

Equations (4.18) and (4.46) allow the substitutions 

[.4][C„, ]{«„,}., = and [0]] {A w }, = 

So, equation (4.82) can be rewritten as 


(4.82) 


(4.83) 


(4.84) 


The form of the integral in equation (4.84) is easily recognized as the product of the first-order 
static bending stiffness and the dynamic bending nodal degrees of freedom. Thus, the final 
equality becomes 

[»4,„ !,,{»,»), = [»4],M,, < 4 - 85 > 

4) The product of the first-order dynamic coupling stiffness and the static bending nodal degrees 
of freedom has the form 

(4.86) 






but recall from equation (4.13) that [#][( V]{ r 'fc} — [^]{4 }• So, the analogous product in equation 
(4.86) can be written in terms of out-of-plane displacement derivatives: 



«t’rf,x 0 

ft,- ) 1 

J «',/-x W s . x | 

= m«). = 

0 «•,/•-/ 


, "',/•« j 


Jl'd U'd'T _ 


«'s,r +«V-x "’s -y ) 


(4.87) 
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This expression can be manipulated into the following equivalent form 




0 


«-W irj-x ) 

[V*]{«ft} ( / = l°}A ( '}j = 

0 

Ji-s.fi 

Ws-y 

M- 

l ) 

{ tt's-y "V-r +«W IVJ.y ) 


which provides the final needed equality. 


1 mft] ( /{ ft ft}. s — [ n i mft]. s { f, ft },/ 


(4.88) 


(4.89) 


The four equalities given in equations (4.72), (4.81), (4.85), and (4.89) can be written in the 
matrix form 


^ I b.il “t” ( ^ i ft ) B.il ti 1 bin, (I 
H^mb.il 0 


«b.s 

(I III .» 


» lft,.s + (»ift)tf.* 11 1 bin .s 
» 1 mb.* 0 



(4.90) 


from which it is easily seen that the assembled system matrices share the analogous equality: 


[-VILML = 


(4.91) 


The cross-terms generated by the expansion of the second-order nonlinear stiffness into static 
and dynamic components need manipulation as well. The product of each cross term with the 
corresponding displacement vector, on the element level, can be manipulated as follows in the 
discussions designated by numbers l') through 4'). 

1') The product of the second-order bending stiffness cross-term [n2<,] srf and static bending nodal 
degrees of freedom {<//,}, is written as 

- / KtrwTMWitrvKswj, (4.92) 

.s' 

but recall from the previous manipulations designated by 2b) and 4) that 

[0]»K'}j = =* [^L/K’tflU'ftL = [0UQ]{«fc}<* (4.93) 

Thus, the following substitution can be made. 

[ n2 ftL{«ft} s = [»2ft].,{afc} rf (4.94) 

2') The product of the second-order bending stiffness cross-term [n2ft] srf and dynamic bending 
nodal degrees of freedom {«&},/ has the form 

["2ft]. S( /{ rt ft} ( / = I [O] 7 [#]r[-4][0],/[( 'ol^Siob}^ (4.95) 

.s' 

but recall from numbers 1) and 2a) that 


[•4][0] rf [G]{«ft} ( / = {M»} d and = [ A H/KV]{«ftL (4.96) 
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Thus, the expression can be rewritten as follows. 




(4.97) 

s 

3') Similarly, the product of the second-order bending stiffness cross-term and static 

bending nodal degrees of freedom {///,}, has the form 

(4.98) 


>2»], ,.{«»),= / 


With the following substitutions 

[,4][((] s [r g ]{r/ 6 } s = {A*}, and rf{V#}.s = [Y'flUQ] M,/ 

the expression can be rewritten as 

["-/>], t.A a b}* = / 


(4.99) 


(4.100) 

S 

= [»-?»]., {«&},/ 

4') Finally, the product of the second-order bending stiffness cross-term [n2/,] (/s and dynamic 
bending nodal degrees of freedom {///,},/ can be written as 

(4.101) 




(4.104) 


but from number 1') it is known that 

[<9] I {(Tr'} y = =4- [^].s-[^^] {"*}(/ = [^](/[^]{° fc }. s (4.102) 

Thus, this term can be rewritten in the following form. 

W'hiis {«<>},/ = [" 2 &]</{ f, fc}.s (4.103) 

So, on the system level, implementation of the manipulations designated by 1') through 4') gives 

mU{A} s = 

munh = rnA^h 

[A r 2] ( , s {-4}. s = [N2)Mh 

[*2U^l} rf = [Y%{-4}, 

Expanding the system equations of motion, equations (4.68), results in 

[A/]{,4} ( / + [K]{A} S + [K]{A} d - [A'at]{4},. - [K AT ]{Ah + 

+ i[.Vl],{-4}, + i[A 1],{4}. + j[VI|,,).4}„ + i[ J Y2],M), + ^1, 

+ i[A’2]w{.4}, + + 3|A’2U(.4). + + J 1*2)4 'O. 

+ kv2],,).4},, = {/’</)} + (/■•(()} + {/’a? ) - 

5 (4.105) 
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and making the substitutions from equations (4.91) and (4.104) gives the following expression. 

+ [/'•]{, 1}, + [/>](. I),, - - [A-ArH-tb + i[.vi]J U. 

+ j[A'l]J.4},, + i[Al],,(,l),,+ i[.V2L{.b, + j[A’2L(.4) rf 

+ i(-V2].{,l},, + i[.Y2], ,{,!)» + {'!},, + I l|. + t), 

+ = {P(/>} + {/■■(/)! + )/'aji 

(4.106) 

Separating the terms involving the static displacements { A } s from those involving the dynamic 
displacements! .1 } (/ results in the following static and dynamic governing equations. 

([A] - [A'a/ 1 + j|.VI], + i[.V2], + [.V2], ,){,!), = { I’at ) 

W I i|. + ([A] - [ A" ax] + I-VIL+ j[.Vl], ( + [A’2], (4.107) 

+ j[-V2U{.4b = {/>(/)) + {F(/)} 

The matrix [_Y’2] rf , appearing in both equations, can be neglected because it is second order with 
respect to the small dynamic displacements The matrix [.V I ] (/ , appearing in the dynamic 

equation, can be neglected also because the product [.Vl] (/ {.4} (/ is second-order in {A} d as well. 
Therefore, the final static and dynamic governing equations are as follows. 

(|A) - [A'a,] + j[.Vl], + ::Y2].){.1|. = \l \ T I Static 

</ { »(..• + ([A] - l A' a/1 + [VI], + [.Y21.H.4), (4 - 108 > 

= {P(/)} + {FlO} - Dynamic 

In the applications being considered for this technology, there is often a combination of a 
flexible structures with a rigid boundary and combined thermal and mechanical loads. Thus, 
there is an opportunity for thermal buckling and an analysis capable of predicting such an event 
is necessary. Often it is also desirable to predict the post-buckled deflection of a structure 
subsequent to such an event. It will be shown later in this document that certain combinations 
of constituent materials in a SMAHC structure can exhibit a thermal buckling event and return 
the post-buckling deflection to zero. In order to capture and better understand this behavior, 
thermal buckling and post-buckling analyses are necessary. 

Thus, three types of analyses are required to study the response of the structures considered 
in this study: (1) thermal buckling analysis, (2) thermal post-buckling analysis, and (3) dynamic 
analysis. Each of these analyses are discussed in the following two chapters. Acoustic radiation 
analysis is presented in Chapter 7 to study the structural-acoustic response and noise transmission 
characteristics of panel-type structures. A flowchart of the overall thermomechanical analysis 
procedure is given in Figure 4.2. 
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• Perform thermal buckling analysis 

• Store bending eigenvector and/or 
membrane deflections 



• Perform thermal post-buckling analysis 

• Store bending and membrane deflections 



Perform vibration and 
acoustic radiation analysis 


Figure 4.2 Overall thermomechanical analysis flowchart. 
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Chapter 5 
Static Response 


The static response considered in this study is assumed to be solely due to a thermal load. 
Two analysis types are potentially needed to determine this response: (1) thermal buckling 
analysis and (2) thermal post-buckling analysis. Each analysis will be discussed in this chapter. 
Thermal buckling and post-buckling considerations necessitate the investigation of static stability. 
Therefore, the fundamental equations governing this behavior will be derived first. 

5.1 Static Stability and Nonlinear Analysis 

The idea of static stability commonly arises in the context of a structure subjected to 
compressive mechanical loads in an axial direction. The goal is typically to determine the 
critical load for which the structure loses stability and may be deflected significantly by a 
lateral load, from which it will not return upon lateral load removal (commonly referred to as 
buckling). In practice, a structure in such a condition often deflects (buckles) as a result of 
imperfections, without the application of lateral loads. Whether buckled by mechanical loads 
and/or imperfections, it is often desirable to calculate the deflected shape of a structure subsequent 
to buckling (often termed post-buckling deflected shape). Such an analysis usually involves 
geometrically nonlinear analysis as the deflections are typically large. 

In the case of thermoelastic response analysis, such as that considered here, a compressive 
load is generated by thermal expansion of a structure in the presence of mechanical boundary 
constraints. Thus, the goal is to determine the critical temperature at which the structure becomes 
unstable (also called the buckling temperature). Analysis of the thermal post-buckling deflection 
will also be of particular value. 

Stability of an equilibrium configuration can be investigated through the use of the adjacent 
equilibrium method [ 46 , 47 ], Application of this method entails imposing small increments to 
the displacements of a system in static equilibrium and examining the resulting adjacent con- 
figuration. Let {^(.4)} denote an equilibrium configuration for a structure with a displacement 
field given by {.4}. Now, let the static displacements be given an increment 

{A} — ► {,4} -I- {6,4} such that {6,4 J } ~ {0} (5.1) 

where the {6 .4} are infinitesimal displacements. The adjacent static configuration \ ^{A + 6.4)} 
can be approximated by a truncated Taylor Series expansion: 

{*(.4 + M)) = {*(.4)} + +0({«,4-'}) (5.2) 
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One could proceed to determine the derivative of the static configuration with respect to the 
static displacements, but the necessary ingredients for an equivalent procedure has already been 
developed in the previous chapter and can be employed here. 

Recall that the equilibrium configuration { ( 4 ) } is governed by the static system in 
equations (4.108): 


W'D) = ([/'] - [A'at] + j[A'l] + }[.Y2]){,1) - {/’ A , } = {(I) (5.3) 

where the static subscript .* has been omitted for simplicity. The expression for the increment in 
the static displacements, equation (5.1), is analogous to the separation of the static and dynamic 
components of the total displacements given in equation (4.58). Substitution of {1} + { <*> .1 } into 
equation (5.3) gives the adjacent static configuration that is not necessarily in equilibrium: 

{*(.1 + Mil = ([A] - [Aat] + i[ Vl(.l + M)| + ^Y2(4 + 64)])({4} + {64})- {Pat} 

(5.4) 

The manipulation of terms in the discussion following equation (4.58) allows the first- and 
second-order nonlinear stiffness matrices for the static system in equation (5.4) to be written as 

[A 1] = [A 1]„ 4- [A 1 h 

[A -] = [ V2] (l + [X2] M + [A'2] w , + [A-2], 

where the subscripts 0 and 6 indicate the influence of the displacements { A } and { A .4 } , respec- 
tively. 

Expanding [ VI] and [ Y2] in equation (5.4) gives the following expression for the adjacent 
configuration. 


{^(4 + 6 A ) } — [[A ] — [/v at] + -([-Y i ]o + [A 1 ]<s ) 


+ :j([A 2] 0 + [A 2 ] oa + [A2]^ (l + [A 2]^ )] ( {4} + {64}) - {/{at} 


(5.6) 


Further manipulation allows several terms in equation (5.6) to be written in forms analogous to 
equations (4.91) and (4.104): 

[A’l]*{4} = [A 1]„{64} 

1*2],* {4} = [A r 2]„{64} 

[A -] 0/ s {64} = [A'2]/, { 1} (5.7) 

[A2]„ 0 {4} = [A 2] 0 {64} 

[A r 2]^ ffl {64} = [.Y2],{4} 

Expanding equation (5.6) gives 

{*(.-1 + M)} = [A]{.4) + [A]{A4) - [AatK-4) - [A'at]{A1) + {[.VI], ,(.4) 

+ + }[.VI],{.4) + i[.Vl| t {«.4) + i[.V2]„{,l) 

+ + i(A-2] M {-4) + }[-V2]„»{A4) + j[.V21, 0 {4} 

+ 4v 4„{(./1) + k\’2|„{.4} + }(-V2],{A4) - {/' A7 | 

(5.8) 
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(5.9) 


and making the substitutions from equations (5.7) results in 

{ty(.4 + (V 4)} = [A'] { .4 } + [A']{<S.4} - [A'^r]{.4 } - [A avH^'U + ^[A 

+ ^Vl]„{*4} + i[- vl ]ii{ <s -'} + r,l YI U <s1 } + x -]<> I M 

+ ^[-V*2] 0 { A.4} + -],){(' -4} + -[A 2]^ {.4} + —[A 2] 0 {A,4} 

+ ^[A r 2]^ { .4 } + i[.V2]^{.4} + ^A"2],{M} - {/ j a/ } = {0} 

Separating terms in {.4} and { A .4 } , neglecting second- and higher-order terms in {0.4), and 
omitting the subscript 0 gives 

{'I' (.4 + 0.4)} = ([A ] - [A at] + -[A 1] + -[A 2] ) { .4 } — {I at} (5.10) 

+ ([A] — [A A7 ] -I- [A 1] + [A 2]) {0.4} 

In order for this adjacent static configuration to be in equilibrium we must have {^(.4 4- 0.4)} = 
{0}. Then, the form for the incremental static system of equations becomes evident when the 
relations in equations (5.2) and (5.3) are compared to that in equation (5.10). 

([A] - [A a7 ] + [A' 1 ] + [A 2]) {0.4} = {0} (5.1 D 

as the original equilibrium configuration {^(.4)} is identically satisfied. It can also be seen 
that the first-order derivative of the original equilibrium configuration with respect to the static 
displacement vector is given by 

= [A-,„] = [A] - |A a T) + [VI] + [,V2] (5.12) 

d{A\ 

which is called the tangent stiffness matrix. 

It will be shown in the following sections that Equations (5.10)— (5.12) are useful for thermal 
buckling and post-buckling analyses. The incremental static system of equations will be useful for 
investigation of thermal buckling as the expression must be satisfied for an adjacent equilibrium 
to exist. The tangent stiffness matrix will be useful in calculation of static deflections due to 
thermal moments/bending-membrane coupling or post-buckling deflections, where a Newton- 
Raphson approach will require the derivative of the stiffness coefficients. These procedures will 
be described subsequently. 

5.2 Thermal Buckling Analysis 

A thermal buckling analysis is carried out to determine the critical (buckling) temperature 
distribution of a structure r rr (.r,ty) = A T rT (.v.y) + 7 0 , where T 0 is the uniform reference 
or ambient temperature distribution. This type of analysis is useful for determination of a 
need for thermal post-buckling analysis for a particular configuration and loading condition. 
Although buckling analyses are inherently linear, and the system in concern is decidedly 
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nonlinear (temperature-dependent properties), there is at least one way to seek a solution within 
the framework of the present formulation. Furthermore, a thermal buckling solution is very 
desirable because of the insight it affords, the computational efficiency for cases when thermal 
post-buckling analyses are not warranted, and the ability to provide an initial trial solution for 
thermal post-buckling analyses when they are necessary. All of these points will be made more 
evident in the following discussions. 

Equations governing thermal buckling can be obtained from the static equations of equi- 
librium, equations (4.108), and the incremental static system of equations, equation (5.11). As 
stated previously, thermal buckling is a stability problem, indicating a condition for the on- 
set of buckling. Therefore, the nonlinear stiffness terms due to out-of-plane deflection can be 
neglected. The stiffness terms involving bending-membrane coupling can also be neglected be- 
cause any increase in temperature will result in out-of-plane deflection. An exception to this 
exists for antisymmetric angle-ply laminates which is a special case that will not be treated 
here. Moment resultant terms are also neglected for similar reasons, i.e., the thermal load is 
assumed constant through the thickness of the structure. Cases for which there is either bending- 
membrane coupling or thermal moment effects should be treated with a thermal post-buckling 
procedure, described in the next section, to determine the buckling temperature and static thermal 
deflections. The static governing equations from (4.108), with the aforementioned restrictions, 
simplify to the following. 

){£} = UU 

The term [A'4] is retained because it is dependent upon the membrane displacements {.4,„} only. 


'Kb 

0 


A'av 

o' 

1 

+ 2 

A 1 b 

O' 

0 

Am 


0 

0 

0 

0 


5.2.1 Constant Matrix Material Properties 


The proposed solution to the thermal buckling problem, discussed subsequently, entails 
separation of the SMA and matrix material thermal effects. Recall the SMA hybrid composite 
lamina constitutive relations in principal material coordinates, after separating the SMA and 
matrix thermal effects, given in equation (2.29). It is noted that the material system constituents 
are not fully separated in that equation because of their combined effects in the reduced stiffness 
matrix. However, the coefficients C|„ and n w essentially render them separate in the 1-direction 
and will lead to a useful form for the governing equations. Incorporation of those constitutive 
relations into the finite element formulation results in a system of static governing equations 
analogous to those in equations (4.108). The aforementioned assumptions inherent to thermal 
buckling analysis reduce this static system to the following. 


(\Kb 0 

VL 


A AT' 0 
0 0 


A A T 
0 


+ ■ 


A 1 b 
0 


0 

0 




0 } 

nAl J m 

(5.14) 


The subscripts a and in indicate SMA and matrix material, respectively. It is obvious that the 
bending and membrane equations of the static problem are now de-coupled. The membrane 
displacements due to an applied thermal load that varies only in the plane of the structure can 
thus be determined from the equation 


[A »( ] { A m } — { Pfn AT },, 4“ { f niAT } ik 


(5.15) 
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For a given applied temperature distribution AY f.r. //), equation (4.55) can be used to form 
the thermal force vectors {P w at} u and {P m AT and equation (5.15) can be solved for the 
membrane displacements: 

{ '!„,} = [A m ] ([^nAt},, + { ^ ih AT } ) (5 [A] 


The stability of this thermoelastic equilibrium configuration can be investigated through the use 
of the incremental static equations (5.11). 

The incremental static system of equations can be simplified as above by neglecting nonlinear 
stiffness terms due to out-of-plane displacements, bending-membrane coupling stiffnesses, and 
moment resultants: 


(\K h 0 

^ [ 0 Am 


A'a'/ 0 

0 0 


A'a 7 
0 


0 

0 



VI b 
0 



+ 


-VI ft 
0 




(5.17) 


Note that the first-order nonlinear stiffness matrix due to membrane displacements has been 
separated into components evaluated with membrane displacements caused by SMA and matrix 
material thermal effects, respectively. Thus, the equations governing the stability of the out-of- 
plane displacements can be written as 


([Kb] - [A A7’]« - [A'a7 L + [- V ><-]«. + = {0} 


(5.18) 


Consider a temperature distribution change of the following form 

AT(x.y) = XATi(x.y) (5.19) 


where A is a scalar. For a matrix material with constant material properties, the thermal force 
vector and geometric stiffness matrix due to the thermal effects are linear functions of the 
temperature change: 

{PmAT} mAT{r.y) ^{P»‘A1 } mATi(J',y) (5 70 ) 

[A'a:tLa7’(x,<,) = "MA A7’]r,i AT) (x.;y) 

A similar relation can be written for the corresponding first-order nonlinear stiffness matrix 

[A H] ( /,A7(x.vl = '^[4 l/)] m .A7’i(/-,;i/) (5—1) 

because equations (5.15) are linear and [Nl b ] m is linearly dependent upon the membrane 
displacements due to matrix material thermal effects. 

Analogous expressions cannot be written for the thermal force vector {Pm at) aAT( 
geometric stiffness matrix [A at] ( ,a7'(x,;v)’ anc ^ firsVorder nonlinear stiffness matrix [A 1 t>] a AT\ r ,u} 
associated with the SMA thermal effects over the same temperature change A T(x.y) = 
AA7’i (x.y) as they are all nonlinear functions of temperature. The stiffness matrix [A fc ] is 
also nonlinear with temperature because of the influence of the SMA constituent. At the present, 
let all of the terms involving the SMA thermal effects be simply evaluated at the imposed 
temperature distribution and taken as constant; the rationale for this approach will be explained 
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subsequently. In this case, equation (5.18) can be rewritten for the temperature distribution 
change AT(.r.y) = Wl\(x.y) as 

([A?,] - [Aa7']„ ~ -M A A7 1 ],,, ± [A 4] ( , + A[A 4],„){A4/,} = { 0 } (5.22) 

which can be expressed as an eigenvalue problem 

([A*] - [A at],, + [-Y4DM = A( [A’av ]„, - [-Vl fc ] J{o} (5.23) 

A critical temperature change exists if there is a solution to equation (5.23) for which { c>} is 
nontrivial. The first eigen-solution gives the critical temperature change from the expression 
AT CI (.r, //) = Ai A7’i(.r. //) and the corresponding buckling mode shape {o}j. 

An iterative procedure can be employed to account for the nonlinear functionality of the 
terms involving the SMA thermal effects in equation (5.23). The iterative procedure entails 
imposing any desired temperature distribution 7) (.<•.//) = 7b + A7i(.r.;/) to evaluate the 
material properties. Then, equations (2.37) and (4.55) an be used to evaluate the SMA and 
matrix material thermal force vectors { P,„ a 7 }„ and { A„,at}„, and the corresponding geometric 
stiffness matrices [A’a 7’] (I and [A' a E quations (5.15) can be solved for the membrane 
displacements and {.4„, } m , which can then be used in equations (4.18) and (4.55) to 

form the first-order nonlinear stiffness matrices [.Yl/,] (( and [ATj,] m . Subsequently, the bending 
eigenvalue problem, equation (5.23), can be formed and solved for an estimate of the buckling 
temperature /_>(./■. y) = A{7i(.r. //) = Aj[7b ± AT) (.;•,/;)] and mode shape { o } | , where the 
superscript indicates iteration number. This process is repeated until convergence is achieved 
(A" = 1.000 ±.001), at which point usage of the linearized eigenvalue problem has been 
rendered appropriate and accurate because A7’„(.i\ y) is such a small quantity that the material 
properties are essentially constant. 

5.2.2 Non-Constant Matrix Material Properties 

The assumption that the matrix material properties were constant with temperature was 
implicit to the equations and solution methodology developed in the previous section. Under 
such conditions, a rationale was developed to separate the linear and nonlinear components of the 
governing equations and use an iterative approach to compensate for the nonlinearity. Recall that 
during the linear/nonlinear separation process, a comment was made that the reduced stiffness 
matrix disallowed complete separation of the constituent material effects. This point was also 
made in Chapter 2. Thus, some nonlinear effects (e.g., modulus of the SMA) do appear on 
the right-hand side of equation (5.23) and the solution method of the previous section must be 
applicable for at least some forms of temperature dependence in the matrix material. 

This has, in fact, been observed for all of the cases considered in this study, which include 
cases for which all of the matrix material properties vary significantly with temperature. The 
validity of the thermal buckling solution for such cases was established by comparing with 
a thermal post-buckling method. The thermal post-buckling method serves as another means 
of determining a buckling temperature, as the temperature at which the thermal post-buckling 
deflection departs from zero. In the event that a thermal buckling solution could not be achieved, 
an initial guess would be formulated from the thermal buckling problem and used to start a 
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thermal post-buckling solution. More detail on this approach will be given following the post- 
buckling solution method description at the end of the next section. A flow chart exemplifying 
the thermal buckling solution procedure is shown in Figure 5.1. 



Figure 5.1 Thermal buckling calculation flowchart. (Continued) . . . 
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1 


If AT(x,y) = AT p (x,y), 
store membrane deflections 


Calculate absolute 
norm error, E 


AT(x,y) = AT cr (x,y) 



Dynamic analysis 



Increase maxiter or 
choose another AT p (x,y) 


Store membrane deflections {A m } C r i 
and bending Eigenvector {<[>} 


Post-buckling analysis 


Figure 5.1 Thermal buckling calculation flowchart. 


It will be shown in Chapter 8 that two buckling temperatures can exist for some SMA hybrid 
composite structures. Details will be deferred to that discussion, but some general comments 
are useful here. This effect is a manifestation of a material system with a SMA constituent that 
has an austenite start temperature .4 S greater than ambient temperature and is not to be confused 
with higher-order buckling modes occurring at temperatures above initial thermal buckling. A 
typical sequence of events in such a case can be described as follows. At temperatures less than 
.4 S , the SMA is not activated and behaves like a conventional linear thermoelastic material. So 
the SMAHC structure will buckle if a critical temperature exists below .4,. The SMA becomes 
activated at .4., and the recovery stress rapidly increases with temperature, which overcomes the 
compressive thermal stress and reduces the thermal deflection (potentially eliminating it). At 
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higher temperatures, the rate of increase of the recovery stress with temperature diminishes and 
may be surpassed by the rate of increase in matrix compressive thermal stress. Further increases 
in temperature will cause net-compressive stress development, causing the thermal deflection to 
grow again or lead to a second incidence of thermal buckling for a structure that was rendered flat. 

5.3 Thermal Post-buckling Analysis 

Once the buckling temperature is determined, the next step will depend on the temperature 
for the particular application of concern. For temperatures greater than the critical temperature 
(T(.r.y) > //)), a post-buckling analysis must be performed to determine the thermal 

deflection. The thermal post-buckling solution is started at the buckling temperature T rr {.t\ y) 
for expediency and ease of solution. Starting at the ambient temperature can be much more 
computationally intensive (for cases with large A T cr (.v,y)) and is numerically ill-conditioned. 
Further remarks about initiating the post-buckling solution will be given later in this section. For 
the cases described above, where two buckling temperatures can exist, it is advisable (although 
not necessary) to initiate a thermal post-buckling analysis at the first buckling temperature to 
assess the evolution of the static deflection in arriving at the desired operating temperature. The 
governing equations for thermal post-buckling deflection are given in equations (4.108) 

A/i] - [/IAV'1 + j[.VI] + A -’ ) ( 'll = {''Al l (5-24) 

where the static subscript * has been omitted. Note that there is no need to separate the SMA 
and matrix material thermal effects in thermal post-buckling analysis. One effective approach 
for solving this nonlinear system of equations is the Newton-Raphson iterative method. 

Consider a static system described by equation (5.24) at the buckling temperature distribution 
T cr (.r.y) and denote this equilibrium configuration by {'J'(zl)}. A small increment in the 
temperature distribution T(x , y) = T cr (.i\ y)-\-bT{.i, y) will result in the formation of an adjacent 
equilibrium configuration denoted by {$(.4 + 74)}. This adjacent equilibrium configuration can 
be derived from a truncated Taylor series expansion as shown in equation (5.10): 

{ 4* ( 4 4- 74)} = ( [A ] — [A ^yf] + -[A 1 ] + -[A 2]) {4} — {Pat } (5.25) 

+ {[!<} - [Pat] + [A r l] + [AA>]){74} = {0} 

The first collection of terms on the right-hand side of this equation is no longer a statement 
of static equilibrium 

{^(4)} ± ([A'] - [Pat] + ^[A r l] + ^[.V2]){4} - {Pat} + {0} (5-26) 

as the stiffnesses and thermal load vector are evaluated at T(.i\y ) = T cr (x.y) + 77’(.r, y), while 
{4} is associated with the system at temperature 7’ rr (.r. //). So, the expression for the adjacent 
equilibrium configuration can be rewritten in terms of the tangent stiffness matrix, equation 
(5.12), and an incremental load vector 

[A,«„]{74} = {*Pat} < 5 - 27 > 
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where the incremental load vector is defined by 

{W'a7 ' } = - ([A] - [Aa,-] + i[:VI] + j[A : 2]){7|) (5.28) 

Although a small increment in temperature may render the static system essentially linear with 
temperature, the new equilibrium configuration is also governed by the geometric nonlinearity 
of the first- and second-order nonlinear stiffness matrices, which are linearly and quadratically 
dependent upon the out-of-plane displacements. Thus, an iterative procedure is necessary to 
converge on the new equilibrium configuration. 

The incremental form of equation (5.24) for the j-th iteration can be written as 

[/Han ]j{6A}j = {SPat}j (5.29) 

The tangent stiffness and incremental load vector for the subsequent iteration 
[/Han] ;+ i and Way }, + j are updated by using the updated thermal deflections 

{,1} = {.4} + {M} ; (5.30) 


The iteration procedure is continued, with the previous iteration’s deflections as an estimate, 
until the deflection estimate approaches the true equilibrium configuration. As the deflection 
estimate approaches the equilibrium configuration, the expression for the incremental load vector 
approaches an expression of equilibrium, at which point the incremental load vector {bl^r } ) 
and, consequently, the incremental displacements {£.4} • approach zero. Convergence is achieved 
when either incremental values are reduced to a specified small quantity defined by the error 
criterion. A maximum norm error criterion was used in this study, as described by Bergan and 
Clough [48], 

As mentioned previously, the thermal post-buckling analysis is sensitive to the starting 
temperature and the corresponding trial solution. An initial trial solution for the thermal post- 
buckling procedure, starting at the buckling temperature, can be constructed from the thermal 
buckling solution. The membrane deflections, determined from equations (5.16) evaluated at 
the (first) buckling temperature, are used directly. The buckling mode shape, determined from 
equation (5.23), is normalized to its maximum displacement entry and to a fraction of the 
structure thickness (half-thickness was found to perform well) to give the out-of-plane portion 
of the initial trial solution a physically scaled value: 


f A b \ f0.5//{o}/o l( . mu , 1 

\ -4 m J \ { A '« } rr J 


(5.31) 


Convergence on the displacements due to a rather large temperature increment AT(x.y) 
can typically be improved by splitting the large temperature increment into a number of smaller 
increments: 


ATi(x.y) 


AV(-r.//) 

• V A7 


(5.32) 


Iteration is performed at each of the smaller temperature increments in succession until the 
converged equilibrium configuration (post-buckled static solution) is achieved for the desired 
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temperature. The converged displacement vector at each of the smaller temperature increments 
is scaled and used as a trial vector for the next temperature increment : 


{Ah = {A},_ : 


Tj(.v,y)-T 0 (j'.u) 
Ti-\{.v, y) - T 0 {.r. y) 


(5.33) 


Note that a mechanism for re-initialization of the out-of-plane displacements is necessary because 
of the potential for elimination of the post-buckling deflection due to SMA recovery forces for 
the cases with two buckling temperatures. A flow chart demonstrating the thermal post-buckling 
solution procedure is shown in Figure 5.2. 

One final remark concerns any potential cases for which convergence on a thermal buckling 
solution is not possible. In such cases, a thermal post-buckling solution should be initiated 
from the reference temperature and should proceed at temperature increments until a thermal 
post-buckling deflection is developed. At each increment in temperature, the in-plane problem 
of the buckling solution should be solved for the in-plane deflections. The in-plane deflections 
will be used to form the thermal buckling eigenvalue problem. The eigenvector resulting from a 
single eigen-solution would then be used, along with the in-plane deflections, to form an initial 
guess for the post-buckling solution at that temperature. Once a thermal post-buckling deflection 
is achieved, a refined estimate of the buckling temperature can be determined by restarting the 
post-buckling solution in the vicinity of the first solution with a smaller temperature increment. 
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Figure 5.2 Thermal post-buckling calculation flowchart. (Continued) 
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5.4 Analytical Validation 

In the previous sections, a formulation and approach for prediction of the static stability 
(thermal buckling) and large-deflection (thermal post-buckling) of SMA hybrid composite panel- 
type structures was derived. The formulation and implementation will now be validated through 
comparison with fundamental cases and experimental measurements. Comparisons will be made 
with known solutions for the cases of a simply-supported isotropic beam, simply-supported and 
clamped isotropic plates, and a simply-supported symmetrically-laminated composite plate to 
benchmark the analysis. Comparisons will be made with experimental measurements for the 
case of a clamped-clamped SMAHC beam to fully validate the thermomechanical model for 
static response. 

Note that the effects of initial stresses, initial deflections, thermal moments due to a 
temperature gradient through the thickness, and bending-membrane coupling (non-symmetric 
composite laminates) are not included in the analysis, although solution with such conditions 
is quite tractable by the thermal post-buckling solution technique. Therefore, the thermal 
buckling cases presented are true bifurcation problems and, as such, the resulting post-buckling 
deflection could be positive or negative. In practice, the predicted buckling direction is typically 
predetermined by the sign of the initial guess. The sense of the buckling deflection in a particular 
experimental configuration is usually influenced by some combination of imperfections including 
initial stresses, initial deflections, and anomalies in boundary conditions. The positive sense will 
be taken in the purely numerical results to be presented and the sense will be chosen to match 
the experimental behavior when comparing with experimental measurement. 

5.4.1 Simply-Supported Isotropic Beam 

Consider a simply-supported aluminum beam with dimensions of 0.5588x0.0254x0.0023 m 
(18x1x0.090 inches) and subject to a uniform temperature increase. The material properties for 
the specimen are shown in Table 5.1. 


Table 5.1 

Material 

Young’s Modulus 

68.95 

Shear Modulus 

26.55 

Poisson’s Ratio 

0.3 

Coeff. of Thermal Exp. 

22.3 /. 


properties for aluminum specimens. 


GPa 

10 Msi 

GPa 

3.85 Msi 

<c/°F 

12.4 // f /°F 


A classical solution for the critical temperature change to induce buckling can be easily 
derived based upon beam-column theory, as shown by Boley and Weiner [49] for example. For 
a uniform temperature distribution, the axial load is given by 


P T = AEoAT 


(5.34) 
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where A is the beam cross-sectional area, E is the Young's modulus, and o is the thermal 
expansion coefficient. Thermal buckling occurs when the axial load reaches the Euler buckling 
load, defined by the familiar expression 


where A is the area moment of inertia, and /, is the beam length. Thus, the critical temperature 
change is found to be 


A 7 = 


77-1 


(5.36) 


This gives a numerical value of l.658°F for the predicted critical temperature change from 
beam-column analysis as compared to 1.643°F for the present finite element analysis with a 
36x2-element mesh. The slight discrepancy is attributable to the fact that the classical analysis 
does not include transverse effects because it is a one-dimensional analysis. 


5.4.2 Simply-Supported and Clamped Isotropic Plates 

A classical solution for the critical temperature change to induce buckling for an simply- 
supported isotropic plate is also given by Boley and Weiner [49]. For a uniform temperature 
distribution, the in-plane thermal force is given by 

h_ 

2 

N r = nE f NT (h = <\EhXI' (5-37) 


where h is the thickness of the plate. Thermal buckling occurs when the in-plane thermal load 
reaches the reaches a critical value, given by 


N T cr = (1 


< 7 “ 

l + lA 


'[) 


a- 


(5.38) 


where v is the Poisson’s ratio, a and b are the panel length and width, and D is the flexural 
rigidity given by 


D = 


Eh :i 

12(1 


(5.39) 


Thus, the critical temperature change is given by the following expression 


A7V 


— v 



--I) 
a 2 a Eh 


(5.40) 


As an example, consider a simply-supported aluminum plate with dimensions of 
0.3556x0.254x0.001 m (14x10x0.040 inches) and the material properties given in Table 
5.1. The classical solution for this case gives a numerical value of 1.233°F for the predicted 
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Figure 5.3 Normalized maximum deflection versus dimensionless thermal load for a clamped, 
0.102x0.102x0.001 m (4x4x0.040 inch), aluminum plate with uniform thermal load. 


critical temperature change, which is matched exactly by the present finite element analysis 
with a 14 x 10 element mesh. 

A classical solution for the thermal post-buckling deflection of a clamped isotropic plate 
was developed by Paul [50]. He used a twenty-five term series approximation in a classical 
formulation with von Karman deformation assumptions. One case he analyzed consisted of 
a clamped aluminum plate with dimensions of 0.102x0.102x0.001 m (4x4x0.040 inches) with 
material properties as given in Table 5.1, except with a slightly different thermal expansion 
coefficient o=22.5 //e/°C (12.5 //e/°F) and a uniform thermal load. A plot of the normalized 
center deflection versus dimensionless thermal load for this case is shown in Figure 5.3. It 
can be seen that the present finite element analysis is in excellent agreement with the classical 
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solution. The fine resolution Sat = 100 for the finite element results was used for presentation 
purposes only. Convergence on the solution at the final temperature was readily achieved 
with only one temperature step. In fact, the solution was much more expedient albeit not 
conducive to conveying information about the deflection behavior at intermediate temperatures. 
The dimensionless thermal load shown as the abscissa of the plot is given by Paul as 


oA77r 

Load = 12(1 -I- v )- — 77 - 7 - 

7 ! — // “ 


(5.41) 


The range on the load from approximately 5.3 to 12 is equivalent to a A 7 from zero to 15.6 C 
(60° F). 


5.4.3 Simply-Supported Composite Plate 

Meyers and Hyer [5 1 ] employed a Rayleigh-Ritz method to developed a classical solution 
for the buckling temperature and the thermal post-buckling deflection of simply-supported, 
symmetrically-laminated composite plates. They analyzed plates with fixed and sliding simple 
supports. For example, Meyers and Hyer reported a buckling temperature of 20.8°C (69.4° F) 
for a graphite-epoxy plate with dimensions of 0.152x0.152x0.001 m (6x6x0.040 inches), quasi- 
isotropic lamination of (±45/0/90) s , fixed simple supports, a uniform temperature distribution, 
and material properties as shown in Table 5.2. The corresponding solution from the present 
finite element analysis is 20.6°C (69.05°F) with a 6 x 6 -element mesh. 


Table 5.2 Material properties in principal material coordinates for a graphite-epoxy lamina. 


Young’s Modulus, Ej 

155.1 GPa 

22.5 Msi 

Young’s Modulus, E 2 

8.07 GPa 

1.17 Msi 

Shear Modulus, G 12 

4.55 GPa 

0.66 Msi 

Poisson’s Ratio, v \2 

0.22 


Coeff. of Thermal Exp., 0 1 

-0.07 fit PC 

-0.04 fttPF 

Coeff. of Thermal Exp., 02 

30.1 ftcP C 

16.7 fitP F 


Corresponding thermal post-buckling deflections predicted by the classical and finite element 
analyses are shown in Figure 5.4. In this case, normalized center deflection is plotted against 
normalized temperature change. Excellent agreement is achieved between the classical and finite 
element analyses. Again, a fine resolution ft at = 100 was used in the finite element analysis 
for presentation purposes, but more expedient calculation of the final thermal post-buckling 
deflection was achieved with only one temperature step. 

Recall that both cases for which thermal post-buckling deflections were calculated are 
classified as bifurcation problems. The analogous negative deflection results could be attained by 
using a thermal buckling eigenvector of opposite sign in the formation of an initial trial vector. 
Note that this ambiguity is of no consequence to the thermal buckling solution as it is the relative 
magnitudes in the buckling eigenvector that are important. 
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Figure 5.4 Normalized maximum deflection versus normalized AT for a simply-supported, 
0.152x0.152x0.001 m (6x6x0.040 inches), (±45/0/90) s graphite-epoxy laminate. 

5.5 Experimental Validation 

The thermal buckling and post-buckling solutions were validated against known solutions 
for conventional structures in the previous section. Attention will now turn to experimental 
measurement of the thermoelastic response of a SMAHC beam. The experimental results will 
be used to demonstrate some of the phenomena discussed earlier in this chapter and to validate 
the formulation and solution methods for the case of a SMAHC structure. 

5.5.1 Clamped SMAHC Beam 

One of the SMAHC beam specimens is shown mounted in an aluminum fixture designed to 
provide clamped boundary conditions for beam specimens with this geometry in Figure 5.5. The 
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beam has an unsupported length of 0.4572 m ( 18 in) within the mechanical grips. The fixture was 
also designed to provide electrical connection and mechanical restraint for the Nitinol ribbon leads 
protruding from the ends of the SMAHC beam specimens. The mechanical and electrical grips 
are indicated in the figure by the symbols “m” and “e”, respectively. The electrical connections 
are isolated from the fixture by a layer of fibrous ceramic insulation. The beam is also thermally 
isolated from the fixture by a layer of this insulation on both sides. The beam is heated by DC 
electrical current controlled by a thermocouple measurement located approximately on the beam 
centerline and 0.0127 m (0.5 in) from the right mechanical grip, indicated in the figure by “tc”. 
The thermal controller consists of a process controller (Omega model CN77353-PV), and a DC 
power supply (HP model 6652A). Thermocouples also monitor the temperature at each end of 
the fixture, indicated in Figure 5.5, and one measures the ambient air temperature. 



Figure 5.5 Front view of the SMAHC beam specimen mounted 
in clamping fixture used for static and dynamic testing. 


A rear view of the experimental configuration is shown in Figure 5.6, which shows the 
installation of a laser displacement transducer. The laser displacement transducer (Aromat 
model LM100) was attached to a three-axis optical positioning traverse, and the assembly was 
mounted in the vicinity of the beam mid-span. The beam mechanical grips and the ribbon 
mechanical/electrical grips were torqued to 13.56 N»m (120 in*lbf). The three-axis traverse 
was used to position the laser impingement position at the beam mid-span and to adjust the 
laser transducer/beam surface distance to produce zero offset. Changes in beam and fixture 
temperatures during the gripping process and prior to starting the thermal controller were 
insignificant. The beam was heated from ambient temperature to 121.1°C (250°F) over a duration 
of 5 minutes and 1 1 seconds, then allowed to cool naturally to ambient temperature again. The 
overall duration of the test was 13 minutes and 26 seconds. Measurements of the ambient air, 
right/left fixture end, and beam temperatures and the laser transducer output were recorded at 
approximately 2.5 second intervals. 

The beam was modeled with a 36x2 mesh of the 24 degree-of-freedom rectangular plate 
elements described in Chapter 4. This test constituted thermal cycle number four (4) on the 
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Figure 5.6 Experimental configuration for measurement of 
the static thermoelastic response of a SMAHC beam specimen. 


Nitinol ribbon in this particular beam specimen. Element properties were constructed from 
the appropriate measured material properties described in Chapter 3 and given in Tables 3.5 
and 3.8. Although the beam ends are held firmly within the mechanical grips, the boundary 
conditions are not well approximated by perfectly clamped conditions. It will be shown in the 
Experimental Validation section of Chapter 6 that it was necessary to alter the clamped boundary 
conditions with torsional springs, with a spring rate of 1.98 N*m/rad (17.5 in*lbf/rad), applied to 
the longitudinal-slope degrees of freedom to match the fundamental natural frequency. It will 
be assumed that the boundary conditions for the thermal post-buckling solution are essentially 
the same as that for the dynamic analysis. In previous experiments, the temperature distribution 
generated under these conditions was measured with an infrared thermal imaging camera. It was 
found that the temperature distribution was remarkably uniform along the unsupported length of 
the beam. Therefore, the thermal load was modeled as a uniform temperature distribution. 

A comparison between the measured and predicted thermal post-buckling deflection (nor- 
malized to the beam thickness) at the beam mid-span is shown in Figure 5.7. Note that the 
measured deflection does exhibit the first part of the complex thermoelasitic behavior described 
earlier in this chapter; buckling at a temperature below .4,, achieving a maximum deflection 
where the Nitinol recovery stress starts to dominate, and returns to a flat configuration for the 
remainder of the test. The reason that a second instance of buckling is not exhibited is that the 
volume fraction of Nitinol is high enough in this case to drive the predicted second buckling 
temperature well beyond the glass-transition temperature for the matrix material. 

The predicted post-buckling behavior is in excellent qualitative agreement with the experi- 
mental measurement. The main differences between the predicted and measured responses are 
the slope of the curve from initial buckling to the maximum deflection point, the magnitude of 
the maximum deflection, and the temperature at which the post-buckling deflection is returned to 
zero. These discrepancies are attributable to the fact that these systems are extremely sensitive 
to a multitude of factors including initial deflection, initial stress, boundary conditions, etc. For 
example, initial deflection in the physical beam is the reason for the more gradual growth of 
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Figure 5.7 Measured and predicted thermal post-buckling deflection 
(normalized to the beam thickness) at the mid-span of the SMAHC beam. 


thermal deflection in the experimental results. The discrepancies in peak deflection and deflection 
elimination temperature are attributable to boundary conditions. 

Note that the in-plane boundary conditions are modeled as immovable. Recall that in the 
physical system the beam continues through the mechanical grips and the Nitinol leads are 
terminated beyond that in the electrical grips. Only the ribbon/matrix interface supports the 
recovery force within the mechanical grip area. The main restraint for the Nitinol is the 
mechanical termination at the electrical grips. Thus, the compliance of the entire support 
structure, from one electrical grip to the other, needs to be considered. It was found that 
the layer of electrical insulation between the electrical and mechanical grips was the cause for 
the discrepancies. A compression test was performed on the insulating material and the modulus 
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was found to be less than 137.9 MPa (20 ksi). While the exact stiffness of the boundary would 
be difficult to quantify, there is no doubt that strain in the insulation is partially relieving the 
recovery stress and causing a diminished stiffening effect in the physical system. Changes 
in the model to reflect the compliant in-plane boundaries would move the predicted buckling 
temperature, peak deflection, and deflection elimination temperature to better agreement with 
experimental measurement. 

In this chapter, the static stability of a structure described by the governing equations 
developed in the previous chapter was analyzed and solution procedures were developed to 
predict the thermal buckling and post-buckling behavior of SMAHC structures. The formulation 
and procedures were validated against known solutions prior to attempting correlation with 
experiment. Experimental results from a SMAHC beam specimen confirm complex thermoelastic 
behavior that was predicted by the model. Excellent agreement was achieved between prediction 
and experimental measurement and reasons for existing discrepancies were discussed. A similar 
approach will be taken in the next chapter, where the dynamic response performance of a 
SMAHC beam will be demonstrated and compared with numerical prediction. 
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Chapter 6 
Dynamic Response 


The analytical techniques and computational methods used to solve the equations governing 
the dynamic response of SMA hybrid composite panel-type structures subjected to mechanical 
and thermal loads are treated in this chapter. The thermal load is assumed to be steady-state 
and uniform through the thickness of the structure, i.e., T = T(.r.y). Mechanical loads are 
assumed to consist of surface pressure, concentrated forces, and inertial (base acceleration) 
loads. Continuous and discrete frequency domain analyses are presented. Validation of the 
analyses with known solutions for fundamental cases and with experimental results for several 
cases including a SMA hybrid composite beam are presented. 

6.1 Equations of Motion 

The equations of motion for the system described above were derived in Chapter 4, equation 
(4.108), and can be rewritten in the compact form 

[M] { A} + [A] { .4 } = { P( t )} + { F(t ) } - [M] { V g } (6-1) 

where [A] is the stiffness matrix for the flat or buckled/deflected structure. For a flat (unbuckled) 
structure 

[A] = [A] - [A'at] + [A r 1] unbuckled (62) 

[A] = [A',„„] = [7\] - [Aa 7'1 + [A r l] + [.Y2] buckled 

Note that the only nonzero entries in [AT] for the unbuckled case are due to [A H] resulting 
from the in-plane displacements, see equation (5.16), and the stiffness matrix for the thermally 
buckled structure is precisely the converged tangent stiffness matrix from the thermal post- 
buckling solution. For a thermally buckled/de flee ted structure 

[A'] = ( [A ] - [A',v at ] + [A r 1 ] + [-V2] ) (6.3) 

where the nonlinear stiffness matrices [AT] and [N2] are evaluated using the thermally buck- 
led/deflected displacements {.4},. Although it is not necessary, the three dynamic load vectors 
will be considered separately for clarity in the dynamic response solution and because practical 
experimental configurations dictate separate loading cases. 
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Figure 6.1 Schematic of acoustically excited panel. 


6.1.1 Acoustic Excitation 

Recall the form for the surface pressure excitation load vector from equation (4.56) 


<m ) 



(6.4) 


The bending components for the system are assembled from the element contributions given in 
equations (4.55), which can be rewritten as 


{i>b(f)} = //>(£. ?/. 0{£«’(£^/)} (l A (6-5) 

where />(£, i/, I ) is the spatial and temporal description of the pressure, {/,„.(£.//)} is the vector 
of out-of-plane interpolation functions, (£.//) are the coordinates in the plane of the structure, 
and A ( is the area of the element. 

In this study, it is assumed that the surface pressure is generated by acoustic excitation. Then 
the pressure loading on the structure consists of the incident pressure, the reflected pressure, and 
the transmitted pressure as shown schematically in Figure 6.1: 

/>(£,//.() = pi (£.//. /) +p r {Z,t),t) -Pttt-V-I) (6-6) 

The pressure loading may also be written as the sum of the blocked pressure and the radiated 
pressure, as suggested by Roussos [52]. The blocked pressure is the pressure on the incident side 
when the structure is considered rigid, i.e., />*,(£, ?/, () = 2 />,(£, 7/,/), and the radiated pressure 
is that due to structure motion only. The acoustic radiation problem can be solved subsequent 
to the forced response solution by assuming that the radiated pressure is negligible compared to 
the blocked pressure. This assumption results in inaccurate predictions only at frequencies near 
the panel fundamental resonance. Thus, the load vector {Pfc(0} is formed using the blocked 
pressure (/>(£. // .t) = ?/. / ) = 2p, (£. //./)) and the transmitted pressure (pressure radiated 

to the free field) can be determined from the resulting structural motion, without the need for 
a fully coupled analysis. 
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Furthermore, it is assumed in this study that the acoustic pressure excitation consists of 
obliquely incident plane waves of arbitrary temporal variation. Then, the general form for the 
acoustic pressure along the direction of propagation given by a Cartesian coordinate - s has the 

form ,r x \ ,, / .\ 

/>)(£.//,/) — s /0 (6 7 ) 

= />(/- h • .v/c) 

where /? is a unit normal aligned with the direction of propagation in the coordinate system 
designated by .r. So, the incident pressure (i.e., c = 0) can be described by 

;),(£,?/. /) = p(t -£ sin 0; cos - >1 sin 0, sin 6,/c) (6.8) 


where 0, and are the polar and azimuthal angles of incidence, respectively. 

The incident pressure is modeled by nodal pressures because of the spatial discretization of 
the finite element approach. Then, the blocked pressure element load vector in equation (6.5) 
becomes 

{/i 6 (/)}= / 2 {!«.({. >/)}/>(/- £ si" 0j cos Oj/c - I/ sin 0j sin o,/c) dA ( (6.9) 

A, 

A more simple form for non-normal blocked acoustic pressure loading will be developed in the 
discrete spectral response section. In the case of normal acoustic pressure excitation, equation 
6.9 simplifies to the following 

{/*(/)}= [ * {LwiClMn dA < 

J ( 6 . 10 ) 

= {ib}p( f ) 

So, the acoustic load vector for the entire system is assembled from element contributions given 
in either of the latter two equations. The system load vector for the normal acoustic excitation 
case takes a simple form 

wm = {i xo < 611 > 


6.1.2 Concentrated Force Excitation 

The concentrated force load vector from equation (4.56) has an analogous form as follows 


{m} 



(6.12) 


The out-of-plane displacement (bending) components for the system are assembled from the 
element contributions given in equations (4.55), which can be rewritten as 


Mt)} = tj 


(6.13) 


'A. 
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where /'(f) is the temporal description of the concentrated force(s), <*>({£} — {£},) is a spatial 
Dirac delta function, and n / is the number of forces within the element area. Note that the 
concentrated forces can be arbitrarily correlated. 

It is assumed that the concentrated force) s) are applied at nodal locations in the model. 
Otherwise, the effect of the concentrated force would be distributed to nodes of the element 
containing the force, which is an undesirable effect. Evaluation of the element area integral 
in equation (6.13) simplifies the form for the element contributions to the concentrated force 
load vector: 

a r 

{/*( 0} = y /'(0U«'(6 -'/.)} 

71 (6.14) 

= ■/?(*) 

Assembly of the element contributions in the system concentrated load vector simply positions 
each concentrated force at a location in the system load vector corresponding to the associated 
out-of-plane degree of freedom. Then, the resulting concentrated load vector for the entire 
system can be written as follows 

{F b (t)} = {/'(/)} (6.15) 

where the only nonzero components in {/'(/)} correspond to forced nodal locations. 


6.1.3 Inertial Excitation 


The inertial loading can be expressed in an analogous load vector form from equation (4.56): 


{(-'(/)} 


U'b(t) 

l 0 


Mb o i r/'V/ji 

0 M m J \ o j 


(6.16) 


The out-of-plane displacement components for the system are assembled from the element 
contributions: 

{</&(0} = -M{«0(O} 

= -K]{l }«*(/) (6.17) 

= {lb}Ug{t) 


where {1} is a vector with unit entries and ii g (t) is the temporal description of the spatially 
uniform base acceleration. So, the system inertial load vector can be written as follows 


{6(fc(0} = {T}u y (t) 


(6.18) 


Solutions to the system of equations given in equation (6.1) can be sought using physical 
or modal degrees of freedom. The modal approach has the advantage of vastly improved 
computational speed without significant loss in accuracy in many cases. The modal approach 
is also quite useful, particularly when comparing with experimental measurements, as powerful 
tools such as experimental modal analysis can be used for incremental validation of the analytical 
model. The modal approach will be used in this study and is discussed briefly here. 
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6.2 Eigen-solution and Modal Coordinates 

Natural frequencies uj r and mode shapes {©, } of the vibration about static equilibrium can 
be obtained from the dynamic eigenvalue problem 

<[.!/]{©, } = [A'] {Or} (6.19) 

A set of uncoupled modal equations with reduced degrees of freedom can thus be obtained from 
equation (6.1) as 

\l\m\ + [2<ru.vj {?(/)} + r^’rJ {</(()} = {/(')} (6.20) 

where r = 1.2 \ (/ , the delimiters [J indicate a diagonal matrix, and [7J is the identity 

matrix, and G is the modal critical damping ratio. The modal force vector {/(/)} is given by 

one of the following expressions as the time dependent loads are considered separately 

{/'(/)} = [©] 7 {/>(/)} acoustic pressure excitation 

= [©] 1 { T}/>( / ) normal acoustic pressure (6 ^ 1 ) 

{/(/)} = [©] 7 { /'( / ) } concentrated force excitation 

{/'(/)} = [o ] 1 { 1 } tig ( / ) inertial/base excitation 

The physical degrees of freedom can be recovered from the truncated modal transformation: 

•v, 

{.4(0} = Yi ■(() = [ < ?]{ f /( / )} \ < -W (6*22) 

r=l 

where [oj is the modal matrix, N, t is the number of modes retained, and A ,/ 0 y is the total number 
of degrees of freedom. 

The modal response, solution to the modal equations of motion (6.20 ), can be determined by 
time domain or frequency domain analysis. In either case, a continuous or discrete approach can 
be employed. In this study, dynamic analysis is done in the frequency domain through continuous 
and discrete techniques. The continuous frequency domain approach is included mainly for 
demonstrative and comparative purposes. The discrete frequency domain approach employed 
in this study is the main spectral solution tool because of the ability to handle complicated 
frequency distributions in a straightforward manner. 

6.3 Continuous Spectral Analysis — Random Response 

The time response of the modal coordinates can be obtained from the convolution 
(Duhamel’s) integral 

I t to 

{(/(/)} =/[/»(/ -t)\{J'(t)} (It = I \h(T)\{f(f - -)} (It (6.23) 

i o « 
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where |7/(/)J is a diagonal matrix of impulse response functions and {/(/)} is the modal force 
vector. The modal response correlation matrix [( ' qq {1 \. f_>)] can be formed from the following 
expected value 

t \ ~ to t > /o 

i:i{<l(h)}{<i(f>)} T } = j f r/dr,)J£[{/(/i-ri )}{/(/, -T>)} l }\h(-,)\<lT l( b, (6.24) 

0 0 

where the expected value in the integrand is the correlation matrix of the excitation 

[( 7/( / i - - T ~ )] • 

The correlation matrix for stationary > excitation becomes a function of the time difference 
only and will be denoted as follows 

£[{/(/ 1 - Ti )}{/(/•.» - rj)} 7 ] = [Cfj(t> - 1 1 + Ti - r > )] (6.25) 

For a load with spectral density s 7/(^.’) or Gjj(f), the Wiener-Khinchine relations for stationary 
random processes gives 

OC 

[Cff(h - /, + r, - r,)] = j rf- 

(6.26) 

= I [Gf f (f)]t ,2Tf{h -' ,+ri - T2] df 

o 


where Gf/U) = for = 2k f (/ > 0). The modal response correlation matrix, 

equation (6.24), then becomes 


OC ti~1 o f 2 ^0 

[('w(ti.h)] = I I [//(nJJr'^'dn J 

— oc 0 0 

OC 

— oc 


(6.27) 


where [7f(u,\/)j is a diagonal matrix of transient frequency response functions. 

For stationary response, let to — * — oc and H(~ . I ) — > H (u:) so the modal displacement 
correlation matrix becomes 

OC 

K’ w (r)]= I r//V)J [%M] ( 6 . 28 ) 

— 'OC' 


where r = t i — t\ and 


H r( u -’ ) 


0 



1 

-f- 2/C r u.V^' 


(6.29) 


no 



is a stationary frequency response function. Assuming mean-square differentiability, Correlations 
of the modal velocity and acceleration responses can be determined by assuming mean-square 
differentiability: 

no 

K *(r)| = f -WMJ 

— 2>C 

and 

OC 

[( ',/</( r)] = I ^ [ // + Mj [• S 7/M] ( '~ T (loJ 


(6.30) 


(6.31) 


where the relations ^ ^ and have been employed. The three modal response 

correlation matrices can be written in the following concise form. 


{[( 


m 


X 

(Ol-KWrU-KWrni = f J"- A) \iru )J [%(.')] r«(.-)J<“ r 


(6.32) 


Note that an arbitrary element (/•, .s) of the modal response correlation matrices has the following 
form 

(E [q r (l )q s (/ + r)]. E[q r (1 + t)]. E[<j, (i)q $ {f + r)]) = L[C qr<if ( T )< ( tj r iA T )' ( </r^( r )] 

no 

—no 

oc 

(6.33) 

Recall the expression for the modal force vector {/(/)} for each load type from equations 
(6.21). Then, the modal force vector correlation and spectral density matrices become 

[Cff{ T )} = [0] T [Cp P ( T )}[(?] => [£//(*;)] = 

[Cff{ T )] = [o} T [Cff’( T )} [<?>] => [>7/M] = [^] 7 [£/'/' M][<£] 

[Cf f (T)} = H T {r}6X, i9 (r){r} T [<(>] => [%M] = [^] r {T}s^){T) T [<g^ 

for acoustic pressure, concentrated force, and base acceleration excitation, respectively. For the 

case of normal acoustic pressure, the corresponding relations above reduce to 

[Cff(r)] = [o] 7 ’{r}r w (r){r} r [ 0 ] => [%M] = [of {r}/>M{r} 7 [6] (6.35) 

The main simplification in these expressions is the presence of a single auto-correlation/power- 
spectral-density function rather than a covariance/spectral-density matrix describing the spatial 
variation of the acoustic load. Note that for the cases of normal acoustic, single concentrated 
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force, or base acceleration excitation, the spectral density function in the integrand of equation 
(6.33) is the same for all elements (r.,s) of the modal response correlation matrices: 


( /'„ [</ r ( / )f/. s ( / +7 )]. h [</ r ( / )</.s(( + 7")]. L [(jr{f + T )] ) — I ^ ^H r (a,’)//. s (a. , )>* + (u.’)( (lu.' 

— oc 
oc 

= f (2 xf) i0 '- A] li;(f)lI,(f )(!"(/)< l2TfT <(f 

it 

( 6 . 36 ) 

Examples of the evaluation of these integrals for determination of the dynamic response quantities 
are given in Appendix F. 

The responses in the physical coordinates can be recovered from the modal expansion, 
equation (6.22), i.e., 

)T 


[C-u(r)] = H[G /(/ (r)](o] J 


( 6 . 37 ) 


Thus, the physical coordinate displacement, velocity, and acceleration correlation matrices 
become 


([(' aa ( t )}. [<”. U (r)].[r /U (r)]) 


OC 

= J J iUA] [o}\ir(^)\ [%(w)] \HU)\{o} T ( ^ r duj 

( 6 . 38 ) 


6.4 Discrete Spectral Analysis 

A discrete spectral analysis approach will now be formulated as an alternative to the 
continuous approach because continuous approaches are limited to rather simple loading cases, 
which are not very realistic for practical applications and not easily generated in an experimental 
situation. In practice, there is always a finite amount of data. Thus, one would employ a Finite 
Fourier Transform to a sampled time history as observed through a time window. Furthermore, 
for discretely sampled time data, a Discrete Fourier Transform (DFT) can be used to estimate 
the Finite Fourier Transforms, etc. A particularly efficient means of performing a DFT is known 
as the Fast Fourier Transform (FFT). Of greater importance is that most experimental data is 
processed via the FFT. A means of using experimentally measured data has obvious beneficial 
implications. 

A discrete spectral analysis can be employed by assuming that the modal force vector, and 
thus the modal response, can be written as the sum of harmonic functions: 

{/(/)} = {/„}(*"•'"' ( 6 . 39 ) 

H=0 

al,d A’. — I 

<«<')}=£<$»}'“■' < 6 - 40 > 

11=0 
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where is the number of discrete frequencies in the transformation and {./„ } and {</» } are the 
complex spectral amplitude vectors, one for each frequency (see Doyle [53]). Recall the modal 
expansion for the physical coordinate response from equation (6.22): 

A2--1 

{Mi)} = Y 

«=° (6.41) 

= Y [<?}{(! 

w=0 


Substitution of equations (6.39) and (6.40) into the modal equations of motion, equation (6.20), 
gives v 

m + w. raCr-vJ + = E (642) 


M=0 


//=(! 


This relation can be written for any one of the discrete frequencies u>„ as 

( r-^j + [/2(, r u.vu.’„j + r^'rj ){</»} = {./»} 


(6.43) 


This relation is simply a set of uncoupled algebraic equations for the modal displacement spectral 
amplitudes. Thus, the modal displacement spectral amplitude vector becomes 

{<ln} = \H(*u)\{f«} (6 ' 44 > 


where the [//(-„ )j is a diagonal matrix of frequency response functions given by 


Hr 


— os~ + / 2( r u: r 


(6.45) 


Recall the expression for the incident acoustic pressure in equation (6.8). In this discrete 
spectral analysis, the incident pressure is written as a sum of harmonic components: 


Pi(M V'l) = p{t — £ sin &i cos <j>i — // sin 0, sin <?, ) 
— l 


= Y 


A T ,£ sin 0 t cos<t>,— k n r/sin t), sin <j> t 


(6.46) 


n=0 


where k„ = lo„/c is the acoustic wave number. Then, the blocked pressure element load vector 
from equation (6.9) becomes 


{/»(<)} = Y 


17 ) 


w=() 
N u > — 1 


A. 


= Y {~n»>}i , » < ^ nt 


(6.47) 
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where h„ = h n sin 0, cos 0 , and A„ = k„ sin 0, sin 0 , and { 7 /,„} contains the integrated spatial 
correlation information of the incident acoustic pressure. The assembled system load vector due 
to acoustic excitation would be written as 

AT. — 1 

{IW)} = Y {M/V"' <6.48) 

»=i) 

The modal force vector is related to the mechanical excitation load vectors from the Galerkin 
weighted average of the modal system (modal system diagonalization), equation (6.21): 

{/(<)} = [o ] 1 {/HO} => {/«} = [of {■ 

{W)} = [o\ T {/'(>)} =* {/„} = [ofifj (6.49) 

{/(0} = [&f {•0'7y(0 =* {./»} = H 1 

where f >„ , and u gn are the complex spectral amplitudes of the acoustic pressure, concentrated 
force, and base acceleration, respectively. Note that, in contrast to the acoustic pressure modal 
force spectral density matrix in the continuous spectral analysis, the expression for the modal 
force spectral amplitude due to acoustic pressure is valid for any angles of incidence. 

If the excitation is deterministic, the dynamic response solution in terms of modal dis- 
placement, velocity, and acceleration time responses can be computed by an inverse FFT from 

equations (6.40) and (6.44): 

AT-) 

{$(')})= Y l ( 6 . 50 ) 

H=0 

where the modal force spectral amplitude has one of the forms given in equations (6.49). Then 
the responses in physical coordinates follow from equation (6.41): 

AT— 1 

({.4(/)}.{.4(/)}. {,!(/)}) = Y R/(^)J {A (6.51) 

«=o 

However, for stochastic excitation, a statistical approach must be used. Additionally, dynamic 
response cross-correlation quantities are of interest for the acoustic radiation discussion to follow 
in the next chapter. Therefore, the necessary relations for stochastic dynamic response will be 
developed next. 

The modal displacement response correlation matrix can be formed from equation (6.40) as 
the following expected value 

[(«(*!' * 2 )] = £[{ 9 (*i )}{</(*-')} 7 ] 

AT— 1 AT—] (6 5°) 

= Y . E £[<&}<«») v -""-”"" 1 

?M = () W = 0 

Substituting from equation (6.44) gives 

AT - 1 AT— 1 

tr ;,)]= £ 5 T r //*(-•... > j ^ [{/,•„}{/., I r ] r //<-., >j 1 1 1 (6.53) 

r/f=() n =0 
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where £[{ /*, } {/„ } 7 ] is a generalized spectral matrix. Note that if the excitation is deterministic, 
the ensemble average indicated by the expected value in equation (6.53) is not necessary. If the 
excitation is stochastic, the necessary ensemble averages can be performed ahead of time. 

For a stationary process, the expected value of the excitation complex spectral amplitudes 
takes a special form [54] 

EK&H/.f] = CKAH/.I’R <«•»> 

where S wn is the Kronecker delta. Then, the stationary modal displacement correlation matrix 
results ( 

[f’„(r)] = £ {<!,,}*)<■“”* (6.55) 

«=0 

where r = f > - l\. The stationary modal velocity and acceleration correlation matrices follow 
from the relations = ~jf anc * rfh ~ Jf as 


A' 2 -i 

[r,„(r)]= Y. 

n=0 

= y Hi } 7 ]-“ 

w=« 


(6.56) 


The full expressions can be written in the following concise form for convenience 


A : oi — 1 


{[(-'qq( T )]* )]- K W T )D = E 4 “ ' \ H* (- » )J E[{ /« } { /» JIT 11 ) J ^ ^ 


w=0 

A^.-l 


(6.57) 


= E 4 0, [^//(^- >] )j ^ 


w=o 


where appropriate quantities from equations (6.49) are employed to evaluate the single-sided, 
stationary modal force spectral matrix [6// (-;„)] for circular frequency 


[f»7/(^ )] — [^] 7 { } f ) nPti {I n } [°] 

[r I - // (^.„)l = [<2.] 3 '{/~>{/i} T '[o] 

\(i / )] = [o] { r« } u (/» ii gn { } [^] 


acoustic pressure 
concentrated forces 
base acceleration 


(6.58) 


Again, the responses in physical coordinates are obtained from the modal expansion 


A’. - 1 




lo} 1 ( t ' c,,T 


»=0 


(6.59) 


These expressions for the modal and physical coordinate response correlation matrices are the 
discrete analogs to those given in equations (6.32) and (6.38). The modal/physical response 
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mean-square values (variances) are determined from the diagonal elements of the corresponding 
correlation matrices in the above equations with r = 0. 

In general, the full expressions in equation (6.57)/(6.59) must be used for calculation of 
the modal/physical responses. However, significant computational savings can be realized 
by recognizing a simplification for an excitation class. Note that for the cases of a single 
concentrated force or multiple contemporaneous forces of the same amplitude, the expression 
for the single-sided, stationary spectral matrix can be simplified: 

[«//<“.•,)] = M T {n/i;/i,)r) , i»] (6.60) 

So, the general form for the spectral matrix for acoustic pressure, single or multiple equivalent 
concentrated forces, or base acceleration can be written as 

[OffM] = <6.61) 


Substituting into equation (6.59) for the physical coordinate correlation matrices results in 

A"„. — 1 

([f ,,!(»]. [r.. u (r)],[r iU (r)])= Y J! )JA) [o}\ir(^m ! {r:}ai 

{l’» } 7 [o][//U'„ )J 
Y w>L 0 ’-' 4) {Hl(a.’„)}{HA(^)} J 




;V U -1 


I ( l^nT 


M= 0 


where 


{hip,,)} = [o][frp„)jH 7 {r;K/x/2 


(6.62) 


(6.63) 


is the conjugated physical frequency response vector to the excitation amplitude determined by 

(l! = (6-64) 


at any discrete frequency os n . 

It is clear that calculation of these physical response correlation matrices has been reduced 
to calculation of a frequency response vector. Mean-square (r = 0) physical responses are 
determined from 

A'*— 1 

K“} = Y 4°- 2 ' 4) |{HaP, i )}|- (6.65) 

n=0 

Any excitation that can be put into the form given in equation (6.61) can take advantage of this 
efficient means of response calculation. Examples of excitations not satisfying this requirement 
include temporally random planar acoustic waves with random incidence, turbulent boundary 
layer loading, multiple non-contemporaneous concentrated forces, etc. 

In the previous sections, forms for the time-dependent load vectors in equation (6.1) 
were derived and a modal approach for continuous and discrete frequency domain analysis 
of the vibration response of panel-type structures subjected to those loads was developed. 
The formulation and implementation will now be validated through comparison with classical 
solutions and experimental measurements in the next two sections. 
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f(t) 


t Ug(t) 




j Ug(t) 


Figure 6.2 Schematic of simply-supported beam subject to planar 
acoustic, concentrated force, and base acceleration excitation. 


Table 6.1 Material properties for aluminum specimens. 


Young’s Modulus 
Shear Modulus 
Poisson's Ratio 
Mass Density 
Coeff. of Thermal Exp. 


68.95 GPa 10 Msi 

26.55 GPa 3.85 Msi 

0.3 

2.767e+3 kg/m 3 2.588e-4 lbf s 2 /in 4 

22.3 iu/° C 12.4 /n/°F 


6.5 Analytical Validation 

In this section, comparisons will be made with classical analysis for the cases of a simply- 
supported aluminum beam and a simply-supported aluminum panel to benchmark the analysis. 
All three dynamic, mechanical loads are considered; acoustic pressure, inertial, and concentrated 
force. Both the continuous and discrete frequency domain approaches are employed. 


6.5.1 Simply-Supported Aluminum Beam 

Consider a simply-supported aluminum beam shown schematically in Figure 6.2 with 
dimensions of 0.4572x0.0254x0.0023 m (18x1x0.090 inches). The material properties for 
the specimen are given in Table 6.1. It is assumed that the various loads are mutually exclusive, 
and results will be presented as such. Results for continuous and discrete spectral analysis from 
the present formulation will be presented. A classical solution for the random response of this 
structure subjected to each of the loads was derived for comparison. 

The natural frequencies from the classical solution are governed by the equation 


2 xfn — w’jj 




(6.66) 


where m is the mass per unit length and a is the length of the beam. The beam was modeled by 
a finite element mesh of 36x2 elements. A comparison of the natural frequencies as predicted 
by the classical analysis and the finite element analysis are shown in Table 6.2. Excellent 
agreement between the two solutions is apparent. 

The three excitation types are each modeled by white noise with spectral amplitudes which 
give consistent loading. For example, a lg RMS loading over a 10-500 Hz bandwidth is 
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Table 6.2 Comparison of aluminum 0.4572x0.0254x0.0023 m (18x1x0.090 inch) 
simply-supported beam modal frequencies as predicted by classical and finite element analyses. 


Mode Number 

Classical, Hz 

FEM, Hz 

1 

24.76 

24.76 

2 

99.04 

99.11 

3 

222.84 

223.20 

4 

396.15 

397.20 


Table 6.3 Comparison of mid-span RMS responses for an aluminum 0.4572x0.0254x0.0023 m 
(18x1x0.090 inch) simply-supported beam as predicted by classical and finite element analyses. 


Analysis 

RMS Displ. (in) 

RMS Vel. (in/s) 

RMS Accel. 


m 

in 

m/s 

in/s 

m/s 2 

in/s 2 


Base Acceleration 

Classical 

7.26e-4 

0.0286 

0.1138 

4.482 

— 

— 

FE Continuous 

7.21e-4 

0.0284 

0.1135 

4.469 

— 

— 

FE Discrete 

7.21e-4 

0.0284 

0.1135 

4.469 

26.52 

1044 


Acoustic Pressure 

Classical 

7.26e-4 

0.0286 

0.1138 

4.482 

— 

— 

FE Continuous 

7.21e-4 

0.0284 

0.1135 

4.469 

— 

— 

FE Discrete 

7.21e-4 

0.0284 

0.1135 

4.469 

26.52 

1044 


Concentrated Force 

Classical 

5.74e-4 

0.0226 

0.1068 

4.205 

— 

— 

FE Continuous 

5.69e-4 

0.0224 

0.1063 

4.184 

— 

— 

FE Discrete 

5.69e-4 

0.0224 

0.1063 

4.184 

83.64 

3293 


considered to be a representative inertial loading. This overall load corresponds to a spectral 
level of 0. 1966 (m/s 2 ) 2 /Hz (304.704 (in/s 2 ) 2 /Hz). One can arrive at an equivalent surface pressure 
excitation by equating the product of the acceleration level and mass of the beam, which results 
in a load of 0.721 N (0.162 lbf), with the product of the unknown surface pressure level and the 
surface area of the beam. The resulting equivalent acoustic pressure excitation has an overall 
value of 62.1 Pa (0.009 psi) RMS or a spectral level of 7.859 Pa 2 /Hz (1.65307e-7 (psi) 2 /Hz) over 
a 10-500 Hz bandwidth. Similarly, a concentrated force of 0.721 N (0.162 lbf) RMS or a spectral 
level of 1.0596e-3 N 2 /Hz (5.35595e-5 lbf 2 /Hz) over 10-500 Hz gives a consistent load when 
applied to the beam at mid-span. To make things more interesting, however, the concentrated 
force was applied at node 71 of 1 1 1 nodes with coordinates of (15, 0.5) inches. Comparisons of 
the classical and finite element predictions of the RMS mid-span responses are shown in Table 
6.3. A value of 0.02 was used in all cases for the modal critical damping ratio to desensitize the 
solutions to frequency resolution. The classical and FE continuous spectral analysis predictions 
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were obtained using the assumption of white noise excitation with a summation over the first four 
modes. FE discrete spectral analysis predictions were obtained with a 10-500 Hz band-limited 
white noise model, a 0.5 Hz frequency resolution, and summation over the first four modes. The 
classical and finite element solutions are in excellent agreement in all cases. 


6.5.2 Simply-Supported Aluminum Panel 

Consider a simply-supported aluminum panel shown schematically in Figure 6.3 with 
dimensions of 0.3556x0.254x0.001 m (14x10x0.040 inches). The material properties for the 
specimen are the same as those given in Table 6.1 for the beam case. Again, it is assumed that 
the various loads are mutually exclusive. A classical solution for the random response of this 
structure subjected to each of the loads was derived for comparison. 



Figure 6.3 Schematic of simply-supported panel subjected to planar 
acoustic, concentrated force, and base acceleration excitation. 


The natural frequencies from the classical solution are governed by the equation 


‘2/t/nm — 



(6.67) 


where p is the mass per unit area, a and b are the panel length and width, and I) is the flexural 
rigidity given by 


I) = 


Eh* 

I2( 1 - ) 


(6.68) 
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Table 6.4 Comparison of aluminum 0.3556x0.254x0.001 m (14x10x0.040 inch) 
simply-supported panel modal frequencies as predicted by classical and finite element analyses. 


Mode Number 

x-index 

y-index 

Classical, Hz 

FEM, Hz 

1 

1 

1 

56.4 

56.4 

2 

2 

1 

113.7 

113.7 

3 

1 

2 

168.6 

168.6 

4 

3 

i 

209.0 

209.0 

5 

2 

2 

225.8 

225.8 

6 

3 

2 

321.1 

321.2 

7 

4 

i 

342.5 

342.6 

8 

1 

3 

355.4 

355.6 

9 

2 

3 

412.7 

412.8 

10 

4 

2 

454.6 

454.8 


The panel was modeled by a finite element mesh of 14x10 elements. A comparison of the natural 
frequencies as predicted by the classical analysis and the finite element analysis are shown in 
Table 6.4. Negligible discrepancy exists between the two solutions. 

The inertial excitation level used in the previous example will also be used in this case, i.e., 
lg RMS or 0.1966 (m/s 2 ) 2 /Hz (304.704 (in/s 2 ) 2 /Hz) over a 10-500 Hz bandwidth. The overall 
load of this inertial excitation is 2.49 N (0.56 lbf). The equivalent surface pressure excitation 
for this case is 27.58 Pa (0.004 psi) RMS or a spectral level of 1.5524 Pa 2 /Hz (3.265327e-8 
(psi) 2 /Hz) over a 10-500 Hz bandwidth. Finally, a concentrated force of 2.49 N (0.56 lbf) RMS 
or a spectral level of 12.66e-3 N 2 /Hz (6.4000409e— 4 lbf 2 /Hz) over 10-500 Hz is applied at node 
116 of 165 nodes with coordinates of x=0.254, y=0. 178 m (x=10, y—1 inches). Comparisons of 
the classical and finite element predictions of RMS center-point responses are shown in Table 
6.5. A value of 0.02 was used in all cases for the modal critical damping ratio. The classical 
and FE continuous spectral analysis predictions were obtained using the assumption of band- 
unlimited white noise with a summation over the first ten modes. FE discrete spectral analysis 
predictions were obtained with a 10-500 Hz band-limited white noise model, a 0.5 Hz frequency 
resolution, and summation over the first ten modes. It is clear that excellent comparisons are 
achieved in all cases. 

Comparisons of the displacement power spectral densities PSDs at the coordinates x=0. 102, 
y= 0.076 m (x=4, y=3 inches) predicted by the classical and finite element analyses for the three 
load cases above are presented in Figures 6.4—6. 6 . In all cases, the classical and finite element 
PSDs are indistinguishable. Note the presence of the non-symmetric modes in the concentrated 
force response case. PSDs for a random acoustic pressure case with angles of incidence (#j=60°, 
e>j=45°) are shown in Figure 6.7. This case further demonstrates the excitation of the non- 
symmetric modes and will be useful for reference in Chapter 7. 


120 





















Table 6.5 Comparison of center-point RMS responses for an aluminum 0.3556x0.254x0.001 m 
(14x10x0.040 inch) simply-supported panel as predicted by classical and finite element analyses. 


Analysis 

RMS Displ. 

RMS Vel. 

RMS Accel. 


m 

in 

m/s 

in/s 

m/s 2 

in/s s 


Base Acceleration 

Classical 

2.692e-3 

0.0106 

0.098 

3.845 

— 

— 

FE Continuous 

2.692e-3 

0.0106 

0.098 

3.845 

— 

— 

FE Discrete 

2.692e-3 

0.0106 

0.098 

3.844 

49.96 

1967 


Acoustic Pressure 

Classical 

2.692e-3 

0.0106 

0.098 

3.845 

— 

— 

FE Continuous 

2.692e-3 

0.0106 

0.098 

3.845 

— 

— 

FE Discrete 

2.692e-3 

0.0106 

0.098 

3.844 

49.96 

1967 


Concentrated Force 

Classical 

4.216e-4 

0.0166 

0.157 

6.167 

— 

— 

FE Continuous 

4.216e-4 

0.0166 

0.157 

6.167 

— 

— 

FE Discrete 

4.216e-4 

0.0166 

0.157 

6.167 

92.63 

3647 




Displacement PSD in 2 /Hz, Af=0.5 Hz 



Figure 6.4 Classical and finite element displacement response PSDs at x=0.102, 
y= 0.076 m (x=4, y=3 inches) versus frequency for inertial loading case. 
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Figure 6.5 Classical and finite element displacement response PSDs at x=0.102, y- 0.076 
m (x=4, y=3 inches) versus frequency for normal acoustic pressure loading case. 
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Figure 6.6 Classical and finite element displacement response PSDs at x=0.102, y 
0.076 m (x=4, y=3 inches) versus frequency for concentrated force loading case. 
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Displacement PSD in 2 /Hz, Af=0.5 Hz 

— — k. L 



Frequency, Hz 

Figure 6.7 Classical and finite element displacement response PSDs at x=0.102, y= 0.076 m 
(x=4, y=3 inches) versus frequency for ^=60°, Oj=45° acoustic pressure loading case. 
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It is physically intuitive as to which mechanical loads will excite the non-symmetric modes. 
However, some observations about the numerical modeling of such loads is worth mentioning 
here. In the discussion that follows, attention is limited to the vector and matrix components 
associated with the out-of-plane degrees of freedom only as all other components of the flexural 
eigenvectors are zero. 

In general, a non-symmetric mode will not be excited by any excitation for which the 
corresponding mode shape is orthogonal with respect to the excitation spectral density matrix: 

{o}f[f.V)]{o} ( = 0 (6.69) 

Normal acoustic pressure and inertial loading are special cases of this situation as they result in 
a real, uniform (full with all entries identical) excitation spectral density matrix. Obviously, the 
product of any antisymmetric modal vector transpose and a uniform matrix is the null vector 
transpose. The resulting modal force spectral density matrix is sparse with nonzero entries only 
at indices corresponding to the symmetric modes. 

A concentrated force excitation is modeled by a null spectral density matrix with the 
exception of a single nonzero entry located by the index of the degree of freedom at the point 
of force application. The product of any modal vector with this matrix is nonzero except for 
the case when the nonzero matrix index corresponds to a zero in the modal vector. Hence, all 
modes are excited except those for which the point of force application falls on a nodal line. 
In this case, the modal force spectral density matrix is real and full, only the modes with nodal 
lines containing the point of force application are omitted. 

Obliquely incident acoustic pressure excitation gives a complex, full but non-uniform spectral 
density matrix. The product of any modal vector with this matrix is nonzero except for cases 
when the angles of incidence direct the propagating pressure along an axis of the structure. In 
this special case, portions of the excitation spectral density matrix are uniform and these portions 
have indices corresponding to the antisymmetric nature of even-odd modes (e.g., dj=60°, Oj=0° 
does not excite the *,2 modes). So, with the exception of the modes excluded due to the latter 
situation, the modal force spectral density matrix is complex and full. One final note should be 
made regarding the extent to which the antisymmetric modes are excited. Modes occurring at 
lower frequencies are excited more weakly by excitations of this type because the wavelength 
of the propagating pressure may be long relative to the extent of the structure. This causes the 
incident pressure to appear to be close to uniform. This is exemplified in the case presented in 
Figure 6.7, where the pressure wavelengths range from 34.29 to 0.686 m (112.5 to 2.25 feet) 
at 10 Hz and 500 Hz, respectively. In this case, all antisymmetric modes are excited weakly 
as is evident in the figure. 

The finite element analysis has been shown to produce very accurate results in comparison 
with classical analysis for all three loading types for simply-supported beam and panel structures. 
Now, attention will be turned to experimental configurations for validation of the analytical tool 
with experimental measurement. 

6.6 Experimental Validation 

Comparisons will be made with experimental measurements for the cases of a clamped- 
clamped aluminum beam and a clamped-clamped glass-epoxy beam to assess model and analysis 
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tool accuracy incrementally, prior to making comparisons with measurements from a clamped- 
clamped SMAHC beam. Results will be shown for the SMAHC beam at several temperatures 
indicative of various states in the beam thermomechanical response. Finally, comparisons will 
be made for the case of a fully-clamped aluminum panel as an introduction to the next chapter. 

All of the experimental beam vibration results presented in this study result from base 
acceleration (inertial) excitation. A picture of the overall apparatus is shown in Figure 6.8. 
The shaker system consists of a large electromechanical shaker configured horizontally and 
attached to a magnesium slip table, which oscillates in a horizontal plane during operation. The 
specimen is mounted in the fixture with a vertical orientation to the width-direction in order to 
avoid complications due to gravity. The excitation acceleration is measured by accelerometers 
(PCB model T352M92) placed at each end of the specimen fixture. The normal velocity and 
acceleration responses of the beam are measured with a laser Doppler velocimeter (Ometron, 
VPI) and accelerometers (Endevco model 2250A-10). It is obvious that the beam response 
transducers measure the absolute motion of the beam specimens (including the base motion). 
Most structural dynamic prediction tools, in fact all of those used in this study, are formulated 
to provide relative motion response for such configurations. Therefore, it is imperative that 
experimental procedures allow for the removal of the base motion contributions. The importance 
is not so much for comparison of response amplitude as the base motion is typically much less 
than that of the flexible structure, but for comparison of frequency response or response spectra. 



Figure 6.8 Overall configuration for base acceleration tests on beam specimens. 


To allow the most flexibility in the procedure of base motion removal, a redundant accelerom- 
eter was attached to each end of the support fixture and connected to signal conditioning with 
hardware integration. The redundant reference accelerometer signals were integrated once to 
result in velocity measurements for subtraction from the laser vibrometer output. Subsequently, 
the base motion can be removed from the acceleration or velocity responses in either the time 
domain, by simple subtraction, or in the frequency domain as follows. This discussion will 
focus on the accelerometer measurements, but it equally valid for any response quantity. The 
relative acceleration can be written as 

a = o , — a r (6.70) 
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where the specimen and fixture accelerations have been denoted by«* and a f . For stationary 
stochastic response, the covariance can be obtained from the following expectation 

Fi[a(t)(i{l - r )] = £’[(«.,(/) - «f(t)){(is(t ~ t ) - a j{t - r))] (6.71) 

This can be written in the frequency domain as 

£ > uh( u -') "F Ou fa f (u. - ) (6.72) 

where the o are power and cross spectral density functions. Thus, it is clear that frequency 
domain measurements (i.e., spectral matrix) can be post-processed to obtain the relative motion 
spectra. In this study, the time data was collected and post-processed to remove the base motion 
and generate relative response spectra. 

6.6.1 Clamped Aluminum Beam 

A picture of the apparatus for the clamped-clamped aluminum beam configuration is shown 
in Figure 6.9. The fixture is the same as that described in Chapter 5, but the electrical connection 
blocks and electrical/thermal insulation layers have been removed. The overall beam dimensions 
are 0.5588x0.0254x0.0023 m (22x1x0.09 inches) with an unsupported length of 0.4572 m (18 
inches). The overall beam mass was measured to be 89.79e-3 kg (0.198 lbf). The elastic 
properties of the material were also measured and are given in Table 6.6. 



Figure 6.9 Aluminum beam specimen mounted in base acceleration test configuration. 


The beam was excited by random base acceleration over a 20-400 Hz bandwidth with a 
RMS value of 0.25 g. The beam response was measured by the laser vibrometer and two 
accelerometers at locations 0. 10 m (4 inches) and 0.25 m (10 inches) from right grip, indicated 
by position numbers 1 and 2 in Figure 6.9. Time records with an overall length of 200 seconds 
were collected with a sampling rate of 1024 Hz to allow 50-frame averages with a bandwidth 
of 0-400 Hz and a frequency resolution of 0.25 Hz. The time data were post-processed to 
construct relative response quantities, frequency response functions between the beam response 
transducers and the fixture transducers, and power spectral densities (PSDs) of the acceleration 
input and relative beam responses. 
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Table 6.6 Material properties for aluminum beam specimen in dynamic response test. 


Young’s Modulus 
Shear Modulus 
Poisson’s Ratio 
Density 

Mode Number 
1 

2 

3 


74.47 GPa 10.8 Msi 

27.99 GPa 4.06 Msi 

0.33 

2.767e+3 kg/m 3 2.588e-4 lbf«s/in 4 

Modal Critical Damping Ratio 
0.0025 
0.001 
0.00086 


The beam was modeled with a 36x2 mesh of the 24 degree-of-freedom rectangular plate 
elements described in Chapter 4. Modal damping estimates from measured frequency response 
functions, see Table 6.6, were used in the predictions. The beam grips were torqued to 13.56 
N»m (120 in«lbf), which gave a measured fundamental frequency of 55.5 Hz. The fundamental 
frequency with perfectly clamped boundaries was predicted to be 58.4 Hz. The boundary 
conditions were modeled as clamped with the exception of torsional springs, with a spring 
rate of 113 N«m/rad (1000 in*lbf/rad) applied to the longitudinal-slope degrees of freedom to 
match the fundamental frequency. Only two symmetric modes were expected in the excitation 
bandwidth, but a total of 10 modes was used in the predictions to account for the residual effects 
of the higher order modes. A uniform modal critical damping ratio of 0. 1 % was used for modes 
4-10. The influence of the accelerometers was modeled as lumped masses (~0.4e-3 kg, 2.283e-6 
lbf»s 2 /in 4 ) at appropriate nodes. The input acceleration PSDs, as measured by one of the control 
accelerometers, was taken as input for the predictions. 

Comparisons of predicted and measured displacement PSDs for measurement locations 1 
and 2 are shown in Figures 6.10 and 6.1 1. Note that the “measured” displacement PSDs were 
generated by integrating the acceleration responses in the frequency domain. It can be seen that 
excellent agreement is achieved, except in overall RMS response level. The discrepancy in RMS 
level is due to the fact that the damping is so light in this case that the frequency resolution of 
0.25 Hz is insufficient to adequately capture the peak values in the FRFs or spectra. Of course, 
the spectral analysis conserves the total response energy by spill-over into adjacent frequency 
bins. This spill-over results in slightly broader response peaks with lower amplitude. Curve 
fitting techniques for estimation of modal parameters misinterpret this as higher damping values. 
The main culprit in this case is the damping value at the fundamental frequency, where the 
response peak is extremely sharp. A damping value of approximately 0. 1 % would probably be 
more representative of the actual damping value in the first mode because it is consistent with 
that at mode 3. With this change in damping, the displacement response RMS values at the two 
measurement locations are 1.24e-4 m and 2.54e-4 m, respectively, in nearly perfect agreement 
with measurement. 
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Displacement PSD, m 2 /Hz 



Frequency, Hz 

Figure 6.11 Measured and predicted displacement PSD 
comparison for position 2 on the aluminum beam (Af=0.25 Hz) 
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Table 6.7 Material properties for glass-epoxy beam specimen in dynamic response test. 


Ei 

E 2 

Gj2 

/ 7 12 

P 

Mode Number 
1 

2 

3 

4 


49.30 GPa 
20.00 GPa 
9.65 GPa 
0.29 

2.03 le+3 kg/m 


7.15 Msi 
2.90 Msi 
1.40 Msi 
0.29 

0. 19e-3 lbf»s/in 
Modal Critical Damping Ratio 


4 


0.0055 

0.005 

0.0048 

0.005 


6.6.2 Clamped Glass-Epoxy Beam 

The performance of the FE analysis will now be tested against a baseline composite beam. 
A glass-epoxy beam with dimensions of 0.5588x0.0254x0.0018 m (22x1x0.072 inches) with an 
unsupported length of 0.4572 m ( 1 8 inches) was tested in the same experimental configuration as 
that for the aluminum beam. The overall beam mass was measured to be 53.9e-3 kg (0.1 19 lbf). 
The elastic properties of the material system were described in Chapter 3. The material properties 
at ambient temperature were used to generate the predicted responses. Those properties are 
reiterated in Table 6.7 for convenience. 

The beam excitation and data acquisition configuration was the same as that for the aluminum 
beam with the exception of the excitation bandwidth, which was changed to 10-400 Hz because 
of the lower fundamental frequency in this case. The overall excitation level of 0.25 gRMS was 
maintained. The measurement locations will again be referred to as positions 1 and 2. Time data 
were collected and post-processed to construct relative response quantities, frequency response 
functions between the beam response transducers and the fixture transducers, and power spectral 
densities (PSDs) of the acceleration input and relative beam responses. 

An analogous FE model was also developed for this case. Modal damping values were 
estimated from the measured frequency response functions, see Table 6.7, and used in the 
predictions. In this case, a torque of 13.56 N»m (120 in*lbf) on the beam grips resulted in 
a measured fundamental frequency of 33.5 Hz. The fundamental frequency with perfectly 
clamped boundaries was predicted to be 35 Hz. Torsional springs with a rate of 19.21 N*m/rad 
(170 in»lbf/rad) were applied to match the fundamental frequency. Again, the input acceleration 
PSDs, as measured by one of the control accelerometers, were taken as input for the predictions. 

Predicted and measured displacement PSDs for measurement locations 1 and 2 are shown 
in Figures 6.12 and 6.13. Again, excellent agreement is achieved, including the RMS response 
levels. There is a slight over-prediction of the second symmetric modal frequency, but the 
comparison is within experimental uncertainty. 
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Displacement PSD, itvVHz 



Frequency, Hz 

Figure 6.12 Measured and predicted displacement PSD comparison 
for position 1 on the glass-epoxy beam (Af=0.25 Hz). 
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Displacement PSD, m 2 /Hz 






Figure 6.14 Front view of the SMAHC beam specimen 
mounted in clamping fixture for dynamic testing. 


The dynamic analysis has been validated by comparison with the measured response of 
conventional structures. Attention will now turn to demonstration of the dynamic response 
performance of a SMAHC beam specimen and validation of the formulation by comparison of 
the measured and predicted responses. 

6.6.3 Clamped SMAHC Beam 

The SMAHC beam specimen mounted in the support fixture and configured for base 
acceleration testing is shown in Figure 6.14. The configuration is essentially the same as 
that described in the Experimental Validation section of Chapter 5. The beam is supported 
by the mechanical grips indicated in the figure by “m”. The grips designated as e are 
electrical connections for the Nitinol leads, but also provide mechanical restraint. The electrical 
connections are isolated from the fixture by a layer of fibrous ceramic insulation. The beam is also 
thermally isolated from the fixture by a layer of this insulation on both sides. The beam is heated 
by DC electrical current controlled by a thermocouple measurement located approximately on 
the beam centerline and 0.0127 m (0.5 in) from the right mechanical grip, indicated in the figure 
by “tc”. The thermal controller consists of a process controller (Omega model CN77353-PV), 
and a DC power supply (HP model 6652A). Thermocouples also monitor the temperature at 
each end of the fixture, indicated in Figure 6.14, and one measures the ambient air temperature. 

A rear view of the beam/fixture assembly is given in Figure 6.15, which shows additional 
detail on the dynamic experimental setup. The placement of the control thermocouple is more 
apparent in this view as is the presence of the electrical/thermal insulation. The base acceleration 
excitation was controlled with the average of the two control accelerometers. The “integrated 
accelerometers” were utilized in the calculation of relative velocity response from the laser 
vibrometer data, as described previously. The positions of the two response accelerometers 
are consistent with the descriptions of the previous tests and are again designated by position 
numbers 1 and 2. 

The beam was excited by random base acceleration over a 10-400 Hz bandwidth with a 
RMS value of 0.25 g. Time records with an overall length of 200 seconds were collected with 


Figure 6.15 Rear view of the SMAHC beam specimen 
mounted in clamping fixture for dynamic testing. 


a sampling rate of 1024 Hz to allow 50-frame averages with a bandwidth of 0-400 Hz and 
a frequency resolution of 0.25 Hz. The time data were post-processed to construct relative 
response quantities, frequency response functions between the beam response transducers and 
the fixture transducers, and power spectral densities (PSDs) of the acceleration input and relative 
beam responses. Dynamic data were collected at ambient and at elevated temperatures from 
32.22°C (90°F) to 121.1 1°C (250°F) in 5.56°C (10°F) increments. The beam mechanical and 
ribbon electrical grips were torqued to 13.56 N*m (120 in»lbf) just prior to dynamic excitation 
to minimize effects of changing environmental conditions. The fixture temperature increased by 
approximately 2.22°C (4°F) during the course of data collection, which was in almost direct 
proportion to the increase in ambient air temperature. 

First, some general observations will be made concerning the experimental behavior and 
results. Because of the material system characteristics and test configuration, the beam exhibited 
a thermal post-buckling deflection over range of temperatures, as described in detail in Chapter 
5. The beam exhibited some interesting dynamic behavior (such as intermittent snap-through) 
within this temperature range, particularly at the point of loss of stability in the flat configuration 
and then at the point of regaining stability in the flat configuration, as expected. No attempt will 
be made to predict the dynamic behavior within this range, although it is possible with the present 
analysis for the portions where the dynamic response is geometrically linear about a nonlinear 
static deflection. Although the beam was returned to a flat configuration at approximately 
54.44°C (130°F), the presented results will focus on higher temperatures. The forced vibration 
amplitude was visually noticeable in the inactive state (ambient temperature), but was virtually 
undetectable at temperatures above 71.1 1°C (160°F). No evidence of delamination or other flaws 
were detected subsequent to testing, despite the large thermal post-buckling deflection and rather 
high temperature developed during the tests. 

Other observations can be made from a summary of RMS displacement responses and modal 
parameter estimates shown as a function of temperature in Table 6.8. Temperatures for which the 
beam was buckled are excluded from the table. The RMS displacements correspond to location 
2 and were calculated in the frequency domain from the corresponding accelerometer PSDs. It 
can be seen that, although a majority of the stiffening effect was achieved at approximately 
93.33°C (200°F), performance remains good far beyond that temperature with a maximum 



Table 6.8 Summary of RMS measured displacements at position 2 
on the SMAHC beam and modal parameters versus temperature. 


Temp., °C 

RMS Displ., m 

fi.Hz 

Cl, % 

f 3 , Hz 

C3, % 

23 

0.322e-3 

26.4 

0.89 

160.3 

0.95 

— 

— 

— 

— 

— 

— 

60 

0. 15e-3 

83.1 

0.10 

279.3 

0.45 

65.56 

0.11e-3 

94.4 

0.13 

309.4 


71.11 

0.05e-3 

107.0 

0.50 

342.7 

0.16 

76.67 

0.063e-3 

118.3 

0.19 

372.5 

0.28 

82.22 

0.068e-3 

126.7 

0.15 

— 

— 

87.78 

0.063e-3 

130.9 

0.14 

— 

— 

93.33 

0.062e-3 

134.0 

0.14 

— 

— 

98.89 

0.066e-3 

136.0 

0.11 

— 

— 

104.4 

0.067e-3 

137.6 

0.11 

— 

— 

110.0 

0.073e-3 

139.2 

0.10 

— 

— 

115.6 

0.072e-3 

138.7 

0.10 

— 

— 

121.1 

0.067e-3 

140.0 

0.09 

— 

— 


increase in the fundamental frequency of 430% at 121.1 1°C (250°F). The variation of the modal 
frequencies with temperature is more easily seen in Figure 6.16. There, it can be seen that the 
fundamental frequency starts to reach an asymptotic value at approximately 93.33°C (200° F) 
and mode three escapes the excitation bandwidth at temperatures above 82.22°C (180°F). Note 
the peculiar changes in the mode— 1 damping near 71.1 1°C (160°F). This effect is reflected in the 
measured displacement response RMS values and can be seen in Figure 6.17, which shows the 
measured displacement response PSD at ambient and three elevated temperatures. The four PSDs 
correspond to the filled symbols in Figure 6.16. It can be seen that the maximum attenuation in 
the RMS displacement response is achieved at 71.1 1°C (160°F), where there is maximum peak 
response reduction due to higher damping, although the stiffening effect continues to increase^ at 
higher temperatures. One possible explanation is that there may be an optimal state near 71.11 C 
(160° F) where there is adequate austenite and dynamic stress and low enough temperature to 
allow enhanced damping by stress induced martensitic (SIM) transformation hysteresis. At 
higher temperatures, the dynamic stresses are insufficient to drive the SIM transformation as it 
becomes increasingly difficult to induce. 

Recall that the beam dimensions are 0.5588x0.0254x0.0019 m (22x1x0.078 in) with a 
lamination of (45/0 sma /-45/90) 2s . Five ribbon widths are embedded within each 0° layer. Each 
ribbon has a cross section of 2.29e-3xl.52e-4 m (0.09x0.006 inches). This configuration was 
modeled as a Nitinol volume fraction of 0.5538 in the 0° layers, for an overall volume fraction 
of 13.8%. The mass densities of the glass-epoxy and Nitinol ribbon (prestrained 4%) were 
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Figure 6.16 Measured mode 1 and mode 3 natural 
frequencies of the SMAHC beam versus temperature. 
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Displacement PSD, m 2 /Hz 


Measurement, 74°F 
Measurement, 1 60°F 
Measurement, 200°F 
Measurement, 250°F 







200 

Frequency, Hz 


Figure 6.17 Measured displacement PSDs for position 2 on the 
SMAHC beam at ambient and elevated temperatures (Af=0.25 Hz). 


139 



found to be 0. 19e-3 lbf*s 2 /in 4 and 0.5349e-3 lbf»s 2 /in 4 . With the beam/ribbon dimensions, the 
trimmed ribbon length of 0.66 m (26 inches) and the measured mass densities, the calculated 
mass of the beam is 75.5e-3 kg (0.166 lbf). This is within 1% of the measured overall mass 
of the SMAHC beam, which was found to be 74.92e-3 kg (0.165 lbf). The thermomechanical 
properties of the constituent materials were described in Chapter 3 and are given in Tables 3.5 
and 3.8. Recall that the composite cure cycle constitutes thermal cycle number one (1) on the 
Nitinol ribbon. This was the first test of the SMAHC beam specimen, which corresponds to 
thermal cycle number two (2). The appropriate recovery stress data from Table 3.8 were used 
in the model to predict the dynamic response. 

The element properties for a 32x2-element-mesh model of this beam were evaluated as a 
function of temperature from these material properties. Modal damping estimates were taken 
from the experimental data in Table 6.8. A total of ten modes was used in the predictions in all 
cases to account for residual effects. The damping was assumed to be 0.1% on all antisymmetric 
modes and modes out of bandwidth. In this case, the torque of 13.56 N«m (120 in*lbf) on the 
beam grips resulted in a measured fundamental frequency of 26.4 Hz at ambient temperature. 
The fundamental frequency with perfectly clamped boundaries was predicted to be 37.6 Hz. 
The boundary conditions in this test are not as close to clamped conditions as in the previous 
tests because of the insertion of the thermal insulating material between the beam and grip 
surfaces. Torsional springs with a rate of 1.98 N» m/rad (17.5 in»lbf/rad) were applied to match 
the fundamental frequency. The influence of the accelerometers was modeled as lumped masses 
(~0.4e-3 kg, 2.283e-6 lbf*s 2 /in 4 ) at appropriate nodes. Again, the input acceleration PSDs, as 
measured by one of the control accelerometers, was taken as input for the predictions. 

Measured and predicted displacement PSDs for measurement location 2 on the SMAHC 
beam at four temperatures are shown in Figures 6.18-6.21. The figures show results for 

ambient temperature, 71.1 1°C (160°F), 93.33°C (200°F) and 121. 1°C (250°F), respectively, 
and are representative of the comparisons at all temperatures. The correlation at ambient 
temperature is excellent with less than 9% difference between measured (3.23 mRMS, 0.0127 
inRMS) and predicted (2.96 mRMS, 0.01 16 inRMS) RMS response level. The comparisons at 
elevated temperature show some discrepancy, but are in excellent qualitative agreement. The 
predicted responses consistently show slightly higher modal frequencies and lower RMS levels 
than the corresponding experimental measurements. This trend is due to an over-prediction of 
the recovery force stiffening effect and is attributable to boundary condition modeling issues 
described in the Experimental Validation section of Chapter 5. 

Recall from that discussion that the in-plane boundaries of the physical system are somewhat 
compliant due to the layer of insulation between the mechanical and electrical grips. The 
boundary compliance partially relieved the recovery stress and resulted in a reduced stiffening 
effect. This effect was captured in the model by applying axial springs with a spring rate of 
1.97e+6 N/m (11,250 lbf/in) to the in-plane degrees of freedom at the mechanical grips. The 
effect of fixture expansion with temperature was also modeled by applying appropriate enforced 
displacements at the beam ends. New predicted responses were made with these changes. The 
measured and predicted displacement PSDs for measurement location 2 on the SMAHC beam 
are shown in Figures 6.22-6.25. Excellent agreement is achieved between the predicted and 
measured responses at all temperatures. 
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Figure 6.18 Measured and predicted displacement PSD comparison for 
position 2 on the SMAHC beam at ambient temperature (Af=0.25 Hz). 
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Figure 6.19 Measured and predicted displacement PSD comparison 
for position 2 on the SMAHC beam at 71.1°C (160°F, Af=0.25 Hz). 
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Figure 6.20 Measured and predicted displacement PSD comparison 
for position 2 on the SMAHC beam at 93.3°C (200°F, Af=0.25 Hz). 
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Figure 6.22 Measured and predicted displacement PSD comparison for 
the SMAHC beam at ambient temperature with boundary effects (Af=0.25 Hz). 
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Figure 6.24 Measured and predicted displacement PSD comparison for 
the SMAHC beam at 93.3°C with boundary effects (200°F, Af=0.25 Hz). 
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Receiving room Source room 



Figure 6.26 Transmission loss apparatus schematic. 


The formulation and solution procedures have now been fully exercised and validated against 
a variety of experimental results, most importantly with the static and dynamic response of a 
SMAHC beam. Further work was done to validate the acoustic radiation model. Experimental 
validation was attempted for the case of a clamped aluminum panel. Dynamic responses of the 
panel due to acoustic excitation were measured during the acoustic radiation tests and validation 
of the dynamic response prediction capability for this case will be discussed in the next section. 

6.6.4 Clamped Aluminum Panel 

Experiments were conducted to measure the structural-acoustic response of an aluminum 
panel subject to broadband, normal acoustic excitation. The experiments were performed in the 
Transmission Loss Apparatus (TLA) at NASA Langley Research Center, shown schematically 
in Figure 6.26. The results of the experiments were reported by Turner and Rizzi [55], but are 
briefly reiterated here. Vibration results are reported in this chapter and acoustic radiation results 
are reported in the Experimental Validation section of Chapter 7. 

The panel was installed in a picture-frame-type clamping fixture such that the unsupported 
dimensions were 0.305x0. 178xl.02e-3 m (12x7x0.040 inches). The entire assembly was 
mounted in the TLA partition baffle and subjected to normal-incidence, acoustic excitation 
provided by a speaker positioned approximately 0.914 m (36 inches) from the panel. The 
speaker was driven by a random signal, band-limited to a frequency range of 100-500 Hz. A 
centered 10x5 inch grid of 1 1 by 6 locations (total of 66 locations) was defined on the panel 
surface for vibration response measurements. Power spectra of the panel normal velocity at these 
locations were acquired with the laser vibrometer and a spectrum analyzer. 

The incident acoustic pressure was measured at several locations in the plane of the panel, 
with the panel and clamping fixture removed from the baffle, under identical source conditions 
to that of the response measurements. These source measurements were collected to assess the 
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uniformity of the incident pressure across the panel and to provide input for the computational 
procedure. The overall sound pressure level in the plane of the panel at the opening center was 
approximately 109 dB (ref 20 //Pa, 2.9E-9 psi). 

A 12x7-element mesh of the 24 degree-of-freedom rectangular plate elements was used to 
model the structure. Torsional boundary springs of 15.82 N*m/rad (140 in*lbf/rad) were employed 
to match the fundamental measured frequency. The source pressure spectrum obtained in the 
plane of the panel at the center of the baffle hole, when the panel and fixture were removed, was 
doubled and used as the blocked pressure input to determine the predicted panel response and 
acoustic radiation. The predicted results were generated using the following material properties 
for aluminum: E=10 Msi, //=0.33, and />=2.5756e— 4 lbf s 2 /in 4 . 

The power spectral density of the pressure incident upon the panel is shown in Figure 6.27. 
Comparisons of the predicted and measured panel normal velocity are shown in Figures 6.28 
and 6.29 for two measurement locations; near the panel center and in the lower-right quadrant. 
The measurement location coordinates indicated in the figure are relative to the panel’s lower 
left comer. The agreement is excellent with the exception of the appearance of additional peaks 
in the measured spectra. These peaks are attributable to the first and second non-symmetric 
modes of the panel. A number of factors could cause this type of phenomenon. The most likely 
culprit is the acoustic loading, but the source measurements indicated that the acoustic load was 
very nearly uniform and normal. Furthermore, it is expected that antisymmetric modes would 
be excited very weakly by acoustic pressure loading anomalies as previously discussed in the 
simply-supported aluminum panel section. It was later found that the panel was prestressed 
during fixturing, and it is believed that this is the reason for the significant antisymmetric modal 
response. 
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Figure 6.28 Comparison of predicted and measured panel normal velocity PSD at 
coordinates of x=0.152, y=0. 102 m (x=6, y=4 inches) relative to panel lower left corner. 
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Figure 6.29 Comparison of predicted and measured panel normal velocity PSD at 
coordinates of x=0.229, y=0.051 m (x=9, y=2 inches) relative to panel lower left corner. 


153 


The dynamic load models and spectral analysis method used in this study was presented in 
this chapter. The formulation and procedures were validated against known solutions prior 
to attempting correlation with experiment. The analysis was incrementally validated with 
experimental measurements for the cases of inertially loaded beams including an aluminum 
beam, a glass-epoxy beam and a SMAHC beam. Experimental results from the SMAHC beam 
specimen show enormous control authority over the resonance frequencies and RMS response 
level. Excellent agreement was achieved between prediction and experimental measurement, and 
reasons for existing discrepancies were discussed. The chapter is concluded with comparisons 
of the measured and predicted responses of a fully-clamped aluminum panel in preparation for 
acoustic radiation analysis and validation, which will be presented in the next chapter. 
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Chapter 7 
Acoustic Radiation 


The equations derived in Chapter 4 governing the static and dynamic response of panel- 
type SMA hybrid composite structures subjected to thermal and mechanical loads are applicable 
to a variety of geometric configurations and arbitrary natural boundary conditions. However, 
acoustic radiation and transmission loss analyses are only practical for panel-type structures 
with constraints around the entire perimeter of the structure. Furthermore, it is assumed that the 
structure is supported in such a way that flanking paths for the propagation of acoustic energy are 
prohibited, i.e., acoustic energy can only be transmitted by the flexible structure. Thus, consider 
a rectangular composite panel supported in an infinite rigid baffle (infinitely stiff/massive wall) 
separating two semi-infinite acoustic media. This configuration is typically simulated in practice 
by a transmission loss facility, shown schematically in Figure 6.26. Acoustic radiation analysis 
will be limited to the cases of acoustic pressure or concentrated force excitation on the source 
side of the panel. The panel radiates acoustic pressures to the receiving acoustic field due to its 
dynamic response. The equations governing acoustic radiation responses under these conditions 
will be developed in this chapter using continuous and discrete frequency domain approaches. 

Consider a panel supported in an infinite rigid baffle and radiating to a semi-infinite acoustic 
field, shown schematically in Figure 7.1. The acoustic pressure radiated by the panel to a point 
(x.y.z) in the acoustic field can be related to the dynamic response of the panel through the 
use of Rayleigh’s integral [56]: 
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(7.1) 


where p is the fluid density, w is the out-of-plane acceleration response of the panel, £ and ?/ are 
spatial coordinates in the plane of the panel, r is the distance from (£. //) on the plate to (x, y.z) 
in the acoustic field, c is the acoustic wave speed, l - r - is the retarded time between the two 
points, and A is the surface area of the panel. The subscript / indicates transmitted pressure to 
be consistent with the nomenclature in Chapter 6. 
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Figure 7.1 Schematic of a panel radiating to a semi-infinite acoustic field. 


7.1 Continuous Spectral Analysis 

The cross-correlation of the radiated pressure can be formed from the following expected 
value 
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(7.2) 


where the expected value expression in the integrand is the cross-correlation function of the 
plate normal acceleration i].t\ — >/ . /_> — ^), ?‘i is the distance from (£,?/) on the 

structure to ,r = (x,y,z) in the acoustic field, and r > is the distance from (£'. ?/) on the structure 
to x 1 = (.)■'. y'.z 1 ) in the acoustic field. For stationary response, the time lag can be defined as 
t = l > — 1\ so that Ajj jjjr = and the assumption of mean-square differentiability allows 
equation (7.2) to be rewritten in the form 


{■i-.!).z,x'.y'. z', r') = j j 2— (£, //. f'. //'. r) dA dA' 


(7.3) 
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where r' = r+ r - - ~— . Because the plate normal deflection cross-correlation function is related 
to the cross spectral density function through the Wiener-Khinchine relation 
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and the response is stationary, the following expression can be derived 
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Substituting in equation (7.3) results in the following expression for the cross-correlation of the 
radiated pressure. 
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Recall that the displacement field of the panel is known on the elemental level due to the 
finite element spatial discretization, shown in equation (4.5): 


tr, 


n)} M 


(7.7) 


where {/,„,(£. //)} is the vector of out-of-plane interpolation functions and {<■'/,} is the out-of- 
plane displacement vector, both for element e. Then, the displacement cross-correlation function 
between portions of the plate belonging to elements i and j can be approximated by 

('ww(t-v) ~ {Lw(Z'*1)}Je WhWlJ 


(7.8) 


where ] = K «« } is the physical displacement correlation matrix between elements 

/ and j. Consequently, the panel area integrals in equation (7.6) can be replaced by a double 
sum over element area integrals to obtain 
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(7.9) 

where Ay is the number of elements. It follows that the mean square (r = 0) radiated pressure 
at the point (x.y.z) in the acoustic field has the form 
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where r 

r, = yj (. 1 - - £ )■’ 4- (y - //)' + ~ 2 

and , 

r\ = + (y - v'f + =- (7 ' 12) 

are distances from the same point in the acoustic field to different points on the panel and 
( ri _ r \ )/c is the zero-time-lag retarded time difference (the difference in the time it takes for 
an acoustic disturbance to travel from (£.?/) to {x,y, z) versus (£',//) to (-f, //,-)). 
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Note that the spectral density matrix of the physical displacement response in the integrand 
of the spectral integral in equation (7.10): 



(7.13) 


is a subset of the system physical displacement spectral density matrix. The form for the system 
matrix is inferred from equation (6.38) 

[•‘vi.-tU’)] = [o] [N//M] [f/(u;)J [<?] r (7.14) 


where [o] is the modal matrix, [//(w)J is a diagonal matrix of stationary frequency response 
functions, and the modal force spectral density matrix is given by the appropriate form in 
equations (6.34). So, an arbitrary element {in.n) of the spectral integrand in equation (7.10) 
requires evaluation of an improper integral of the form 
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where the frequency response functions H m {u;) are given by equation (6.29). Note that these 
integrals have the same form as those required for determining modal correlation functions of 
the structure’s acceleration response for nonzero time lag (r / 0), shown in equation (6.36). 
Refer to Appendix F for sample cases of the evaluation of these integrals. 

The element area integrals in equation (7.10) can be evaluated via Gauss-Legendre numerical 
integration, which precludes the necessity for simplifications to the integrand. However, in the 
discussions to follow, it will be useful to introduce acoustic far-field assumptions. See Kinsler 
et al. [57] for a general description of these assumptions and Appendix D for details of this 
particular implementation. These assumptions entail making the approximations 
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where R is the distance from the center of the panel to the point (x.y.z) in the acoustic field, 
K = sin 0 cos p, and A - sin Osin o (see Figure 7.1). 

Equation (7.10) can be rewritten for the acoustic far-field as 
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The out-of-plane interpolation functions L w (£,y) are defined in terms of element local coordi- 
nates in equation (4.4). Therefore, the complex exponential function in the integrand of equation 
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(7.18) must be separated into two components; one is constant for each element index pair (/, j), 
the other must be included in the integration over the element areas in local coordinates. The 
element area integrations can be performed numerically, via Gauss-Legendre numerical integra- 
tion (see Appendix B), but significant improvements in computational efficiency can be realized 
by performing the integrations in closed form (see Appendix C). Other acoustic quantities, such 
as acoustic intensity and total radiated power, can be related to the mean square pressure pj in 
equation (7.18), shown in Appendix E. 


7.2 Discrete Spectral Analysis 

Recall the relationship between the radiated acoustic pressure at a point in the 

acoustic free-held and the dynamic response of the panel as described by Rayleigh s integral 
in equation (7.1). The physical degrees of freedom are known on an elemental basis from 
equation (4.5) 

uy (Z.vJ)*{L u .(t.n)} T {a b } (7.19) 

where {/.„,(£, ;/)} is a vector of out-of-plane interpolation functions and {07,} is the vector of 
out-of-plane degrees of freedom. So, the area integral in the expression for the radiated acoustic 
pressure, equation (7.1), can be rewritten as the sum of element area integrals. 
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where N c is the number of elements and .4, is the element area. 

A discrete spectral analysis can be employed by writing the panel dynamic response and 
radiated acoustic pressure as a sum of harmonic functions given by 
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where is the number of discrete frequencies and {07,,, } and [>,, are complex spectral amplitudes 
of the out-of-plane nodal displacement response and radiated acoustic pressure, respectively. 
Substitution of equations (7.21) and (7.22) in Rayleigh’s integral, equation (7.20), results in 
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This expression can be written for any one of the discrete frequencies as 
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If the excitation is deterministic, the expression for the radiated pressure complex spectral 
amplitudes can be used in an inverse FFT, equation (7.22), to determine the radiated pressure 
time response: 
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The vector { a } A is an element-k subset of the physical displacement spectral amplitude vector 
given by equations (6.41) and (6.59): 

{A bll } = [o}\H(^ n )\{f n } (7.26) 


However, for stochastic excitation, a statistical approach must be used. 

The cross-correlation function of the radiated acoustic pressure can be formulated from 
equation (7.22) through the following expectation 
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The cross-correlation function can be expressed in terms of the panel dynamic response by 
substituting from equation (7.24) as 
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where n is the distance from (£.;/) on the panel to .r = (x.y.z) in the acoustic field and r_> 
is the distance from (£',//) on the panel to = (.r'.y'.z 1 ) in the acoustic field. The term 
£[{«£ m }a {^6»i } / 1 conta i ns the element A- and / components of the generalized spectral matrix 
of the physical coordinate displacement response. 

Recall that the expected value of the spectral amplitudes takes a special form for stationary 
response, so the generalized spectral matrix simplifies to 
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where 6 m „ is the Kronecker delta. Then, the stationary radiated pressure cross-correlation 
function results ( ’ w >(.?i. i\. ,r_>. /_>) — » C vp (x \ . r ), which can be written as 
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where t — i > — t\ is the time lag. This expression can be further simplified for zero time lag 
(r = 0), which should not be confused with zero retarded time lag. Thus, the mean square 
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radiated acoustic pressure at a point in the acoustic field becomes 
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(7.31) 
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where rj and r\ are distances from the same point (r. y , ) in the acoustic field to different points 
on the panel, (/•] — r\)/c is the retarded time difference, and )] contains the element 

k and / entries of the single-sided, stationary, physical displacement spectral matrix [G.-u (“-'»)] 
inferred from equation (6.59): 

)] = [0] \ 11* U'v )J )] \ H(*'v )J [&] 1 (7.32) 


Although the following discussion is analogous to that for acoustic radiation using the 
continuous spectral approach, a brief description will be reiterated here for convenience and 
because the resulting relations are slightly different for the discrete spectral analysis. The element 
area integrals can be evaluated via Gauss-Legendre numerical integration, which precludes the 
necessity for acoustic far-field approximations. However, in the discussions to follow, it will 
be useful to introduce acoustic far-field assumptions for acoustic intensity calculations. See 
Appendix D for a detailed discussion. These assumptions entail making the approximations 
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where R is the distance from the panel center to the point (x.y.z) in the acoustic field 
K = sin 0 cos 4>, and A = sin# sin d, shown in Figure 7.1. 

The mean square radiated acoustic pressure can be written for the acoustic far-field as 
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Note that the interpolation functions £„.(£, ?/) are usually known in terms of element local 
coordinates. Therefore, the complex exponential function in the integrand of equation (7.35) 
must be separated into two components: one is constant for each element index pair (k.l), the 
other must be included in the integration over the element areas in local coordinates. The element 
area integrations can be performed numerically, via Gauss-Legendre numerical integration (see 
Appendix B), but significant improvements in computational efficiency can be realized by 
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performing the integrations in closed form (see Appendix C). Other acoustic quantities, such 
as the acoustic intensity and total radiated power, can be related to the mean square pressure in 
equation (7.35), shown in Appendix E. 

In general, the expression in equation (7.35) must be used to calculate the far-held mean 
square acoustic pressure radiated by a panel. However, recall the simplifications presented, 
starting with equation (6.61), for the excitation class which includes acoustic pressure, single 
or multiple equivalent concentrated forces, and base acceleration. A general expression for the 
single-sided, physical displacement spectral matrix at discrete circular frequency is inferred 
from equation (6.62): 

[fV, u (w„)] = { (*c„ ) } { H A (u,’„ ) } 7 (7.36) 


Substituting into the expression for the far-held mean square pressure and changing the order 
of the summations results in 
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and )} contains the integrated shape functions with the complex exponential. It is obvious 
that the double element sum has been rendered a single element sum for excitations of this 
type. Enormous computational savings result from this simplification. Additional savings can 
be achieved by forming {^(u,,,)} from the closed form integrations in Appendix C. 

7.3 Analytical Validation 

In the previous sections, forms for the far-held mean square pressure radiated by a flexible 
panel-type structure supported in an inhnite baffle and driven by acoustic pressure or concentrated 
force excitation were developed from Rayleigh’s integral using continuous and discrete spectral 
analysis methods. The formulation will now be validated through comparison with fundamental 
cases and experiment in the following two sections. 

7.3.1 Simply-Supported Aluminum Panel 

Comparisons will be made with the classical formulation developed by Roussos [52] for the 
case of a simply-supported aluminum panel. Consider a simply-supported aluminum panel with 
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Table 7.1 Material properties for aluminum panel specimen. 


Young’s Modulus 
Shear Modulus 
Poisson’s Ratio 
Mass Density 
Coeff. of Thermal Exp. 


68.95 GPa 10 Msi 

26.55 GPa 3.85 Msi 

0.3 

2.767e+3 kg/m 3 2.588e-4 lbf s 2 /in 4 

22.3 //</°C 12.4 fu/°F 


dimensions of 0.3556x0.254x0.001 m (14x 10x 0.040 inches). The material properties for the 
specimen are shown in Table 7.1. 

As a first example, let the excitation be characterized by normally-incident, planar acoustic 
waves. The acoustic excitation had a spectral level of 1.52 Pa 2 /Hz (3.2e-8 (psi)“/Hz) over a 
bandwidth of 10-1000 Hz with a spectral resolution of 0.5 Hz. The panel was modeled by a 
finite element mesh of 14x10 elements. Thirty modes were used in both the classical and finite 
element models with a uniform modal critical damping ratio of 0.02. 

A comparison of the tangential versus normal radiated acoustic intensity predicted by the 
classical and finite element methods is shown in Figure 7.2. The plot is for a radial distance of 
3.05 m (120 inches), polar angles of 6»=[0,tt/2] with a tt/ 18 increment, and an azimuthal angle of 
o=0. The same plot, but only considering the single frequency component at 800 Hz, is shown 
in Figure 7.3. It can be seen that the broadband frequency content tends to make the acoustic 
radiation directivity more diffuse as expected. The classical and finite element analyses show 
excellent agreement in both cases. 

The transmission loss versus frequency for this case is shown in Figure 7.4. The usual 
behavior of a dip at the fundamental frequency separating the stiffness-controlled region from 
the resonance-controlled region can be easily seen. The upper end in frequency is approaching the 
mass-controlled region of the transmission loss. The classical and finite element analyses show 
excellent agreement. The transmission loss curves are indistinguishable up to approximately 
800 Hz. 

Now consider planar acoustic excitation with the following angles of incidence (0^60 , 
0 - l= 45°). The excitation has a spectral level of 1.52 Pa~/Hz (3.2e-8 (psi) /Hz) over a 10—1000 
Hz bandwidth with a spectral resolution of 0.5 Hz. The plot in Figure 7.5 presents the tangential 
versus normal component of the radiated acoustic intensity at a radial distance of 3.05 m (120 
inches), polar angles of 0 =[O,tt/ 2] with a tt/18 increment, and an azimuthal angle of <p=0. The 
same plot for the single frequency component at 800 Hz is shown in Figure 7.6. These results 
further demonstrate the good agreement between the classical and finite element analyses. The 
transmission loss versus frequency for this case is shown in Figure 7.7. The antisymmetric 
modal resonances pointed out in Figure 6.7 are seen as tiny dips in the transmission loss curve. 

The acoustic radiation analysis using the FE method for structural response has been shown 
to produce very accurate results in comparison with classical analysis. The analysis will now 
be correlated with experimental measurement. 
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Figure 7.2 Normal versus tangential radiated acoustic intensity R=3.05 m 
(120 inches), 0 -[ 0,2tt], o= 0 for normal acoustic incidence, 10-1000 Hz. 
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Figure 7.4 Transmission loss versus frequency for normal acoustic incidence, 10-1000 Hz. 
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Figure 7.6 Normal versus tangential radiated acoustic intensity R=3.05 m 
(120 inches), ()=[ 0,2 tt], ©=0 for 0 { =6O, ©,=45 acoustic incidence, 800 Hz. 
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Figure 7.7 Transmission loss versus frequency for 9 {— 60, 45 acoustic incidence, 10—1000 Hz. 
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Table 7.2: Receiving microphone array measurement 
locations, shown schematically in Figure 6.26. 


Mic Position 

Radial 

Distance (m) 

Angle (deg) 

Mic Position 

Radial 

Distance (m) 

Angle (deg) 

R1 

| 


R5 


-30 

R2 

B 


R6 


-30 

R3 


30 

R7 

0.9144 

-45 

R4 

1.2192 

30 





7.4 Experimental Validation 

The experimental case considered in this section consists of a clamped aluminum panel 
subjected to normally incident, random acoustic excitation. Comparisons will be made between 
predicted and measured radiated acoustic pressures in a transmission loss configuration. 

7.4.1 Clamped Aluminum Panel 

Experiments were conducted to measure the structural-acoustic response of an aluminum 
panel subject to broadband, normal acoustic excitation. The experiments were performed in the 
Transmission Loss Apparatus (TLA) at NASA Langley Research Center, shown schematically in 
Figure 6.26. The experimental apparatus is described in some detail in the validation section of 
Chapter 6 and in Turner and Rizzi [55]. Comparisons of the predicted and measured panel normal 
velocity are shown in the previous chapter. Acoustic radiation results will be presented here. 

The transmitted acoustic pressure was measured using a single vertical microphone array 
consisting of three 1.27x 1CT 2 m (half-inch) microphones; one at the height of the panel center, 
one 0.3048 m (12 inches) above the center mic, and one 0.3048 m (12 inches) below the center 
mic. Transmitted acoustic pressure autospectra were collected at 21 locations by positioning the 
three-microphone array in each of seven measurement locations, defined in Table 7.2 and shown 
schematically in Figure 6.26. 

The power spectral density of the pressure incident upon the panel is shown in Figure 6.27. 
Comparisons of the predicted and measured transmitted pressure are shown in Figures 7.8 and 7.9 
for two measurement locations; radially perpendicular to the panel center and at an oblique angle. 
The agreement is very good, particularly near the panel resonances. Theoretically, the pressure at 
the top and bottom microphones of the array should be the same, due to symmetry. However, the 
measured spectra showed an increase in the levels between the peaks with decreasing distance 
to the hard floor in the receiving room. Thus, some of the discrepancies are attributable to 
receiving room reflections. Also note that the non-symmetric modes, shown in the previous 
panel response data in Chapter 6, do not radiate, as expected. 


170 












Pressure PSD, (psi) 2 /Hz 


Experiment 
Finite Element 




,V| vAv 


i i i i 


Frequency, Hz 


Figure 7.8 Comparison of the predicted and measured transmitted pressure PSD 
at the center microphone of array location R1 (see Table 7.2 and Figure 6.26). 
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Figure 7.9 Comparison of the predicted and measured transmitted pressure PSD at 
bottom microphone of array location R3 (see Table 7.2 and Figure 6.26). 




Chapter 8 

Applications and Demonstrations 


The tools developed and validated in the previous chapters will now be employed to 
demonstrate the static and dynamic thermomechanical behavior of SMAHC structures and 
compare performance gains achieved by embedded SMA actuators to conventional response 
abatement treatments. Results will be presented for SMAHC beam and panel structures. Early 
work performed in this study [16] showed that changes in SMA volume fraction have more 
effect on structural control authority than prestrain level. More recent work has corroborated this 
statement and shown that the SMA volume fraction is more easily varied/controlled in fabrication 
than the prestrain level. Finally, it is known that the thermomechanical fatigue performance of 
Nitinol actuators degrades substantially at prestrains greater than 5%. The results shown in this 
chapter will focus on conservative prestrain levels of 4% for the beam specimens and 5% for the 
panel specimens. For simplicity, a glass-epoxy matrix material system with constant material 
properties will be used in all cases. Mechanical property data for Nitinol from Cross et al. 
[28] will also be used. The material properties shown in Table 8.1 were used to generate the 
results shown in this chapter. Note that these predictions, using material properties for Nitinol 
from Cross et al. [28], are conservative relative to the properties measured from the Nitinol 
material used previously in this study, i.e., stable recovery stress for Nitinol used in this study 
at 4% prestrain is approximately 413 MPa (59.9 ksi) as compared to 310 MPa (45 ksi) for the 
corresponding case in Cross et al. [28]. 

Recall Figure 2.3 characterizing the recovery stress for Nitinol from Cross et al. [28]. The 
reverse transformation (austenitic) characteristic temperatures, under free recovery conditions, for 
the material represented in this plot are .4, «35°C (95°F) and A f *50°C (120°F). Note from 
Figure 2.3 that A f is increased by mechanical constraint and increasing prestrain, while .4, is not 
significantly affected by either. Furthermore, the “knee” in the recovery stress curve becomes less 
distinct with increasing prestrain. These are well-known phenomena and are attributable to stress 
increasingly inhibiting the formation of austenite and completion of the reverse transformation. 
These insights will be useful in later discussions. 

8.1 General Observations 

The physical behavior of SMA actuators embedded in a composite matrix is somewhat 
non-intuitive and is easily misinterpreted, so some general comments are deemed appropriate 
here. References will be made to Figures 8. la-8. Id, which show the non-deformed/reference 
(gray-shade) and deformed (mesh) in-plane geometries of 0.457x0.025x0.001 m (18x1x0.040 
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Table 8.1 Material properties for glass-epoxy and Nitinol. 


Glass/Epoxy 

Nitinol 

E, 53.78 GPa (7.8 Msi) 

A s 

32.22° C (90° F) 

E 2 17.93 GPa (2.6 Msi) 

Af 

62.78°C (145°F) 

G 12 8.62 GPa (1.25 Msi) 

17 T 

from Figure 2.3 

0.25 

E 

from Figure 2.4 

P 203 1 .8 kg/m 3 

V 

0.3 

(0.1 9x1 O' 3 lb s 2 /in 4 ) 



f'l 5.4x 10' 6 /°C (3.5 x 10" 6 /°F) 

P 

6450 kg/m 3 (0.6067 x 10 3 lbf s 2 /in. 4 ) 

30.6 x 10‘ 6 /°C (1 1.4x 10* 6 /°F) 

° la 

6.61 x 10 6 PC (3.67 x 10‘ 6 /°F), T<A S 
from equation 2.4, T>A S 


(, 2a 

6.61 x 10 6 PC (3.67 x 10 6 /°F), T<A S 
1 l.Ox HT 6 PC (6.1 lx 10 6 /°F), T>A f 
interpolation, A s <T<Af 


in.) conventional and SMAHC beams with a lamination of (45/0/-45/90) s and cantilevered or 
clamped-clamped boundary conditions. Consider the case of a cantilevered, conventional beam 
(no SMA). As the beam is heated, with uniform temperature, it will expand everywhere except 
at the boundary and will be stress-free at locations away from the boundary, shown in Figure 
8.1a. The thermal force that is generated is thereby relieved by in-plane deformation. If the other 
end is constrained, the thermal force in the axial direction cannot be relieved and a compressive 
in-plane load is generated. This load will increase with temperature until a critical value is 
reached, at which point the beam will buckle. Note that the thermal load in the transverse 
direction is balanced by reaction forces at the boundary and relieved by in-plane deformation in 
the transverse direction at locations away from the boundary (see Figure 8.1b). These physical 
phenomena can also be interpreted from the expression for the resultant internal force given by 
equation (2.38). As the cantilevered beam is heated, the thermal force r is balanced by the 
force due to in-plane deformation N mjr such that the resultant force N r is zero. Conversely, a 
nonzero resultant force exists for the clamped-clamped case because of an imbalance between 
and N mT . The imbalance is maintained after buckling by contributions from . 

An analogous, but inverse, situation exists for a SMAHC consisting of glass-epoxy with 
uniformly distributed SMA in all layers (volume fraction of 30% and initial strain of 4%). 
The thermal force in this case has a sense opposite of that expected because the effective 
CTE is contractile rather than expansive. A cantilever beam of this type would contract 
everywhere except at the boundary and will be stress free at locations away from the boundary, 
exemplified in Figure 8. lc. The axial thermal force is again relieved by in-plane deformation (i.e., 
A\y 7 > = \ „,j ). Imposing a constraint at the other end prevents the axial in-plane deformation 
and a substantial tensile in-plane force results. Note that the thermal forces (tensile) in the 
transverse direction are again balanced by reactions at the boundary and transverse in-plane 
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y+1000*v, 




x+1000*u, m 

Figure 8.1 Non-deformed/reference (gray-shade) and deformed (mesh) in-plane geometries for 
0.457x0.025x0.001 m (18x1x0.040 in.) cantilevered and clamped-clamped glass-epoxy and 
SMAHC beams. Displacements calculated at 65.56 °C (150 °F) and magnified 1000 times. 
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deformation at locations away from the boundary, shown in Figure 8. Id. Development of the 
tensile thermal force in the presence of mechanical constraints is the desirable attribute that is 
being explored for adaptive stiffening of composite structures by embedded SMA actuators. 

It is noted that in the previous discussions the arguments regarding SMAHC beams with 
free boundaries are, strictly speaking, only qualitatively correct because the SMA material 
properties used. Figure 2.3, result from constrained recovery measurements. However, the in- 
plane displacements associated with the tensile recovery forces are small relative to the prestrain 
imparted to the SMA actuators, so the predictions should be reasonably accurate for such cases. 
Also note that a stiffening effect does result in this particular case of a SMA hybrid cantilever 
because the modulus of the SMA is superior to that of the matrix material at temperatures where 
the SMA is in the austenitic state. Additional insights will now be gleaned from a series of 
examples. 


8.2 Thermal Buckling and Post-buckling 

Consider a glass-epoxy beam of 0.457x0.025x0.001 m (18x1x0.040 in.) dimensions with a 
lamination of (45/0/-45/90) s and clamped boundaries at both ends. The critical temperature 
change for this structure with a uniform temperature distribution is 1.68°C (3.03°F). The 
normalized maximum post-buckling deflection versus temperature change from ambient for the 
conventional beam (v a =0%) is shown in Figure 8.2 in comparison with SMAHC beams. The 
SMAHC beam responses correspond to the case of a 30% uniformly distributed volume fraction 
of SMA with an initial strain of 4%. The symbols v a and .l s indicate SMA volume fraction, 
prestrain level, and austenite start temperature, respectively. 

It can be seen that two buckling temperatures exist for the SMAHC beam with A .,=32.22° C 
(90°F); the first occurs at a temperature very near that for the conventional beam and the second at 
a much higher temperature. This phenomenon is attributable to the fact that .4, s , for this particular 
alloy, is greater than the ambient temperature of 21.1 1°C (70°F) and can be explained as follows. 
At temperatures less than . l s , the SMA is not activated and behaves like a conventional linear 
thermoelastic material. So the SMAHC beam buckles at a critical temperature below ,4.,. The 
SMA becomes activated at ,4 S and the recovery stress rapidly increases with temperature, which 
overcomes the compressive thermal stress and returns the beam to a flat configuration. At higher 
temperatures, the rate of increase of the recovery stress with temperature diminishes and may 
be surpassed by the rate of increase in matrix compressive thermal stress. Further increases 
in temperature will cause net-compressive stress development, eventually leading to a second 
incidence of thermal buckling. 

In contrast, the ,4.,=15.56°C (60° F) case in Figure 8.2 exhibits only a single thermal buckling 
event at high temperature. This case was produced using the same SMA properties as before by 
shifting the temperature scale associated with those properties -16.67°C (-30°F). This can easily 
be achieved in practice by increasing the Nickel content of the alloy as little as 0.125 weight % 
(see Buehler and Wang [58]). Note that the buckling temperature is slightly higher in this case 
relative to the .4.,=32.22°C (90°F) case. This result is due to the fact that the shift in the SMA 
property temperature scale for the ,4. S =15.56°C (60°F) case produces a slightly higher recovery 
force at the imposed temperature, which can be seen by referring to Figure 2.3. 
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Figure 8.2 Normalized mid-span post-buckling deflection versus temperature change for 
conventional and SMAHC 0.457x0.025x0.001 m (18x1x0.040 in.) clamped beams. 
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Recall that the former effect, near-ambient buckling of a structure then complete arrestment of 
the thermal post-buckling deflection by SMA recovery stress, was demonstrated experimentally 
and correlated with prediction in Chapter 5 for the case of a clamped-clamped SMAHC beam. 
Proof of this behavior has a numerical basis also, which will be explained subsequently. There 
have been differing theories about the elimination of the thermal post-buckling deflection, which 
results in two buckling temperatures for A s greater than ambient. The formulation by Turner, 
Zhong, and Mei [16], based upon the RSM of the SMA constitutive behavior, corroborates the 
results presented here. A different result was reported in Mei, Zhong, and Turner [59] where 
an incremental approach with many temperature increments was deemed necessary to account 
for the SMA material nonlinearity in arriving at the final temperature. That formulation is also 
based upon the RSM, but does not exhibit the complete elimination of thermal post-buckling 
deflection for cases with l s greater than ambient. The difference between the two methods 
was attributed to a mismatch of recovery stress and out-of-plane deflection terms in a thermal 
post-buckling solution with multiple temperature increments based upon the formulation given 
in Turner, Zhong, and Mei [16]. The ECTEM of the SMA constitutive behavior, presented 
in this work, shows that the use of small temperature increments is not necessary because the 
model inherently takes the material nonlinearity into account. Furthermore, any mismatch in 
terms associated with multiple temperature increments in the thermal post-buckling solution can 
be avoided by using a single temperature increment. The results of this approach are shown in 
Figure 8.3 in comparison to the corresponding multi-step solution from Figure 8.2. The number 
of temperature steps used in each post-buckling solution is indicated in the figure by .Y^y . Each 
discrete symbol represents a separate solution achieved by initiating a post-buckling solution at 
the first buckling temperature with a single temperature increment. These results show that the 
present solution is insensitive to the number of temperature increments and complete elimination 
of thermal post-buckling deflection for cases with A s greater than ambient and a flat, stress free 
initial configuration is a general theoretically valid result. The ability to use a single temperature 
increment can result in a huge computational savings also, particularly for high temperatures 
relative to the buckling temperature. 

Recall that the SMAHC beams considered thus far employed a uniformly distributed SMA 
volume fraction of 30%. It is intuitive that performance gains could be achieved by strategically 
placing the SMA actuators. One example of this technique is shown in Figure 8.3, where the 
thermal post-buckling deflection is shown for a SMAHC beam with the same (45/0/-45/90) s 
lamination but SMA actuators in the 0° layers only. A 50% volume fraction in the two 0° layers 
results in an overall volume fraction of 12.5%. This overall volume fraction is less than half 
that of the beam with uniformly distributed actuators, but the performance is slightly improved. 
Note that this approach is more practical from a fabrication point of view as well, particularly 
for automated fabrication where one could imagine replacing strategic spools of prepreg with 
SMA ribbon in a continuous consolidation fabrication operation. 

The effect of changing the SMA volume fraction within the 0° layers is shown in the thermal 
post-buckling results in Figure 8.4. It can be seen that a volume fraction of 2.5% is insufficient 
to eliminate thermal post-buckling, but its effect in reducing the deflection is clearly observed. 
Thermal post-buckling can be eliminated for a range of temperatures with as little as 7.5% SMA 
and the case of 12.5% SMA exhibits deflection elimination for a broad range of temperatures. 
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Figure 8.3 Normalized mid-span post-buckling deflection versus temperature change for 
SMAHC 0.457x0.025x0.001 m (18x1x0.040 in.) clamped beams with SMA in all/0° layers. 
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Figure 8.4 Normalized mid-span post-buckling deflection versus 
temperature change for composite 0.457x0.025x0.001 m (18x1x0.040 
in.) clamped beams with various SMA volume fractions. 
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Similar behavior is exhibited by panel-type specimens as demonstrated in Figure 8.5. This 
figure shows the normalized center deflection as a function of temperature change for a glass- 
epoxy panel with dimensions of 0.36x0.25x0.001 m (14x10x0.045 in.), (45/0/-45/90/0/90/- 
45/0/45) lamination, a varying SMA volume fraction, and clamped boundary conditions. As 
in the latter beam specimen case, the SMA in the panel specimens is assumed to reside in the 
90° and outer 0° layers only. Although it is not shown, it is noted that slight performance gam 
is achieved under these conditions relative to a corresponding laminate with SMA in all layers 
such that the overall SMA volume fraction is the same (see Turner [60]). This material system, 
i.e., with SMA in only a few layers, is consistent with the thought of automated fabrication made 
in the penultimate paragraph. Also note that the middle 0° -layer is necessary to separate the 
two 90° layers containing a substantial SMA volume fraction, thereby avoiding consolidation 
problems during cure. The results in Figure 8.5 show that a volume fraction of 4% is insufficient 
to eliminate thermal post-buckling in this case, but its effect in reducing the deflection is clearly 
observed. As little as 13% SMA will eliminate thermal post-buckling over a small range of 
temperatures; higher concentrations of SMA (e.g., 22%) can greatly extend the temperature 
range. The importance of these effects is more clear when considering the dynamic response 
of these structures. 


8.3 Dynamic Response 


Dynamic response and subsequent sonic fatigue and acoustic transmission analyses will focus 
on panel specimens with dimensions of 0.36x0.25x0.001 m (14x10x0.045 in.) and clamped 
boundary conditions. The dynamic response is assumed to be forced by a normally-incident, 
random acoustic pressure with an overall sound pressure level (OASPL) of 120 dB (ref 20// Pa) 
and a bandwidth of 10-400 Hz. A uniform modal critical damping ratio of 1% was used in 
the dynamic response calculations. The maximum displacement response power spectral density 
(PSD) of the SMAHC panel is shown in comparison to that of the baseline panel in Figure 
8.6. It is clear that the embedded SMA has little effect at ambient temperature. At 82°C 
(180°F,A7 =61°C), however, substantial benefit from the SMA recovery stress is evident in 
shifting the fundamental frequency, shifting other modes out of bandwidth, and reducing peak 
responses. Thus, performance can be enhanced by shifting resonance out of the excitation 
bandwidth and by response reduction within the excitation bandwidth. 


The thermal post-buckling observations have strong implications for the panels’ dynamic 
response. An example is shown in Figure 8.7, where the maximum displacement PSD of 
the SMAHC panel (r„=22%) is shown at ambient and two elevated temperatures. At 82° C 
(180°F, AT=61°C), substantial benefit of the embedded SMA is evident, while further heating 
to 121°C (250°F, A7’=100°C) greatly diminishes the performance (RMS W max =0.0927e-3 m) 
because continued thermal expansion has partially overcome the beneficial effects of the recovery 
stress. This particular material system, glass-epoxy matrix with 22% Nitinol volume fraction at 
5% prestrain, achieves its “optimal” dynamic response at approximately 82 C (1 80 F), which 
corresponds to the “knee” in the recovery stress curve (see Figure 2.3). This result also holds 
for analogous SMAHC panels (r„=22%) with embedded Nitinol having 3% and 4% prestrain, 
which exhibit “optimal” dynamic performance at their respective “knee temperatures of 54 C 
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Figure 8.5 Post-buckling deflection of clamped 0.36x0.25x0.001 
m (14x10x0.045 in.) SMAHC laminates versus AT and v„. 
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Figure 8.6 VI muJ PSD for 0.36x0.25x0.001 m (14x10x0.045 in.) 
baseline and SMAHC laminates with clamped boundary conditions. 
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( 130°F) and 66°C ( 150°F). The optimal temperature for dynamic performance will be designated 
by A„ and is governed by the same thermoelastic balance as that described in the thermal post- 
buckling discussion. In general, with a different material system or at higher prestrain levels (in 
the absence of a distinct “knee”), A„ must be determined numerically or by comparing the slope 
of the recovery stress curve with that of the matrix thermal stress versus temperature. 

Selection of an alloy composition and processing combination that will place A„ near the 
application temperature is obviously recommended for dynamic control applications. An alternate 
approach to achieve similar performance would entail increasing the SMA volume fraction and/or 
prestrain, which will result in a less efficient and/or durable design. Requirements previously 
cited for desirable static and dynamic performance imply potential for a delicate balance between 
avoiding thermal buckling by imposing .4., to be less than ambient temperature and optimizing 
dynamic performance by specifying A 0 to be near the application temperature. Thermal buckling 
can also be avoided by boundary expansion, a common approach for aerospace structures, so .4 , 
can be neglected in such cases and A„ should be used to guide efficient design. 


8.4 Sonic Fatigue and Acoustic Transmission 

The performance of the SMAHC laminate will now be compared with conventional ap- 
proaches including increased stiffness through additional composite layers and a constrained 
layer damping (CLD) treatment. The conventional approaches were analyzed at ambient tem- 
perature to simulate the conditions of boundary expansion, i.e., no thermal buckling. No attempt 
was made to account for changes in matrix properties with temperature. A variety of dynamic 
response performance characteristics are compiled in Table 8.2 for convenience of comparison. 
The SMAHC panel achieves a 77% RMS displacement response reduction, relative to the base- 
line 9-layer glass-epoxy panel, with a 48% mass increase. A conventional composite laminate of 
13 layers (45/0/-45/90/45/0/90/0/45/90/-45/0/45) exhibits a similar mass increase (44%), but only 
achieves a 61% RMS displacement reduction. A RMS displacement response of 0.0474e-3 m 
(nearly equivalent to that of the SMAHC panel) can be achieved with a 16-layer (45/0/-45/90)2 S 
conventional laminate, but at the expense of a 78% mass increase. A typical CLD treatment 
was modeled as a 19.16 N/m 2 (0.2 lbf/ft 2 ) weight addition and a uniform modal critical damp- 
ing ratio of 6%, which are considered to be representative values. This treatment also happens 
to result in a mass increase (42%) similar to that of the SMAHC, but the RMS displacement 
response reduction of 63% is again inferior to the SMAHC. Additional insights are more clear 
from comparisons of the displacement response spectra, as discussed subsequently. 

A comparison of the respective maximum displacement response PSDs is shown in Figure 
8.8. As anticipated, the means by which the three systems achieve dynamic response reduction 
is quite different. The SMAHC achieves a much more weight-efficient stiffening effect relative 
to that of adding layers, while the CLD treatment achieves reduction primarily by damping 
enhancement. The weight-efficient stiffening effect of the embedded SMA is advantageous 
for combination with other passive treatments, which are typically more effective at higher 
frequencies. SMAHCs are also attractive for combination with control approaches involving 
piezoelectric actuators, which typically exhibit actuation authority roll-off below 200 Hz. Such 
combined active-passive approaches could achieve enormous dynamic response reductions by 
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Table 8.2: Results from clamped 0.36x0.25x0.001 m (14x10x0.045 in.) baseline, 
SMAHC, and panels with conventional response abatement treatments. 


Case 

Added 

Mass 

fi.Hz 

. 

RMS W max , 
m 

TL, dB 

Max RMS fit 
@ 120dB 


Baseline 

0% 

87.48 

0.2062e-3 

17.9 

42.0 


SMAHC @ 82°C 



0.0476e-3 

14.3 


205e+6 

13-Layer 


mm 

0.0805e-3 

19.1 



16-Layer 

78% 

158.6 

0.0474e-3 

20.2 

16.8 

23.4e+6 

Baseline-CLD 

42% 

73.4 

0.0768e-3 

27.4 

15.6 

37.8e+6 


using the embedded SMA to shift resonance frequencies and provide peak response control by 
stiffness enhancement, while the active approach would provide enhanced damping at resonances 
still within the excitation bandwidth. 

The dynamic response results have direct implication for sonic fatigue life prediction. 
Consider an analogous acoustic pressure loading condition with an OASPL of 150 dB (ref 20//Pa). 
It is clear that linear dynamic response analysis is no longer valid, but its use at high response 
levels is still common practice in the aerospace industry for sonic fatigue design. This condition 
is more representative of a sonic fatigue environment and will produce more realistic life 
predictions, while it is understood that conservative estimates will be obtained. A representative 
strain-life (S-N) relationship for glass-epoxy laminates was modeled by < = 7385f -6 (A’P°' ] >51 
and the maximum RMS normal strain was taken as the failure criterion. In all cases, the 
maximum strain was the transverse normal strain at the mid-span of the long boundary. Under 
these conditions, the baseline panel has an expected life of 63.7k cycles. The corresponding 
life predictions for the SMAHC, 13-layer conventional, and CLD panels are 205M, 2.84M, 
and 37. 8M cycles, respectively. These estimates are in direct relation with the corresponding 
peak displacements in Figure 8.8. It is intuitive that measures to further attenuate the SMAHC 
peak response would result in excellent sonic fatigue performance. The previously mentioned 
active-passive system, for example, might be a particularly attractive solution. 

A plot of the transmission loss (TL) versus frequency for the SMAHC panel at two 
temperatures is shown in comparison to the baseline panel in Figure 8.9. The main effect 
of the embedded SMA in its most effective state at 82°C (180°F) is a vast increase in the 
stiffness-controlled portion of the TL. This improvement is diminished at higher temperatures 
(e.g., 121°C) for the same reasons as stated above for the dynamic response. Although the 
stiffness-controlled TL for the SMAHC panel shows dramatic improvement, the TL over the 
400 Hz bandwidth is actually diminished because the mass-law effect has been moved to higher 
frequencies. The TL for the SMAHC panel at 82°C (180°F) is shown in comparison to that for 
the two conventional approaches in Figure 8.10. It can be seen that the SMAHC panel is superior 
over a 200 Hz bandwidth, but not as effective at mid-range (resonance -con trolled) frequencies 
due to the shift in mass-law effects. Combination of the SMAHC with other passive and/or 
active approaches is again an attractive alternative for this frequency range. At high frequencies 
(mass controlled) the three approaches will produce nearly equivalent TL. 
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Figure 8.8 W mar PSD comparison of two SMAHCs with conventional treatments. 
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Figure 8.9 TL versus frequency for the baseline and SMAHC laminates. 
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Figure 8.10 TL comparison of two SMAHCs with conventional treatments. 
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Chapter 9 

Summary and Concluding Remarks 


The main objectives of this work were to develop a themomechanical model for predicting 
the behavior of composite structures with embedded shape memory alloy (SMA) actuators and 
to validate that model with experimental measurements. Other objectives were to develop SMA 
actuator processing/characterization and SMA hybrid composite (SMAHC) fabrication techniques 
for manufacturing structures with predictable and repeatable behavior. Finally, the experimental 
results from the SMAHC specimens developed in this work were intended to demonstrate the 
control authority of SMA actuators and demonstrate some of the complex behavior of SMAHC 
structures. 

A fairly extensive review of SMAs and their thermomechanical behavior was given in 
Chapter 2 to introduce the topic and the demands on a constitutive model. A constitutive model 
was then developed by a mechanics-of-materials approach to describe the behavior of a composite 
lamina with an embedded SMA constituent. The thermoelastic model captures the material 
nonlinearity of the SMA thermal strain/stress with temperature and is valid for constrained, 
restrained, or free recovery behavior with appropriate measurements of fundamental engineering 
properties’. The constitutive model is also capable of modeling matrix material nonlinearity and 
is in a form that is amenable to implementation in a commercial analysis environment. Classical 
lamination theory was used to derive the resultant force and moment relations from the SMAHC 
constitutive relations and von Karman strain-displacement relations. 

Fabrication and material system characterization processes for SMAHC specimens were 
described in Chapter 3. The particular SMAHC material system used in this study consists of 
a glass-epoxy, unidirectional, prepreg, tape with embedded Nitinol ribbon. Several practical 
solutions for achieving repeatable results and overcoming practical fabrication issues were 
demonstrated. 

Extensive testing of the glass-epoxy and Nitinol materials was conducted to quantify the 
behavior of a structure fabricated from the SMAHC material system. The mechanical properties 
of the glass-epoxy in principal material directions (E|, E 2 and G i2 , r»i, n 2> were determined 
as a function of temperature from unidirectional and angle-ply tensile specimens. Isothermal 
(modulus) and isolength (recovery force) tests were performed on Nitinol ribbon samples in 
two conditions to determine the material properties as a function of temperature and thermal 
cycle. Thermomechanical properties from the two conditions showed similar trends, but some 
quantitative differences were evident. These differences were attributable mainly to aging effects 
and stressed the importance of testing the actuator material in a condition that closely matches the 
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intended application as the actuator properties are very sensitive to condition and are a dominant 
effect in the SMAHC performance. 

In Chapter 4, the resultant force and moment relations were used along with the nodal degrees 
of freedom and interpolation functions associated with a rectangular 24 degree-of-freedom plate 
element to form a potential energy functional. The extended Hamilton’s principle was used to 
derive the equations governing the response of a SMAHC plate-type structure to steady-state 
thermal and dynamic mechanical loads. The governing equations were then separated into static 
and dynamic components. 

In Chapter 5, the equations governing static stability were derived within the framework 
of the present formulation and used to develop a thermal buckling solution approach. The 
iterative method was shown to produce valid solutions for an arbitrary material system with 
nonlinear matrix material properties in addition to the nonlinearities associated with the SMA. 
A thermal post-buckling solution method was also developed to predict the thermally deflected 
shape subsequent to a thermal buckling event. An iterative Newton-Raphson approach was used 
to account for geometric nonlinearities associated with large thermal deflections. The thermal 
buckling and post-buckling solutions were validated against known solutions for conventional 
structures to benchmark the formulation and implementation. 

A SMAHC beam specimen was tested to determine its thermoelastic response in a clamped- 
clamped configuration with a thermal load generated by resistive heating. The out-of-plane 
deflection was measured as a function of temperature at the beam mid-span by a laser dis- 
placement transducer. The predicted thermoelastic response was in excellent agreement with 
the experimental measurement. Existing discrepancies were attributable to the effects of initial 
deflection and mechanical boundary compliance in the physical system. 

Continuous and discrete frequency domain methods were applied in Chapter 6 to the dynamic 
equations of motion to predict geometrically linear dynamic responses of SMAHC plate-type 
structures to various mechanical loads. The same approach was used in Chapter 7 with a 
Rayleigh integral method to develop expressions for the radiated acoustic pressure resulting 
from the dynamic motion of the structure. These solutions/implementations were validated 
against classical solutions to benchmark the analyses. Further validation was performed with 
experiment for conventional structures to assess the capability of the analyses to capture the 
behavior of a physical system. 

A SMAHC beam specimen was tested in a clamped-clamped configuration with a uniform 
thermal load generated by resistive heating and band-limited random inertial excitation. Dynamic 
response data was collected at ambient and at elevated temperatures from 32.22°C (90°F) 
to 121.11°C (250°F) in 5.56°C (10°F) increments. The fundamental frequency of the beam 
was shifted 430% at 121.1 1°C (250°F), relative to the ambient temperature resonance. The 
maximum RMS displacement response attenuation of 84% was achieved at a temperature of 
71.1 1°C (160°F) rather than at the maximum temperature. This was apparently caused by 
enhanced damping due to stress induced martensitic transformation hysteresis. Nonetheless, 
similar attenuations of around 80% were achieved at higher temperatures. 

The predicted dynamic responses were in good agreement with experimental measurement 
with the exception of a slightly over-predicted stiffening effect at elevated temperatures. This 
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discrepancy was again attributable to the fact that the mechanical boundaries were modeled as 
perfectly rigid, when in fact the presence of thermal/electrical insulation in the gripping apparatus 
resulted in significant boundary compliance. This effect was incorporated into the model, and 
excellent agreement between prediction and measurement was achieved at all temperatures. 

The model and analysis tools developed and validated in this work were then used in Chapter 
8 to demonstrate a variety of thermomechanical response phenomena associated with SMAHC 
structures and to compare the performance of a SMAHC panel-type specimen with conventional 
dynamic response abatement techniques. Two buckling temperatures can exist for SMAHC 
structures with a SMA constituent having an austenite start temperature A s greater than ambient. 
This effect results from the SMA recovery stress “pulling” the structure flat after initially buckling 
at a temperature less than A x and buckling again when continued thermal expansion overcomes 
the recovery stress. The low -temperature buckling phenomenon can be avoided by selecting 
a SMA with A x less than ambient, thereby suppressing buckling altogether for a large range 
of temperatures above ambient. This effect was later supported by experimental observations 
described in Chapters 5 and 6. 

Substantial improvements in performance can be realized by selectively placing the embed- 
ded SMA in strategic locations such that the SMA volume fraction can be used more effectively, 
which is also driven by fabrication considerations. Relatively modest changes in the volume 
fraction in such cases can result in enormous changes in the performance. SMAHC structures 
exhibit an “optimal” operating temperature (1„) for dynamic response performance. A„ is a 
function of the matrix material and SMA composition, volume fraction, and prestrain level. For 
many cases A 0 corresponds to the “knee” temperature in the recovery stress curve. In general, 
A a must be determined numerically or by comparing the slope of the recovery stress curve with 
that of the matrix thermal stress versus temperature. The implied conflict in choosing an alloy 
for aerospace structures with A„ < T„ to avoid thermal buckling and A 0 ~ 7 to “optimize 
dynamic response can be resolved by taking the conventional approach of allowing boundary 
expansion to prevent buckling and using A 0 ~ T to guide structural design. 

The performance of a SMAHC panel was compared to that of a baseline panel and two 
conventional approaches including additional conventional composite layers and a constrained 
layer damping (CLD) treatment. The three response abatement approaches had similar added 
mass relative to the baseline panel, but the SMAHC panel exhibited superior dynamic response 
and sonic fatigue behavior. The SMAHC panel also exhibited superior transmission loss (TL) 
characteristics in the low-frequency, stiffness-controlled region. The weight-efficient stiffening 
effect of the embedded SMA makes it ideal for low-frequency, noise transmission abatement 
and for combination with other passive or active approaches, which are typically more effective 
at higher frequencies. These combined approaches have significant implications for dynamic 
response/sonic fatigue also. 

The major conclusions from the preceding discussions can be summarized as follows: 

1) A practical constitutive model has been developed that accurately captures the mechanics of 
a SMAHC material system. 

2) Reliable SMAHC specimen fabrication procedures were developed and SMAHC beam 
specimens were fabricated. 

3) Experimental measurements from the SMAHC beam specimens demonstrated enormous 
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control over thermal buckling, post-buckling, and vibration responses. 

4) Predicted static and dynamic responses using the model developed in this study were in 
excellent agreement with the experimental measurements. 

5) A SMAHC structure can exhibit two buckling temperatures. This phenomenon was supported 
by experimental measurements. 

6) SMAHC structures can exhibit an optimal operating temperature for dynamic response 
performance. 

7) SMAHC panel structures exhibit superior dynamic response characteristics to the conventional 
approaches considered. 

8) SMAHC structures are ideal for low-frequency TL and for combination with other treatments 
for dynamic response and noise transmission abatement. 

In conclusion, this study takes an engineering approach at considering embedded SMA 
actuators as a means of achieving quasi-static control of composite structures through adaptive 
stiffening in response to environmental (temperature) effects. The model was validated against 
experimental measurement for static and dynamic responses. It is believed that implementation 
of this model in a commercial FE code will lead to a very useful and practical analysis/design 
tool. It is hoped that this study will provide useful information for the research and engineering 
community in regard to fabrication, characterization, and modeling of SMAHC structures. 
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Appendix A 

RSM Hybrid Lamina Constitutive Relations 

The constitutive relations developed in this appendix are referred to as the recovery stress 
model (RSM) as described in Chapter 2. A mechanics-of-materials approach is used to derive 
the thermoelastic constitutive relations governing a SMAHC lamina, based upon observation of 
the recovery stress phenomenon. 

A representative volume element of a SMA fiber-reinforced hybrid composite lamina is 
shown in Figure A.l. The element is taken to be in the plane of the plate. The composite 
matrix, for example graphite-epoxy, has the principal material directions 1 and 2, and the SMA 
fiber is embedded in the 1 -direction. 
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Figure A.l Representative volume element for SMA fiber-reinforced hybrid composite lamina. 


In order to derive the constitutive relation for the 1-direction, it is assumed that a stress a\ 
acts alone on the element (/r_> = 0) and that the SMA fiber and composite matrix are strained 
by the same amount, t \ (i.e., plane sections remain plane). The 1-direction stress-strain relation 
of the SMA fiber can be described as 


^l« = E „ ' 1 +<?,■ • T > A s (A.l) 

or 

aia - L\,(ti -o„AT) . T < A s (A.2) 

where A s is the austenite start temperature and o„ is the thermal expansion coefficient. The 
Young’s modulus E a and the recovery stress a x are temperature dependent. The recovery stress a x 
is also dependent on the initial strain ( r . Description of the recovery stress and modulus behavior 
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for Nitinol is discussed in Chapters 2 and 3. Similarly, the one-dimensional stress-strain relation 
in the 1 -direction for the composite matrix can be expressed as 

(A. 3) 


O I 


i A 7 ’ ) 


V]m — E 1 m ( ( 1 

The resultant force in the 1 -direction {rr 2 = 0) is distributed over the SMA fiber and 
composite matrix and can be written as 

(A. 4) 


(T i j\ i = (T 1 a + <J \ m A 


in 


where (rr,, . 1 , ), (<t 1w , .4„), and (<7„„.,4i,„) are the (stress, cross-sectional area) of the entire 
element, SMA fiber, and composite matrix, respectively. Thus, the average stress oj is 

o’ i = + (T|, (A. 5) 

where r„ = A u /A { and r m = A m /Ai are the volume fractions of SMA and composite matnx, 
respectively. When T > 4 S , the SMA effect is activated and the one-dimensional stress-strain 
relation in the 1 -direction becomes 


rr\ = ( K {1 ( 1 + (7 r ) I'a + F'\m(t\ O [ A / ) l' m 
= L\t \ — h [ rn (\ \ m I'm A I 

(A.6) 

where 

F ' 1 — -£■ 1 m * 1 in H" E a i \t 

(A.7) 

When r < 4 S , the SMA effect is not activated and the stress a\ is 


cr \ — h i ( i — (E a c\ tt v a 4" E\ m c\\ n ,Vtn)^F 
= E l (t l -a i AT) 

(A. 8) 

where „ v 

( riniO imCm i ri a r (l j ( A Q) 

( f - 1 ni 4" 

A similar constitutive relation may be derived for the 2— direction by assuming that the 
applied stress rr> acts upon both the fiber and the matrix (o ] =0). Thus, the one-dimensional 
stress-strain relations in the 2-direction for the SMA fiber and the composite matrix become 

O 2s = <72 = Ea{( >a ~ a a AT) 

(A.10) 

and 

rr ltu = <r> = E-iwit ’m ~ a 2m AT) 

(A. 1 1) 


respectively. The recovery stress does not appear in equation (A. 10), since the SMA fiber initial 
strain t r and recovery stress o T are considered to be a 1-direction effect only. 

The total elongation is due to strain in the composite matrix and the SMA fiber and may 
be written in the form 


.4 1 f j = A m (>m + A u (>a 


(A. 12) 
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Thus, the total strain becomes 


f 2 — ( 2in t 'm "I" ( 2a *\t 


(A. 13) 


Since cr> = — n_> A/'), equations (A. 10) and (A. 11) may be substituted into equation 

(A. 13) to give 
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(A.14) 


Therefore, the modulus and thermal expansion coefficient in the 2-direction become 


and 



(A.15) 

(Ejvt l'a “t" h(ii'nt) 

^ 2m I'm ^d*\i 

(A. 16) 


Expressions for the hybrid composite Poisson’s ratios and shear moduli follow from similar 
derivations. 

The constitutive relations for a thin composite lamina with embedded SMA fibers can be 
derived using a similar engineering approach to give 
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and 



(A. 17) 


(A. 18) 


where [(,)]„, and [Q*] are the reduced stiffness matrices of the composite matrix and the SMA 
hybrid composite lamina, respectively. The [(}*] matrix is temperature dependent and is evaluated 
using the previously derived relations as 


( R\ , //]_> ) — ( E \ m , l'\2m ) t'ni T ( E a . V„ )v a 


(A.19) 


and 


( E > • G i _) ) 


[(£'2 m 


(E 2m E a .Gi 2m G a ) 

' "b (E a , G a )<’/ji 


(A.20) 


where the v ' s are Poisson's ratios and the 6”s are the shear moduli. The thermal expansion 
coefficients oj and <\ > are derived in equations (A.9) and (A. 16). 
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Appendix B 

Gauss-Legendre Numerical Integration 


B.1 Element Matrices 

The element area integrals necessary for evaluation of the element matrices in equations 
(4.55) can be performed numerically using Gauss-Legendre quadrature. The four-node, rectan- 
gular, 24-dof plate element used in this study has interpolation functions given by equations 
(4.3) and (4.4) for the in-plane and out-of-plane degrees of freedom, respectively. Consider the 
element bending mass matrix for example, i.e.. 


["'ft] 


I pli{ //„.(.r. y)}{ //»'(•<'. u)} 1 <IA, 

A, 


b a 

= j j pk{H w (-r-u)}{HM-y)} T 
0 0 




(B.l) 


where (.v.y) are the element coordinates and a and b are the element length and width, 
respectively. 

The element area integrations in equation (B.l) are in a form that is easily adapted to 
Gauss-Legendre numerical integration. The coordinates can be transformed to non-dimensional 
element-centered coordinates (£.?/) with the variable transformation 


.r : [0, a] £ = 2- - 1 : [-1.1] <l£ = &r 

a ~ 

y : [0. b] y = ~ 1 : [- 1 - 1 ] dy = ^ dy 


(B.2) 


Thus, the expression in equation (B.l) can be rewritten as 


["'ft] 


T 

= f J ph{H w [Z,Ti)}{H u ,(Z,v)} \J(Z'V)\<Wv 

-i -1 


(B.3) 


where |./(C '/)| is the determinant of the Jacobian matrix of the transformation: 

m,v)\ 


Or 

dr 


d v 

dy 
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Z 0 
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o i 
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ab 


(BA) 


The element displacement functions are now written in terms of the non-dimensional coordinates 

(£- '/) as 


y) = X 


j=i 


C( t (0 (u(y)u'i 




(B.5) 
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Figure B.l Schematic of plate element in centroidal non-dimensional coordinates. 


where 


<V;(0 = j(-6£ 3 + + 2) 

",(o = + u 2 -t-o) 


(B.6) 


and (O- 'li) are the coordinates of node i. A schematic of an element in the non-dimensional 
coordinates is shown in Figure B.l 

The integral over the element are in equation (B.3) is now in a form that is easily evaluated 
numerically. Gauss-Legendre integration entails evaluating the integrand at a number of equally 
spaced (Gauss) points within the region of integration and summing the product of these 
contributions with corresponding weights. Exact evaluation is possible for polynomial-type 
functions when the order of the function is less than or equal to one less than twice the number 
of Gauss points, i.e., (2N g — 1 > p) where N y is the number of Gauss points and p is the order 
of the polynomial. Employing this method, the integral in equation (B.3) can be rewritten as 


K] = \J(ii-Vj)\ (B.7) 


»=1 J =1 


where 11* are the weights at the Gauss points denoted by and i/ r 


B.2 Mean Square Radiated Acoustic Pressure 

A similar approach can be used to evaluate the element area integrals for acoustic radiation 
calculations in Chapter 7. Recall the form for the mean square radiated pressure from the 
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continuous spectral analysis approach in equation (7.18). In this case, the element physical 
coordinates are denoted by (£.//) and the non-dimensional element-centered coordinates will be 
denoted by (//,/'). Equation (7.18) can be written in terms of the element physical coordinates 

as follows 




V y 

V V / 1 1 1 {£Wf.,»f J o'IJSWl'*'” " ■ ' ~- <l 

{llv ({'■ ’/)} , (111 (!{' fh l 


I) a b a 




* H */ 

»=' J - 1 0 () (I 0 


(B.8) 


Applying the variable transformation given in equation (B.2) yields 


1111 


?uu.o ) . (= yttjjjj 




(L 


'=* 7— 1 — i -l-i-i 


{//„,(,/.//)) |./(//.;/)||./(//V)| dfi <\v <h> ( h ' 

1 ' (B.9) 


where the |./| are determinants of the Jacobian of the transformation as before. Now, applying 
Gauss-Legendre numerical integration gives the final expression for the mean square pressure at 
a point (i.e., for a given set of observer direction angles 0 and ©) in the acoustic far-field: 


//-'(/?. 8 . 0 ) = 


/-l j — 1 A = 1 l=\ m=l n = l 
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■ ^ 


( V rn ~~ 


| fl w (//„,. ) | | Jil'k- u l ) 1 • ,7 II ) I 

(B .10) 


The same procedure can be applied to the expression for the mean square radiated pressure from 
the discrete spectral approach, equation (7.35), which yields 


A ^ A <-■ ^ g ^ 9 ^ 9 ^ i 


T 


t= ] J= i A = i /=l m=l n = l 

~ 1 —f i • (»>„-«).■ ) K + ( ‘-"“‘'I ) A r . 1 , , , 

r=0 




(B.l 1) 


It is apparent that numerical integration for the computation of the radiated pressure can be 
computationally intensive due to the multiply-nested loops. Clever computer coding and/or 
compilation optimizers can diminish this problem, but another approach would be to evaluate 
the necessary element area integrals in closed form. This can save significant computational 
effort and a discussion of this approach is given in Appendix C. 
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B.3 Radiated Acoustic Power 


Determination of the radiated acoustic power requires integration of the normal acoustic 
intensity over a hemispherical surface, given by equation (E.6). This integral must be evaluated 
numerically. Gauss-Legendre numerical integration can be employed through selection of a 
suitable change of variables. The following change of variables results in the desired integral 
limits: 


0:[O.ar/2] u 


K'-t) 


[- 1 . 1 ] 


<10 = - 


<ht 


(B.12) 


o : [0. 27r] v = — (o - 7r) : [— 1. 1] do = it dr 

7i 


Then the radiated acoustic power becomes 


i j 


11= I I / (r. u, r) /•* sin r)| r/y/ 


dr 


-I -1 


(B.13) 


where the determinant of the Jacobian matrix is as follows 


(B.14) 

as is easily seen in making the substitution for <10 and do in equation (B.13) from the transfor- 
mation. Applying Gauss-Legendre integration results in the following expression for the total 
acoustic power radiated by the panel 



00 


7V 
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dii 
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do 

Ou 

00 

Ov 
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7T 


A', Ny 

n = 5^53 /(r, r- sin (ji/,- + j)|.7(« r . r,)| (B.15) 

r=] s=l 

where the acoustic intensity / is related to the mean square radiated acoustic pressure by equation 
(E.4) and the mean square pressure is given in equation (B.10) or (B. 11). 
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Appendix C 

Closed Form Evaluation of Element Area Integrals 


Recall the expressions for the far-field mean square radiated acoustic pressure, given in 
equations (7.18) or (7.35): 


A'. *V< r r ... f (<'-«)k+(>i'-»i)a 

i?(».«)=(4)'EE// <*'«•*»/ .1 7 * 

i = l .;= 1 y — oc 


(C.l 


or 


A, A\ 


p F '(/(.^.o) = (^)'EE / / 

,=1 '=>.1 r. 


A. -I 




^V- U "° K+( — {^U / -7 / )} j dAidA'i 


n = » 


(C.2) 


The form for the out-of-plane interpolation functions {II „•{(,■ >])} f° r the 24 degree-of- freedom 
rectangular plate element are given in equation (4.4). As stated previously in the respective 
acoustic radiation sections, the complex exponential function in the integrand of the above 
equations must be separated into two components. One is constant for each element index pair 
(i.j ) and the other must be included in the integration over the element areas in local coordinates. 

Consider the expression for the mean-square pressure in the discrete frequency domain 
formulation for example. Splitting the complex exponential into element-constant and element- 
variable components and changing the order of the summations gives 


a) = (^) £ it, £ f '**” 

“ 3 w = ( ) 7=1 j = 1 

f {H w (u.v)}‘e- l ^^dA t 

A, 


(<;'-£■ ) K +('/, ) A 


r : . uK+ 1 'A 
i — 7 o/’ 7J 


(C.3) 




where (a, r) are element local coordinates. There are four cases which must be considered 
in evaluation of these area integrals; (I\ ^ 0. A ^ 0). (k ^ 0. A = 0), (k = 0. A ^ 0) and 
(K = 0. A = 0). Integration of each of the 16 out-of-plane interpolation functions with the 
complex exponential, for each of the four cases involving the argument of the exponential 
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function, was performed. The resulting functions are denoted as follows 
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(C.4) 


where a and b are the element length and width, respectively. The expressions resulting from 
integration of the shape functions along with the various exponentials follow. Note that wave 
number is used in these expressions^- = u ;fc. 

(K^O.A/O) 


/,(*. K. A) = -[(-12/ + 12 h MA + m\ + m\e MA + bH- :l A :i ) 

(-12/ + 12/Y wkK + 6//AK + Gf/A-Kf " ,AK + fl 3 A- 3 K 3 )]/(« ;, 6 3 K 4 A 4 A- 8 ) 


(C.5) 
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7j(A\ K. A) = [( 

— 127 + l2it ,bkA + 6 ft A- A + GftA-A/ MA + ft 3 A- 3 A 3 ) 

( 

— 12/ + 127/ '" AK + 6a A K + 6aAkz'" AK + a 3 A- 3 K V' tAK )]/(a 

7;t( A*, K, A) = — | 

(— 127 + 127 / lhkA + GftA-A + GftA-A/ ,hkA + ft 3 A- 3 A 3 / ,MA ) 
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/ U (A\ K. A) = [(f> - Gz ibkx + iibk A + 2/ftA-Az MA - lrk-\ A 

(g - (>< " ,AK + 2/V/AK + 4/V/A-Kz " ,AK + z/ 2 A- 2 K 2 / i " AK )]/ (// 2 ft 2 I\ 4 A 4 / 

/ A-. K. A) = [((i- Gf *'“' A 4- 2/M- A 4- 4/ftA-A/ ,MA + ft 2 A- 2 A 2 / M A ) 

(fi - Gz"* AK + 2/V/AK + 4/V/A-Kz " ,AK + «-k 2 K-< '" AK )]/(// 2 ft 2 K 4 A 4 / 

/ , fi ( A- . K. A) = [(g - Gz' 7,AA + 2/ftAA + 4/ftA-\z' 7 ' AA + trlr.\ 'c MA ) j 
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/ 8 (A-. K) = //(g - Gz"* ak + 4///A-K + 2/V/A-K / i<,AK - // 2 A- 2 K J ) / ( 2// - A- 4 K 4 ) 
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(K = 0. A / 0) 


/, (A-. A ) = «(l2 - 12e ihkA + (nbkA + bihkAt lbkA + /7> ! *- 3 A' 5 ) / (2 b 4 k 4 A 1 ) (C.37) 


A ) = a ( 1 2 - 1 2t lbkA + biblxA + dibkAt MA + it 4 k 3 A 3 ) / (2 frV A 1 ) (C.38) 

/:5( A, A) = «(— 12 + 1 2( - dibkA - dibkAt ibkA - ib 3 k 3 A 3 t MA )l(2b 3 k 4 A 4 ) (C.39) 

A) = «(-12 + \2c lbkA - dibkA - dibkAt ibkA - rl>H i A i e lhkA )i(2b 3 k 4 A 4 ) (C.40) 

A) -- „ J (l2 - 12< * 6AA + (nbkA + dibkAt ibkA + ib*k\\ : , )/(12ft 3 A- 4 A 4 ) (C.41) 

\ ) = </' (— 12 + 1 2e ibkA CubkA -dibk\ ( ,bkA - //>*/.• l A :t ) / ( 1 2/# V. 1 A 1 ) (C.42) 

/t( A, A) = «-(l2 - l'2< ibkX + dibkA + dibkAc MA + iti'k\\*< ibkX )/(l2lf i k 4 A 4 ) (C.43) 
/ 8 (A, A) = a“ (—12 + I2t ibkA - dibkA - dibkAt tbkA - ib 3 k 3 A 3 i MA ) / (\2b 3 k 4 A 4 ) (CM) 

J 9 (k.A) = a(d - d<‘ bkA + AibkA + 2ibkA< MA - Irk 2 A 2 ) / (2lrk 4 A 4 ) (C.45) 

7 U) ( A-. A ) = o (Vi - 6 < MA + UbkA + 2ibkAt MA - lrk 2 A 2 ^ / (2 tr k 4 A 4 ) (C.46) 

7„(A-. A) = a(d - (it thkA + 2 ihkA + AibkAt MA + lrk 2 A 2 t lbkA ^/{2lrk 4 A 4 ) (CM) 

/ 12 (A-, A) = a (d - 6 c MA + 2 ibkA + libk Ac ibkA + lrk 2 A 2 ( MA )/(2b 2 k 4 A 4 ) (CAS) 

/ 13 (A-. A) = « 2 (6 - dt ,UA + AibkA + 2ibkAt lhkA - lr k 2 A 2 } / {\2tr k 4 A 4 ) (C.49) 

7„( A-, A) = a 2 (— G + 6t ibkA - AibkA - 2ibkAt >bkA + &W)/(12&W) (C.50) 

7 15 (A-. A) = «* (—6 + de ibkA - 2ibkA - AibkAt ibkA - lrk 2 A 2 e' hkA )/(Ub 2 k 4 A 4 ) (C.51) 

7 16 (A\ A) = A." (g - d< MA + 2ibkA + AibkAc ibkA + Irk 2 A 2 e lbkA ) / (I2b 2 k 4 A 4 ) (C.52) 



(C.53) 

(C.54) 
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(C.55) 

a 2 lr 

1 14- ^16 

(C.56) 


These expressions were used as function calls to evaluate the vectors pre- and post-multiplying 
the element-ij displacement spectral matrix in equation (C.3). 
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Appendix D 

Radiation Far-Field Approximations 

Prediction of the radiated acoustic pressure to a free field from a vibrating surface requires 
evaluation of Rayleigh’s integral, equation (7.1). The present work casts the radiated acoustic 
pressure response in three forms: continuous spectral analysis of stochastic response, given by 
equation (7.10); and discrete spectral analysis of deterministic and stochastic responses, given 
by equations (7.24) and (7.31). The existence of the radial distance dependence in the integrand 
requires the use of numerical integration or some simplifying assumptions. The integrand can be 
simplified extensively by considering only far-field observer locations in the acoustic free-field. 
Let the observation point (x.i/.z) in the acoustic field be designated by its distance R from the 
panel center and two direction angles: the polar angle 0 = [0.7 r/2], and the azimuthal angle 
o — [0. 2 tt], Thus, 


n = \A- 1 2 + 7 2 + ~- 

,r = R sin 0 cos 0 
i) = H sin 0 sin o 
= /? cos 0 


(D.l) 


and the distance from a point on the panel to the observation point 
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The acoustic far-field approximation entails neglecting the terms in the expression for r that 
involve 1 / R 2 and letting 


1 




/• R 


The remaining terms in the expression for r may be expanded in a binomial series: 
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Retaining the first two terms in the expansion results in 
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where I\ = sin 0 cos d and A = sin0sin<£. Thus, the exponential term in the integrand of 
equation (7.24) can be simplified as follows. 

( -'-’-‘7 ~ , (D.8) 

Finally, the difference in radial distances encountered in equations (7.10) and (7.31) becomes 

Ktf'-0 + A(,/-,/) (P-9) 

and the corresponding exponential term in the integrand takes the following form. 

% (D.10) 
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Appendix E 

Acoustic Intensity, Power, and TL 

The acoustic intensity at a point is equal to the time average of the product of the perturbed 
pressure and velocity. For a pressure and velocity that are random processes, this temporal 
average must be taken in the limit as the characteristic time T — > oc: 

T 

l(.v, //.;) = (/)( .r,y. zj)) = lim Tp f p(.r.y, z.i)u{.r.y, =.i) di (E.l) 

/ — oc / J 

0 

where ii(.r. /) is the acoustic velocity vector. This temporal average may be replaced by 
an expected value: 

/(•' •'/•-) = {p[-v.y.:J)d(.r, y.z.l)) 

= E[p{ .v. y, z. t)ii(.r. //. c, /)] 

by assuming that the processes are ergodic in correlation. In general, the relationship between 
the acoustic pressure and velocity is not known and must be determined concurrently. However, 
in the acoustic far-held, the radial velocity is dominant and the propagating acoustic waves 
appear planar. Consequently, the relationship between the acoustic pressure and velocity can 
be approximated by 

p(H.e.o.t) 


(E.2) 


u(H.e.o.t) ~ ur(ru).Q.t) 


pc 


(E.3) 


where the spatial description has been changed to spherical coordinates for consistency with 
equations (7.18), (7.35), or (7.37). Thus, the far-held radiated acoustic intensity becomes 


I,{R.0.o) = —E\pj(n. O. 0 J) 1 = ^ 
pc pc 


(E.4) 


where the mean square acoustic pressure is given by one of the equations referenced above. 

The acoustic power incident upon a surface is dehned as the integral of the normal acoustic 
intensity over the surface: 


II 


-/ 


I{r. 0. 0 ) ■ h dS 


(E.5) 


where n is the surface unit normal. Therefore, the total acoustic power radiated by the plate is 
calculated by integrating the acoustic intensity over a far-held hemisphere as follows. 


x/2 


Mr 


// 

0=0 0=0 


I(r. 0. 6) H~ sin 0 dO do 


(E.6) 


where the radiated intensity has the form given in equation (E.4). Note that the H 1 in the 
integrand cancels with that in the expression for the far-held acoustic intensity, equations 
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(7.18)/(7.35)/(7.37) and (E.4). This integral must be evaluated numerically, which can be 
accomplished via Gauss-Legendre numerical integration (see Appendix B). 

Recall that for the case of acoustic pressure excitation the incident acoustic pressure consists 
of obliquely incident plane waves. The incident acoustic intensity is then, from the above 

discussion, __ 

li = li (E.7) 

f)C 

and the incident acoustic power is simply the product of the normal acoustic intensity and the 
area of the plate 


n . _ l>> a,>( ’ os Oi (E.8) 

pc 

Finally, the transmission loss (TL) is defined in units of dB by the logarithm of the ratio of 
incident acoustic power to radiated acoustic power 


TL 


= 10 log 



(E.9) 
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Appendix F 

Continuous Spectral Examples 


A continuous spectral analysis was one approach used in Chapter 6 to derive the expression 
for the dynamic response correlation matrices of a panel-type structure subjected to elevated 
and/or one of the following mechanical loads; acoustic pressure, concentrated forces, or base 
acceleration. A similar approach was used in Chapter 7 to derive an expression for the mean 
square radiated acoustic pressure due to the same structure, but limited to two of the three 
mechanical loads; acoustic pressure or concentrated forces. The dynamic response correlation 
matrices and mean square radiated pressure require evaluation of spectral integrals of the forms 
given in equations (6.36) and (7.15), respectively. It is clear that these equations have a very 
similar form and will be denoted by '"'(r) in the following generic expression 


* m a 


T) 


oc 



— oc 


(F.l) 


where the expressions for the frequency response functions //*(o.’) have been substituted from 
equation (6.29) and >'**(«;) is the power spectral density function of the excitation. The 
superscript (<7, r. a) indicates that the multiplier of corresponds to displacement, velocity, 

and acceleration response, respectively. 

Two sample cases for the evaluation of the integral in equation (F.l) are presented in this 
appendix including harmonic and band-unlimited white noise excitation. The case of harmonic 
excitation is obviously not stochastic, but is included as a simple comparative case. 


F.l Harmonic Excitation 


Harmonic excitation can be described in either of the following forms 


So 
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|u-j — u?() 

otherwise 


.5 ++ (u.’) — .So <S(|cc| 


*■’0 J 


(F.2) 


or 

f**(r) = 2 .S'o cos u?o t (F.3) 

where .% is a constant, <*>(|^j — u;o) is a spectral Dirac delta function, r = i> — t\, and the 
subscript ** indicates the excitation could be due to acoustic pressure, concentrated forces, or 
base acceleration. Substitution of 5'**(u,’) into equation (F.l) results in 
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This equation may be simplified by retaining only even terms, because the integral of odd terms 
over even limits vanishes, as 


llnn ' a \ T ) 



) + 1 Cut { n ^ w ^ ^ 
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Enforcing the Dirac delta function results in 
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(F.6) 

This expression can be used directly for calculation of the dynamic response correlation matrices 
given by equation (6.38). The mean square radiated pressure makes use of the acceleration 
response case, i.e., by making the substitution 

_ (s' ~ £)K + ('?' ~ ? /)-^ (F.7) 

c 

and using the result in equation (7.18). Finally, the mean square dynamic responses can be 
obtained from the diagonal elements of the dynamic response correlation matrices for zero time 
difference (r = 0) which get input from the following simplified form for I,„„ (0) 


tirHo) 


.yq .(0.2.4) 


[(u-7„ — uJo)(q-’n - ) + KmCn^-3)»^»^-n] 
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F.2 Band-Unlimited White Noise Excitation 

Band-unlimited white noise excitation may be described by either of the two following forms. 


= ri’ll 


(F.9) 


C**(t) = 27tSoHt) (F.10) 

where .% is a constant for all frequencies and r = /•> - 1\. Substitution of NtU’i i nto equation 
(F.l) results in 
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Figure F. 1 Schematic of upper-half-plane contour for evaluation of 
dynamic response due to band-unlimited white noise excitation. 


One way to evaluate the integral in equation F. 1 1 is to employ the Cauchy residue theorem by 
integrating over a semicircular contour, in the limit as the radius of the contour becomes infinite, 
in the upper half plane (for r > 0) enclosing the poles of the integrand in the UHP. The case 
for C(>), i.e., acceleration response — O in the numerator, in equation (F. 11) is unbounded. 
Therefore, attention is limited to C* *(r), the cases for O 0 -’' in the numerator. A schematic 
of the contour of integration is shown in Figure F.l. 

The integral around the closed contour can be written as 

fU')t U ' r (l u? 
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n 

= //MO^ + j f{uj)t^ T d u 


(F.l 2) 


-R 


Cu 


where the integrand in equation (F.l 1 ) has been denoted by and Cr denotes the circular 

portion of the contour. In the limit as R — > oc, the first integral on the right-hand side takes 
the desired form, the Cauchy principal value, and the sum of the two integrals can be evaluated 
by employing the Cauchy residue theorem as 
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where the expression on the right of the summation represents the residue of tu!T at the 
pole i'k and A r is the number of poles of /(<*?) enclosed within the contour (which becomes all 
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of those in the UHP as R — > oc). Recall the expression for j M shown in equation (F.l 1 ). This 
expression can be rewritten in terms of simple poles as follows 
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where 
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the residue of ./V) at the pole *>; is given by the coefficient of 1 --«’») 111 the Laurent 

expansion (see Churchill and Brown [61]) given by 
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Thus, the integral on the closed contour becomes 
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It can be shown that the integral on the circular portion of the contour is zero (see Hildebrand 
[62]) so the desired result has been achieved as 
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The terms in this expression can the expanded into a sum of real and imaginary components 
and simplified to the following 

{d) . , 47T(C„, +C»)*o c~ gn r c°s^/„r 


Qr) = 


(u-V/,,, — u/’d,, ) + (Cm + Cn) 


•'dm 


+ ^\ln )’ + (C m + Cm 1 


4i: 

-■<!„ -u 4m ~ ( V"' + ^ 

,s'„e ^” r sin 

Oi 

+ 

+ 

1 

{^dt) 

2 J 

t + ^\hi ) + (s»i “f V« ) 


(F.20) 


219 



and 




+ ClT C0Su,' (iii T 


(“'dm ^ \In ) (£m “I" £» ) 


(^/m dn ) (£m + £,,)■ 


1 7T 

u 'dn ^'dn^'din u ''dn*m -^'dn + “1" ) 

■Si< Z" T si]i^’ ( /„r 


(“'dm “'dn ) (£m “h )” 

(“'dm "t” “'dn ) (vm “1“ V// ) 



(F.2 1 ) 

These expressions can be used directly for calculation of the dynamic response correlation 
matrices given by equation (6.38). The mean square radiated pressure cannot be computed in 
this case as the acceleration response for band-unlimited white noise is unbounded, as stated 
previously. Finally, the mean square dynamic responses can be obtained from the diagonal 
elements of the dynamic response correlation matrices for zero time difference (r = 0) which 
get input from the following simplified forms for I If,!, ( 0 ) and 7,^(0) 


7m,! (0) = 




'dn 


“< Im 


(£ro + £« r 


( u 'dn H“ “'dm ) H“ £n 


(¥22) 


and 


7‘:;i(0) - 


4jt.S 


0 [vm (^T/„ + (.n ) + + Cm )] 


(“^dn dm ) ”t~ (£m "i" ) 


i u '\ln dm ) “1“ (<v/h H“ s n 


(¥23) 


220 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No, 0704-0188 


garnering an d maintaining me data needed, and compietmg and revising ln tormation Operations and Reports, 12tS Jetterson 

as ~ n ' ocgo - — 

un HATCC COVFRFD 


1 . AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

January 2001 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


4. TITLE AND SUBTITLE T1 , ., - 

Thermomechanical Response of Shape Memory Alloy Hybrid Composites 


5. FUNDING NUMBERS 

WU 706-32-31-02 


6. AUTHOR(S) 

Travis L. Turner 


7. PERFORMING 6toAti\Ztf\6ti NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-2199 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L- 18057 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA/TM-2001 -2 1 0656 


11 T^rSmatL^sented in this report was offered as a thesis in partial fulfillment of the requirements for the 
Degree of Doctofof Philosophy in Engineering Mechanics, Virginia Polytechnic Institute and State University, 
Blacksburg, Virginia, November 2000. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b ' DISTRIBUT,0N CODE 

Unclassified-Unlimited 

Subject Category 7 1 Distribution: Standard 

Availability: NASA CASI (301 ) 621-0390 

13. ABSTRACT (Maximum 200 words) 

This study examines the use of embedded shape memory alloy (SMA) actuators for adaptive control of the thermo- 
mechanical 6 re sponse of composite structures'. A nonlinear thermomechamcal model is presented for analyzing 
shape memory Moy hybrid composite (SMAHC) structures exposed to steady-state thermal and dynamic mechan- 
ical loads Also presented are ( 1 ) fabrication procedures for SMAHC specimens, U) characterization of the constit- 
uent materials for model quantification, (3) development of the test apparatus for vll.dat" 

experiments on specimens with and without SMA, (4) discussion of the experimental results and (5) validation of 
thl ScaUnd numerical tools developed in the study. Excellent agreement is achieved between the predicted 
and measured SAMHC responses including thermal buckling, thermal post-buckling and dynamic response due to 
inertial loading The validated model and thermomechanical analysis tools are used to demonstrate a variety of 
static and dynamic response behaviors including control of static (thermal buckling and post-buckling) and 
dynamic responses (vibration, sonic fatigue, and acoustic transmission), and SMAHC design cOTSideraUons for 
these applications. SMAHCs are shown to have significant advantages over conventional response abate 
approaches for vibration, sonic fatigue, and noise control. 


14 Nkinol; Embedded actuators; Hybrid composite fabrication; Nonlinear 


thermoelasticity; Geometric nonlinearity 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 

NSN 7540-01-280-5500 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

239 


16. PRICE CODE 

All 


20. LIMITATION 
OF ABSTRACT 

UL 

Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 




